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Abstract 

Flavor oscillations by itself and its coupling with chiral oscillations and/or spin- flipping are the 
most relevant quantum phenomena of neutrino physics. This report deals with the quantum theory 
of flavor oscillations in vacuum, extended to fermionic particles in the several subtle aspects of 
the first quantization and second quantization theories. At first, the basic controversies regarding 
quantum-mechanical derivations of the flavor conversion formulas are reviewed based on the inter- 
nal wave packet (IWP) framework. In this scenario, the use of the Dirac equation is required for a 
satisfactory evolution of fermionic mass-eigenstates since in the standard treatment of oscillations 
the mass-eigenstates are implicitly assumed to be scalars and, consequently, the spinorial form of 
neutrino wave functions is not included in the calculations. Within first quantized theories, besides 
fiavor oscillations, chiral oscillations automatically appear when we set the dynamic equations for 
a fermionic Dirac- iype particle. It is also observed that there is no constraint between chiral os- 
cillations, when it takes place in vacuum, and the process of spin-flipping related to the helicity 
quantum number, which does not take place in vacuum. The left-handed chiral nature of created 
and detected neutrinos can be implemented in the first quantized Dirac theory in presence of mix- 
ing; the probability loss due to the changing of initially left-handed neutrinos to the undetected 
right-handed neutrinos can be obtained in analytic form. These modifications introduce correction 
factors proportional to m^/El that are very difficult to be quantified by the current phenomeno- 
logical analysis. All these effects can also be identified when the non-minimal coupling with an 
external (electro)magnetic field in the neutrino interacting Lagrangian is taken into account. In the 
context of a causal relativistic theory of a free particle, one of the two effects should be present in 
flavor oscillations: (a) rapid oscillations or (b) initial flavor violation. Concerning second quantized 
approaches, a simple second quantized treatment exhibits a tiny but inevitable initial flavor vio- 
lation without the possibility of rapid oscillations. Such effect is a consequence of an intrinsically 
indefinite but approximately well defined neutrino fiavor. Within a realistic calculation in pion 
decay, including the quantum field treatment of the creation process with finite decay width, it is 
possible to quantify such violation. The violation effects are shown to be much larger than loop 
induced lepton flavor violation processes, already present in the standard model in the presence of 
massive neutrinos with mixing. For the implicitly assumed fermionic nature of the Dirac theory, 
the conclusions of this report lead to lessons concerning flavor mixing, chiral oscillations, interfer- 
ence between positive and negative frequency components of Dirac equation solutions, and the fleld 
formulation of quantum oscillations. 
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I. INTRODUCTION 



Particle mixing and flavor oscillations 



138 



139] continue to stimulate interesting and 



sometimes fascinating discussions on the many subtleties of quantum mechanics involved in 
oscillation phenomena 4-6|. The flavor mixing models the quantum field prescrip tions 



13 



140| and, generically, the quantum mechanics of oscillation phenomena 



14 



19 



142 1 



lave been extensively studied in the last years. In particular, the properties of neutrinos 



20 



2l| obtained in all of these frameworks have become the subject of an increasing num- 



ber of theoretical constructions. Notwithstanding the exceptional ferment in this field, the 
numerous conceptual difficulties in describing accurately the particle mixing and oscillations 
have renewed the interest in understanding the derivation of the flavor conversion proba- 
bility formulas and in overcoming the main physical inconsistencies hidden in the standard 
theoretical approaches. 

The flavor oscillation analysis have been supported by compelling experimental evidences 
which have continuously ratified that neutrinos undergo flavor oscillations in vacuum and in 
matter. One can focus, for instance, on the outstanding results of the Super-Kamiokande 

ligh- 



atmospheric neutrino experiment 22|], in which a significant up-down asymmetry of the 
energy muon events was observed, the results of the SNO solar neutrino experiment 123 



2^, 



in which a direct evidence for the transition of the solar electron neutrinos into other flavors 
was obtained, and also the results of the KamLAND experiment 25j that confirmed that the 
disappearance of solar electron neutrinos is mainly due to oscillations among active neutrinos 
and not due to other types of neutrino conversion mechanisms 26|, |27|. The experimental 



data could be completely interpreted and understood in terms of three neutrino flavors, with 
the exception of the LSND anomaly 
the MiniBoone experiment 



30 



24 



28 



29] . Such anomaly, although not confirmed by 



light neutrino flavor which had to be inert; T 

is largely excluded by the oscillation data 



3l|, led to speculations of the existence of (at least) a fourth 
le presence of such light sterile neutrinos 



32 



33] but depending on their mass scale it 



may influence certain astrophysical and cosmological phenomena ranging from the thermal 



evolution of the Universe 
part of dark matter 



34l | to supernova explosion, pulsar kicks and even a significant 



35|. On the other hand, the hypothesis of mixing between known 



neutrino species (electron, muon and tau) and higher mass neutrinos has much stronger 
theoretical motivation since it may account for the lightness of the active neutrinos through 
the seesaw mechanism 



36 



37] and also account for the matter and antimatter asymmetry 



of the Universe through the mechanism of leptogenesis 38|]. The observation of such mixing 



5 



at low energies, however, could be extremely difficult. 

In parallel, the neutrino spin-flipping attributed to some dynamic external [39] interacting 
process, which comes from the non-minimal coupling of a magnetic moment with an external 
electromagnetic fleld j40|, was formerly supposed to be a relevant effect in the context of the 
solar-neutrino puzzle. As a consequence of a non-vanishing magnetic moment interacting 
with an external electromagnetic fleld, left-handed neutrinos could change their helicity (to 
right-handed) 4l|. The effects on flavor oscillations due to external magnetic interactions 
in a kind of chirality-preserving phenomenon were also studied j42| but they lack a full 
detailed theoretical analysis. It was partially provided by some recent theoretical studies 
where the Dirac/Majorana characteristic of neutrinos becomes relevant 0, Q^iSl. Only for 
ultra-relativistic (UR) neutrinos, however, changing helicity approximately means changing 
chirality. One of our goals is to demonstrate that, in the context of oscillation phenomena 
and in the framework of a flrst quantized theory, the small differences between the concepts of 



chirality and helicity, which 



for massless particles 



39 



lad been interpreted as representing the same physical quantities 



42 



44 



46 



47| . can be quantifled for massive particles. It raises the 



possibility that chirality coupled to flavor oscillations could lead to some small modiflcations 
to the standard flavor conversion formula fl9|. 

On the theoretical background, under the point of view of a flrst quantized theory, the 
treatment of the flavor oscillation phenomena in terms of the intermediate wave packet 
(IWP) approach 

plane- wave formalism |49l . |50|. However, a common argument against the IWPs is that 
oscillating neutrinos are neither prepared nor observed js^. This point was partially clarifled 
by Giunti ll|] who suggested a solution in terms of an improved version of the IWP model 
where the wave packet (WP) of the oscillating particle is explicitly compute d wi th fleld- 



48] eliminates the most controversial points arising with the standard 



14Q|. 



theoretical methods by means of the external wave packet (EWP) approach [8|, 

Such approach, in contrast to IWP approaches, became the customary way to avoid the 
ambiguities involving the question on how neutrinos are created and detected. According to 
Ref. [8:l|, the IWP treatments are the simpler flrst quantized ones treating the propagation 
of neutrinos as free localized wave packets. In contrast, EWP approaches consider localized 
wave packets for the sources and detection particles while the neutrinos were considered 
intermediate virtual particles. 

On the other hand, another classiflcation scheme can be used do classify the various 
existing treatments considering a more physical criterion irrespective of the use of wave 
packets. It refers to the descriptions of neutrino oscillations that (A) include explicitly 
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the interactions responsible for the mixing and those (B) that only treat the propagation of 
neutrinos, i. e., the mixing is an ad hoc ingredient. A more subtle aspect in between would be 
the (explicit or phenomenologically modeled) consideration of the production (and detection) 
process(es). In general, the IWP approaches are of type (B). The EWP approaches are 
of type (A). The Blasone and Vitiello approach jol, although in the quantum field theory 
(QFT), is of type (B) since mixing is introduced without explicitly including the interaction 
responsible for it. The type (B) approaches have the virtue that they can be formulated 
in a way in which total oscillation probability in time is always conserved and normalized 
to unity jo], [5]]-54|. This feature will be present in all first quantized approaches treated 
here (sections [TTl and IIIip and in a second quantized version (section IVAI) . If different 
observables are considered, or a modeling of the details of the production and detection 



processes is attempted, further normalization is necessary jl4j . Il8l . 155] . In such cases, the 
oscillating observable might differ from the oscillation probability. On the other hand, 
type (A) approaches tend to be more realistic and can account for the production and 
detection processes giving experimentally observable oscillation probabilities [56]. Of course, 



57- 



they are essential to the investigation of how neutrinos are produced and detected 
59|. To rigorously derive a flavor /chiral conversion formula for fermionic particles (non- 
minimally coupled to an external magnetic field), we avoid the field theoretical methods in 
the preliminary investigation. As we are initially interested in the Dirac equation properties, 
as a first analysis, the IWP framework is a suitable simplification for the understanding of 
the physical aspects concerning the oscillation phenomena. 

In section |TT] of this manuscript, the spatial localization is included in the formulation of 
flavor oscillations. Quite generally, the analytical description of the dynamical evolution of 
a mass-eigenstate do not involve the wave packet limitations. In particular, the analytical 
properties of gaussian distributions 0.Q enable.. tKefl..^dtKe.econd 
order corrections to the oscillating behavior of propagating particles. We assume sharply 
peaked momentum distributions and then we approximate the mass-eigenstate energy in 
order to analytically obtain the expressions for the wave packet time evolution and for the 
flavor oscillation probability. In the IWP approach a consistent energy expansion is taken up 
to the second order term in the wave packet parameter (7^ ~ {aEi)~^, that satisfies a, ^ 1 
for sharply peaked momentum distributions. The wave packet spreading as well as the loss 
of coherence between the propagating wave packets are quantified in both non-relativistic 
(NR) and ultra-relativistic (UR) propagation regimes. Thus, the preliminary step of our 
study consists in a self-consistent approximation to the mass-eigenstate energy in order to 
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analytically obtain the expressions for the wave packet time evolution and for the flavor 
oscillation probability. We also identify an additional time-dependent phase which changes 
the standard oscillating character of the flavor conversion formula. 

The extension to the Dirac theory is introduced in section llllt where spin and relativistic 
completeness are included into the flavor oscillation formulation. This section is concerned 
with the analytical derivation of a flavor conversion formula where the fermionic instead of 
the scalar character of a propagating mass-eigenstate is assumed. To that end we shall use 
the Dirac equation as the evolution equation for the mass-eigenstates and show that the Dirac 
formalism is useful and essential in keeping clear many of the conceptual aspects of quantum 
oscillation phenomena that naturally arise in a relativistic spin one-half particle theory. More 
particularly, we show that a superposition of both positive and negative frequency solutions 
of the Dirac equation is often a necessary condition to correctly describe the time evolution of 
the mass-eigenstate wave packets. We give, for strictly peaked momentum distributions and 
UR particles, an analytic expression for the Dirac flavor conversion probability. A modified 
formula for the conversion probability is shown and an additional rapid oscillation term, 
coming from the interference between the positive and negative frequency contributions, is 
found. To completely disentangle the influence of the initial wave packet to the phenomenon, 
an analysis independent of the initial wave packet is conducted through the calculation of the 
Dirac time evolution kernel in the presence of flavor mixing. The properties of completeness 
and causality are briefly analyzed. An analogous calculation is performed for relativistic 
spin zero particles to show that rapid oscillations are indeed a consequence of the presence 
of positive and negative frequency solutions (completeness) for relativistic wave equations 
and not of the spin degree of freedom itself. Within first quantized Dirac theory, we also 
establish the inextricable relationship between two phenomena: initial flavor violation and 
rapid oscillations. 

Further consequences of spin structure and relativistic completeness, such as chiral os- 
cillations, are discussed in section IIVI Chiral oscillations naturally enter the discussion of 
neutrino oscillations because neutrinos are produced and detected through weak interac- 
tions that are chiral in nature, more specifically, left-handed in chirality. It is shown that 
the inclusion of chiral oscillation effects, together with the time-evolution of spinorial wave 
packets for the mass-eigenstates, can modify the flavor conversion probability formula. In 
particular, the probability loss due to the conversion of left-handed to right-handed neutrinos 
is calculated. The differences between the dynamics of chirality and helicity for a neutrino 
non-minimally coupled to an external (electro)magnetic fleld are expressed in terms of the 
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equation of the motion of the correspondent operators 7^ and h, respectively. In particu- 
lar, the oscillating effects can be explained as an implication of the zitterbewegung (ZBW) 
phenomenon that emerges when the Dirac equation solution is used for describing the time 
evolution of a wave packet [52] . Due to this tenuous relation between ZBW and chiral oscil- 
lations, the question to be analyzed in this section concerns with the immediate description 
of chiral oscillations in terms of the ZBW motion, i.e., we shall demonstrate that, in fact, 
chiral oscillations are coupled with the ZBW motion so that they cannot exist independently 
of each other. It provides the interpretation of chiral oscillations as very rapid oscillations 
in position along the direction of motion, i.e., longitudinal to the momentum of the parti- 
cle. In a subsequent step, we report about a further class of static properties of neutrinos, 
namely, the (electro) magnetic moment associated to the Lagrangian with non-minimal cou- 
pling. It allows the comparison between the dynamics of chiral oscillations and the dynamics 
of spin-flipping in the presence of an external magnetic field. It is also verified how the in- 
teraction with an external (electro)magnetic field can modify the neutrino fiavor oscillation 



formula 



53| . To summarize, the basic idea is thus to quantify the modifications that appear 



in the fiavor- chirality conversion formula, previously obtained for free propagating particles 
in vacuum |l9|, when an external magnetic field can affect chiral oscillations. 



In section |V] we finally analyze the inclusion of some aspects of field quantization into 
the description. Firstly, a simple second quantized description of the fiavor oscillation phe- 
nomenon is devised based on free second quantized Dirac theory. There is no interference 
term between positive and negative components, but it still gives simple normalized oscil- 
lation probabilities. We also review the central issue distinguishing the general IWP and 



EWP^ 



6 Of approaches: despite its direct unobservability, is the intermediate neutrino a 



real (on-shell) particle propagating freely? The answer is affirmative except for the possi- 
bility of contribution of the antineutrino component in the propagation of virtual neutrinos 



which, nevertheless, is very small js^ . Thus IWP description is a good approximation of 
the oscillation phenomenon 0, |6l|. We also compare the distinct approach of Blasone and 
Vitiello 9I with the first quantized description of neutrino oscillations. 

We have also considered that a central issue of the phenomenon of fiavor oscillations 
discussed along our manuscript is: how the coherent superposition of mass-eigenstate neu- 
trinos, i.e., the fiavor state, is created and detected 58 



59l . |62|? To answer such question, 
it is necessary to explicitly consider the interactions responsible for creation and/or detec- 
tion. Our contribution to such a fruitful discussion is a detailed calculation of the creation 
probability of the neutrino produced through pion decay performed using the full (pertur- 



9 



bative) QFT formalism at tree level, with explicit inclusion of pion localization. As a result, 
it is possible to study how the localization properties of neutrinos follows from the parent 
particle (pion) that decays. The calculation then provides the missing ingredients to quan- 
tify the new effect of intrinsic neutrino flavor violation 63|. Such effect is already present 
in first quantized formulations and it is manifested as an initial flavor violation or flavor 
indefinition but its presence was not mandatory and its magnitude could not be calculated a 
priori jsil- Moreover, we can show that the coherent creation of neutrino flavor states follows 
from the common negligible contribution of neutrino masses to their creation probabilities. 
At the same time, in the strict sense, we can also conclude that neutrino flavor is only an 
approximately well defined concept. 

The manuscript is structured in order to allow the reader to recognize the origin of each 
novel ingredient that can be included in the description of the quantum flavor oscillation 
phenomenon for neutrinos. We draw our conclusions in section | 



II. INCLUSION OF SPATIAL LOCALIZATION 



This section deals with the quantum-mechanical derivation of the oscillation formula. 



as well as the inclusion of spatial localization t 
extensively discussed in the literature m 



irough the IWP framework, as it has been 



50(1 . Our main contribution concerns the 



identification, through analytical expressions, of secondary effects coupled with the wave 
packet decoherence, which introduce small modifications to the oscillation pattern. In par- 
ticular, a preliminary discussion introducing the possibility of initial flavor violation and 
the respective consequences at neutrino creation/propagation/detection, is contextualized 
in the IWP framework. 

Associating a plane wave with each mass-eigenstate 



49 



50] is certainly the simplest and 



probably the most intuitive way to describe the interference phenomenon that gives rise 
to flavor oscillations in terms of an oscillation length and an oscillation probability. For 
oscillations between two different neutrino flavors that we will choose to be z/g and z/^ by 
convenience, the probability for flavor transition is usually expressed in terms of the mixing 
angle 6 and of the relative phase A$ by 



Vive i^fi) = sin'(2(9) sin' 
where the Lorentz invariant phase difference, 

A$ = A{ET -pL), 



(1) 



(2) 
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sets that initial pure flavor states are modified with time and distance. The mass-eigenstate 
phase difference A$ is then conventionally evaluated by setting AT = AL = and consid- 
ering, for ultra-relativistic (UR) particles, T ^ L and j ~ -E'1,2; i- 6. 



A<l> = TAE - LApRi L {AE - Ap) 
One thus gets the well-known expression jsQl 

V{Vf, ^ Vf^;L) = sin^(26') sin- 



Am^ 
2p 



L. 



Am^ 



Ap 



L 



(3) 



(4) 



Considering the plane wave approach, controversial points arises even in the derivation 



of formulas containing extra factors in the oscillation length |64j-|66|. In particular, the ex- 
tra factor of two in particle oscillation phases was discussed and re futed in the context of 
theory and phenomenology of neutral meson-antimeson oscillations 14l| . In addition, the 
simplified view of using wave packets (WPs) allows us to understand the origin of these 
extra factors. It is implicitly assumed that at creation the flavor state is unique even up to 
the phase at all points and times of creation. In the wave packet treatment, at time T and 
at a fixed position in the overlapping region, one experiences the interference between space 
points whose separation at creation is given by AvT and this implies that an additional 



initial phase is automatically included in the wave packet formalism |18|. l5l|. The final result 
contains the difference of phase given in Eq. ([3]). We do not intend here to re-discuss the 
many controversies in the plane wave derivations of the oscillation probability formula. We 
only remark that an approach strictly considering plane waves leads to conceptual difficul- 
ties and fails to explain fundamental aspects of particle oscillations, such as localization and 
coherence length. Wave packets eliminate some of these problems j48]. In fact, the use 
of wave packets for propagating mass-eigenstates (IWP model) guarantees the existence of 
a coherence length, avoids the ambiguous approximations in the plane wave derivation of 
the phase difference and, under particular conditions of minimal loss of coherence, recovers 
the oscillation probability given in Eq. (jl]). Moreover, the coherence necessary for neutrino 
oscillations depends crucially on localization aspects of the particles involved in the produc- 
tion of neutrinos l^^ . This point of view can be supported by quantum field theory (QFT) 



arguments as well 



11 



In practice, the loss of coherence is only relevant for neutrinos traveling cosmological 



distances 



67] . At the same time, it is not easy to determine the size of the wave packets 



at creation and it is not clear whether it makes sense to consider a unique time of creation 
62| . It configures a common argument against the IWP formalism, i.e., oscillating 
neutrinos are neither prepared nor observed. Consequently, it would be more convenient to 



11 



write a transition probability between the observable particles involved in the production and 
detection process. This point of view characterizes the so-called EWP approach js], [ll|. The 
oscillating particle, described as an internal line of a Feynman diagram by a relativistic mixed 
scalar propagator, propagates between the source and target (external) particles represented 
by wave packets. The function which represents the overlap of the incoming and outgoing 
wave packets in the EWP model corresponds to the wave function of the propagating mass- 
eigenstate in the IWP formalism. Remarkably, it could be shown that the probability 
densities for UR stable oscillating particles in both frameworks are mathematically equivalent 







8|. However, the IWP picture brings up a problem, as the overlap function takes into 
account not only the properties of the source, but also of the detector. This is unusual for 
a wave packet interpretation and not satisfying for causality [s^. This point was clarified by 
Giunti who evaluates the problem by proposing an improved version of the IWP model 
where the wave packet of the oscillating particle is explicitly computed with field-theoretical 
methods in terms of external wave packets. Despite of not being applied in a completely free 
way, the (intermediate) wave packet treatment commonly simplifies the discussion of some 



physical aspects associated to the oscillation phenomena |5lL |66j. Thus, it makes sense, 
as a preliminary investigation, to consider a wave packet associated with the propagating 
particle. 

A. The IWP framework 

The main aspects of oscillation phenomena can be understood by studying the two flavor 
problem. In this case, by associating the wave packets 0^ and 02 to mass-eigenstates and 
i^a, flavor wave packets can be described by the Vf,-\\ke state vector 

^{x,t) = 01 (a;, i) 0036*1/1 + 02 (a;, t) sin ^1^2 

= [0i(a;, i) cos^ ^ + 02(a;, t) sin^ ^] Uf, + [4>2{x,t) — (f)-i{x,t)\cos9 s\n6 

= (t)^^{x,t;0)Ue + (l)u^{x,t;e)Vfj,, (5) 

where Ve and i/^ are flavor states that are related to mass-eigenstates and 1/2 by the 
mixing relation 

I u„\ I rns sin ^ \ I u-i \ 

(6) 







f cos 6 


sine 1 


II 








[ — sin 9 


cos 9 , 







The mixing relation ([6]) can be also written 

i^a = Uaii^i, a = e,/Lt, i = l,2, (7) 
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where the mixing matrix U can be easily extracted. 

It is important do emphasize that (f)i{x,t) and (j)2{x,t) are usual normalizable wave 
functions, normalized to unity, while (j)i,^{x,t;6) and (f)^^{x,t;6) are not normalizable wave 
functions because of their time dependent norms (flavor oscillation). What is normalizable, 
as we will see, is the total probability over all flavors. 



(fx 



which is automatic, with U in Eq. (j7j) being unitary, once 0i and 02 are normalized to unity. 



d^x \4>i{x,t) 



1 



d^x\(t)2{x,t)('' 



We should remark that, in general IWP treatments, automatic normalization of total prob- 
ability is not guaranteed if the detection process is modeled by detection wave packets as 
in Refs. fl^. flS], i.e., a further normalization procedure is necessary. Therefore, within IWP 
treatments, we will not model the detection process and only consider idealized measure- 
ments that are consistent with the above normalization conditions. 

In addition to the restriction to two families, substantial mathematical simplifications 
result from the assumption that the space dependence of wave functions is one-dimensional 
{{x,t) — )■ {z,t)). Therefore, we shall use these simplifications to calculate the oscillation 
probabilities. The complementary effects of a 3-dimensional analysis are well explored in 
Ref. jsj. The probability of finding a flavor state j/^ at the instant t is equal to the integrated 
squared modulus of the i^^ coefficient 

sin^ (26) 



dZ \4)y^ 



{ 1 - INT(0 } , 



where iNT(t) represents the mass-eigenstate interference term given by 



iNT(t) = Re 



dz (p[{z,t) (t)2{z,t) 



For mass-eigenstate wave packets given by 

r ^21 

exp [ip, z], 

at time t = 0, where i = 1, 2, the corresponding time evolution is given by 

dp^ 



( 2 V 













27r 



f{p.-P^) exp 



-iE^]t + ip,z 



(8) 
(9) 
(10) 
(11) 



where Ep} = {pi + mf)^ and (f{p^-Pi) = (27ra^)4 exp 



. To obtain the oscillation 



probability, we can calculate the interference term iNT(i) by integrating 



^{Pz-Pi)^{Pz-P2) exp [-i AEp^ t] 



exp 



-{aApf 



dpz 
27r 



^^iPz-Po)exp [-lAEp^ t], 



(12) 
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where we have changed the 2;-integration into a p^-integration and introduced the quantities 
Ap = Pi~P2, Po = + P2) and AEp^ = Ep]^-Ep^]. The oscillation term is bounded by 
the exponential function exp[— (a Ap)^/8] at any instant of time. Under this condition we 
would never observe a pure flavor state. Moreover, oscillations are considerably suppressed 
if aAp > 1. Hence, a necessary condition to observe oscillations is that aAp ^ 1. This 
constraint can also be expressed by 6p ^ Ap where 6p is the momentum uncertainty of 
the particle. The overlap between the momentum distributions is indeed relevant only 
for 6p ^ Ap. Strictly speaking, we are assuming that the oscillation length {t^ ^^^2 ) is 



sufficiently larger than the wave packet width. It simply says that the wave packet must not 
extend as wide as the oscillation length, otherwise the oscillations are washed out 
Turning back to Eq. (fT2l) . without loss of generality, we can assume 



48 



61 



iNT(t) = Re <( / ^ ^\v.-V,) exp [-i AEp^ t] | . (13) 

This equation is often obtained by assuming two mass-eigenstate wave packets described 
by the same momentum distribution centered around the average momentum p = Pf,. This 
hypothesis also guarantees instantaneous creation of a pure flavor state i^e at t = [51]. In 
fact, for 0i(2,O) = 02(^,0), we get, from Eq. ([5]), 



./)^,(z,0,^) = (/.i(z,0) =(/)2(^,0) = ( ^Vexp 



exp [ipo z\ (14) 



^vra 

and (z, 0,0) = 0. Therefore, in what follows, we shall use this simplification. 



1. Spatial localization versus temporal average 

The fiavor conversion probability from fiavor Ue to is basically the squared modulus 
of the probability amplitude (f),y^{x,t) of the state in Eq. ([5]), with its initial localization 
properties determined by (j)^^{x,0), integrated over all space. The probability thus depends 
on time. But in a real experiment, time t is not a measurable variable; just the distance L 
between the source and the detector is known. To rewrite P{t) as P{L), the formula L = vt 
describing the trajectory of a free classical particle is, sometimes, inadvertently invoked. 
To clarify this point, let us calculate the probability of the beam of particles, produced at 
|a;| =0, to reach a physical detector of volume V, at average distance |a3| = L, by integrating 
the corresponding current density of probability }{x,t) over the surface dV enclosing the 
detector and integrating over the time of observation from ti to t2, as 

V{t,<t<t2) = - [ dt [ dS-j{x,t), (15) 

Jt^ JdV 

14 



where the minus sign arises because we want to quantify the flux entering V but dS points 
outwards the surface dV. This procedure, although straightforward and natural, is not 
generally adopted and more complicated methods are used instead. The reason for the 
rejection of Eq. (fT5|) is very simple: there is a difficulty in defining correctly the probability 
current density }e,fi{x,t), for fiavor defined neutrino states Ug or z/^. Indeed, such neutrino 
states have undefined masses, and in general it is not possible to define conserved currents for 
themy|. Nevertheless, it is possible to define an approximately conserved current je,^(a;,t) 
that obeys 



where we simplified the notation by using 0e^^(a;,t) = (p^^ i,^{x,t). The terms violating this 
conservation are proportional to the neutrino mass differences and can be neglected locally. 
The explicit construction for wave functions satisfying the NR Schroedinger equation can be 
found in Ref . . One can show that the results remain valid for Dirac fermions and spinless 
particles if \(j)e^^{x,t)\'^ is replaced by the corresponding probability or flavor charge density 



A typical experiment which tries to observe particle oscillations between two flavors, 
assumed generically to be the flavors z/g and z/^, measures the flux of particles in the 
detector localized at some distance L from the source which produces particles of flavor Ue. 
The time of the measurements is not known. Usually typical experiments last hours, days 
or even years (like the observation of solar neutrinos). So the most appropriate way to flnd 
the probability (or number of particles) to cross the (theoretically) closed surface dV of the 
detector is to integrate the probability current density over the surface and integrate the 
result once more over the duration of the measurements. To that end, we make use of the 
conservation law for the total current, for two flavors Ue and z/^. 



Notice Eq. (fT7|) is exact since je + = ji + j2) despite Eq. (1161) approximate nature. Making 
use of the Gauss theorem and Eq. (fT7|) . we get for the sum of probabilities of Eq. (fT5|) . with 
flavors z/g and z/^. 



— |</)e,^(a;,f)|2 + V-je,A.(a3,t) «0 



(16) 



[see Eq. ffM ]. 



— {\Mx,t)\' + \(j)^{x,t)\') + V-{je{x,t)+J^{x,t))=0, 



(17) 




d'x{\M^,t)\' + \<P^{x,t)\%^ - / d'x{\Mx,t)\' + \<t>^{x,t)\%^,{18) 
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where (f)e,fi{x, t) correspond to the same scalar wave function considered in Eq. ([5]). However, 
it is a fact that the above spatial integration in the volume V is bounded by the extension 
of the detector and, essentially, by the position and localization of the wave packet. If we 
consider ti ^ to be the creation time, the last integral gives zero, and the first integral 
results in a value different from zero only when t ~ ~ -^/^, where L is the source-detector 
distance and v is the average velocity. The above probability expressions could thus be 
written as 

Vu^+,^{t,<t<U)=Vu^+u,iU-^T^L/v)= [ d^x{\Mx,L/v)\' + \^^{x,L/v)\'). (19) 

Jv 

In one-dimension analysis, considering appropriate cylindrical surfaces, we have 

/ dS-i{x,t) = J{z,t), (20) 
JdA 

and the continuity equation is reduced to 

^ +00 



d 
di 



2^ 



+ {Je{z,t) + Jf,{z,t)) =0, (21) 



SO that, from approximate conservation of the flavor current (fT6l) over distances and time 
scales much smaller than the oscillation length. 



+ 00 



Vu^,,^{z = L,U--T^L/v)=Vu^,u,{t^L/v)= I dz\4>eAz,L/v;9)\'= f dt J^^^iL , t) , (22) 

where the ^-integration has been extended from —oo to +oo because we are assuming the 
detector extension D is much larger than the wave packet width a, i. e., D ^ a 144 1. 

In the literature, the change of variables J dz — > j dtv is frequently noticed. For the 
spatial integration of the above expression, it is mathematically acceptable when 2; ~ t; t in 
one-dimension analysis, i.e., 

/ + 00 ~ + oo „ 

dz\^eAz,L/v;6)\' dt\cl>,AL,t;e)\'= dtJeAL,t). (23) 

However, it is important to remark that automatic normalization is only guaranteed for 
space integration. 

To summarize this point, it seems obvious from the above fundamental quantum me- 
chanical calculations that the spatial integration is not in confront with time integration. 
In fact, the spatial integration just concerns the wave packet localization. When one makes 
some assertion about the measurements, additional integrations over the measurement time 
and energy may (or should) be required (in order to obtain time/energy averaged values for 
Pu^,u^-, and as we have observed, it can be definitely adequate to the analysis here performed. 
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Another way of reconciling temporal oscillation with spatial oscillation makes use of 
the usual textbook connection between position and time for a free particle in Quantum 
Mechanics: the Ehrenfest's theorem. A free particle can be found to be located at mean 
position {x{t)) = (v)t at time t with error given by Ax = ((X — (a;))^) . For each mass- 
eigenstate, when dispersion can be neglected, such quantities are well defined and obey 

{xi{t)) = {vi)t , Axi w Cj = constant, (24) 

for initial position {xi{0)) = 0, where (vj) is the group velocity, the expectation value of 
Vi = p/Ei{p) in momentum space. For flavor states, in the regime of total overlap, the 
mean position should be given by vt, where v is the average of (vj), with error no larger 
than max(aj). When neutrinos are detected after traveling a distance L, which is the only 
experimentally known variable, it is licit to replace t in the formulas hj L/v with error given 
by max(aj), as long as such quantity is much smaller than the detector characteristic size 
D. 

2. The analytical approach 

Now we turn back to the IWP framework in order to obtain an analytical expression for 
(f)i{z,t). To evaluate the integral in Eq. (fTTj) . we firstly rewrite the energy Ep] as 

= E,[l + a,ia, + 2v,)]^ , (25) 

where = {m?- +Po)5, v.^ = ^ and = . The use of free gaussian wave packets is 
frequently assumed in NR quantum mechanics because the calculations can be carried out 
exactly for these particular functions and, consequently, the main physical aspects can be 
easily interpreted from the final analytical expressions. The reason lies in the fact that the 
frequency components of the mass-eigenstate wave packets, Ep} = pl/2mi, modify the mo- 
mentum distribution into "generalized" gaussian functions, easily integrated by well-known 
methods. The term pi in Ep] is then responsible for the variation in time of the width of 
the mass-eigenstate wave packets, the so-called spreading phenomenon. 

In relativistic quantum mechanics, however, the frequency components of the mass- 
eigenstate wave packets, Ep] = y/jp~+~rn^, do not permit an immediate analytic integra- 
tion. This difficulty, however, may be remedied by assuming a sharply peaked momentum 
distribution, i.e., {aEi)^^ ~ o"j ^ 1. Meanwhile, the integral in Eq. (fTT]) can be analytically 
solved only if we consider terms up to order af in the series expansion. In this case, we can 



Pz 



1 + 



E^ 
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conveniently truncate the power series 



4') 

Vz 



<7. 



1 + fT,V, + ^ (1 



Tfl ■ 



(26) 



and get an analytic expression for the oscillation probability. The zeroth-order term in 
the previous expansion, E^, gives the standard plane-wave oscillation phase. The first- 
order term, PoO";, is responsible for the slippage (loss of overlap) between the mass-eigenstate 
wave packets due to their different group velocities. It represents a linear correction to the 
standard oscillation phase jsil- Finally, the second-order term, ^o"^^, which is a (quadratic) 
secondary correction, will give the well-known spreading effects in the time propagation of 
the wave packet and will be also responsible for an additional phase to be computed in the 
final calculation. 

For gaussian momentum distributions, all the terms discussed above can be analytically 
quantified. By substituting (l26l) in Eq. (fTT]) and changing the p^-integration into a a- 
integration, we obtain the explicit form of the mass-eigenstate wave packet time evolution. 



where 



and 



exp [ - i{Eit - poz + 9i{t, z))] exp 

2m?t {z-v,ty 



am 



arctan 



2mlt 



am 



1 + 



4m^ 
a^Ef 



(27) 



(28) 



(29) 



The time-dependent quantities ai{t) and 9i{t,z) contain all the physically significant in- 
formations which arise from the second-order term in the power series expansion ( l26l) . The 
spreading of the propagating wave packet can be immediately quantified by interpreting 
a,(t) as a time-dependent width, i.e., the spatial localization of the propagating particle 
is effectively given by a,(i) which increases during the time evolution. In the NR prop- 

jsQl- For times t ^ a^m, the 



agation regime, ai{t) reduces to a^^i(t) 



1 + 



rt^ 



effective wave packet width a^^m becomes much larger than the initial width a. On the 
other hand, the wave packet spreading in the UR propagation regime is approximated by 
= ^\h- ~^ ~ ci- The UR spreading is practically negligible if we consider the 



same time-scale T for both NR and UR cases, i.e., af^(T) <^ af^(r). To illustrate this 
characteristic, we reproduce from [52| the time-dependence of am i^i Fig. [1] where we have 
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assumed a particle with a definite mass value m^. By computing the squared modulus of 
the mass-eigenstate wave function, 



2(z-v,t) 



(30) 



am 

we reproduce from [52| the wave packet spreading in both NR and UR propagation regimes 
in Fig.[2] which is in correspondence with Fig.[T] It confirms that the wave packet spreading 
is irrelevant for UR particles. 

Returning to Eq. (l27j) . we could interpret another second order effect by observing the 
time-behavior of the phase 6^{t,z). By taking into account the wave packet localization, we 
assume that the amplitude of the wave function is relevant in the interval \z — Vit\ < a,{t). 
Due to the 2;- dependence, each wave packet space-point z evolves in time in a different way. 
If we observe the propagation of the space-point z = v^t, the increasing function d,{t, v,t) 
assume values limited by the interval [0, Otherwise, for any other space-point given by 
z = Vit + K ai{t),Q < \K\ < 1, the phase z) does not have a lower limit. We shall show in 
the next subsection that the presence of a time-dependent phase can modify the oscillation 
character of the flavor conversion formula. Anyway, the phase 9i{t,z) is not influent on the 
free mass-eigenstate wave packet propagation as we can see from Eq. (130|) . 



3. The oscillation probability 



By evaluating the integral (fT3|) with the approximation (126|) and performing some math- 
ematical manipulations, we obtain a factored expression for the interference term (l9l). 



(Avty 



2o2 [1 + Sp'(t)]_ 



iNT(t) = DMP(i) X Osc(t), 
which contains the damping term 

DMP(t) = [1 + Sp'(t)]"^exp 
and the oscillation term 

Osc(t) = Re{exp[-iAEt -ie{t)]} 
= cos[AEt + e{t)], 

where 



^ , , t . fm'\ Avt 



Po 



and 



f arctan [SP(t)] - 2p° ^ s 



(31) 



(32) 



(33) 
(34) 
(35) 
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with p= 1- [3+ {^Y] ^, AE = E,- E^, E = y/WK and Av = vi - V2. The time- 
dependent quantities SP{t) and Q{t) carry the second-order corrections and, consequently, 
the spreading effect to the oscillation probability formula. If AE ^ E, the parameter p 

2 

is limited by the interval [1, —2] and it assumes the zero value when ^ ^ ^. Therefore, 
by considering increasing values of p^, from NR (NR) to UR (UR) propagation regimes, 
the time derivatives of SP{t) and Q{t) have their signals inverted when ^ 



AE 



and fixing 
reaches the value i. 

The suppression of the oscillating behavior, which is primarily caused by the separation 
between the mass-eigenstate wave packets, is quantified by the damping term DMP(t). In or- 
der to compare DMP(i) with the correspondent function without the second-order corrections 
(without spreading). 



T>MVws{t) = exp 
we calculate the ratio between Eq. fl32|) and Eq. (j36]) : 

DMP(t) 



(Avt)= 



2a2 



Dmp^s(<) 



1 + n-^ (AE±) 



exp 



2a6 m 



7^ k2 , 



(36) 



(37) 



The NR limit is obtained by setting = 1 and = in Eq. ( l36i) . In the same way, the UR 
limit is obtained by setting = 4 and = E. Between these limits, the minimal deviation 
from unity of Eq. (l37l) occurs when ^ | (p ^ 0). Returning to the exponential term 
of Eq. ( 132|) . we observe that the oscillation amplitude is more relevant when Avt ^ a. It 
characterizes the regime of minimal loss of coherence where the spatial overlap between the 
mass-eigenstate wave packets is significant. In such regime, we always have J^j^^''*^ 1. 

We reproduce from 5^ the Fig.Othat shows the ratio of Eq. f l37p for different propagation 
regimes, where we have arbitrarily set = 10. For asymptotic times, the time- dependent 
term Sp(t) effectively extends the interference between the mass-eigenstate wave packets 
since 

DMP(t) t^co aE _Jpl{/^Etf] ^^g^ 



■ exp 



2a2 



DMP^ys(t) {pAEty. 
but, in this case, the oscillations are almost completely destroyed by DMP(t) (see Fig.|5]). 

The oscillating function Osc(t) of the interference term iNT(t) differs from the standard 
oscillating term, cos [AEt], by the presence of the additional phase Q{t) which is essentially 
a second-order correction. The modifications introduced by the additional phase B(t) are 
discussed in Fig.H] where we have compared the time-behavior of Osc(i) to cos[AEt] for 
different propagation regimes. The effective bound value assumed by 0(i) is determined 
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by the damping behavior of DMP(i). To illustrate this flavor oscillation behavior, we plot 
both the curves representing DMP(t) and Q{t) in Fig.O We note the phase changes slowly 
in the NR regime. The modulus of the phase \Q{t)\ rapidly reaches its upper limit when 
^ > I and, after a certain time, it continues to evolve approximately linearly in time. Es- 
sentially, the oscillation vanishes rapidly. In case of (ultra-relativistic) neutrino oscillations, 
the values considered for a in Figs. HHS] are reproduced from 52| and they can be considered 
unrealistically large. 

By superposing the effects of DMP(t) in Fig.O and the oscillating character Osc(i) ex- 
pressed in Fig.m we immediately obtain the flavor oscillation probability 



sin^ (29) 



1 - [1 + Sp'(t)]~3exp 



(Avt) 



2a2 [1 + Sp'(i)] 



cos [AEt + Q{t)] 



(39) 



which is illustrated in Fig.[6l 

Obviously, the larger is the value of aE, the smaller are the wave packet effects. If it 
was sufficiently larger to not consider the second order corrections expressed in Eq. (12^ . we 
could compute the oscillation probability with the leading corrections due to the decoherence 
effect. 



sin^ (29) 



1 — exp 



(Avt) 
2a2 



cos [AEt] 



(40) 



5l| . Under minimal decoherence condi- 



which corresponds to the same result obtained by 
tions (Avt ^ a), the above expression reproduces the standard plane wave result. 



sin^ (29) 



{1 - cos [AEt]} = sin" {29) sin' 



Am' 
AE 



L 



(41) 



since we have assumed aE ^ 1. 



B. Initial flavor violation 

The two family flavor conversion formula, obtained by simply associating scalar wave 
functions for neutrino mass-eigenstates, was given by Eq. (|8]). We can rewrite such expres- 
sion, in the three-dimensional space, as 

Viue^i^^-t) = Jd^x\Mx,t)-Mx,t)f . (42) 

It is evident from Eq. fH2|) that V{iye t) ^ if 0i(a;, t) ^ (f)2{x, t) as complex functions 

of X. In particular, for t = 0, we necessarily have V^Ve ^n'l 0) 7^ if 0i(a;, 0) 7^ 4'2i^, 0). 
Therefore, we might have an initial flavor violation, without propagation, if the initial wave 
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functions associated to the two neutrino mass-eigenstates are different. This 



problem of initial flavor definition in the formulation of neutrino oscillations 51|. 



eads to the 



Three points of view can be adopted concerning initial flavor definition: (i) the v^- 
flavor state is indeed well defined by Eq. ([5]) and we should have 0i(a3, 0) = 02(2^,0) to 
guarantee (instantaneous) initial flavor definition; (ii) the definition of neutrino flavor should 
be generalized to accommodate more general initial conditions for neutrino creation; or (iii) 
exact flavor definition can not be achieved and initial flavor violation should be regarded as a 
genuine physical effect. Within the simple approach of intermediate scalar wave packets, we 
can only adopt the positions (i) and (ii), and only estimate the consequences of point of view 
(iii). We will see in the following that for a free and first quantized Dirac fermion the strict 
adoption of (i) leads to the inevitable presence of oscillating terms with short oscillation 
length Vq^^ ~ po ~ -^Osc5 addition to usual flavor oscillation characterized by the 
usual oscillation length Lq^^ = To avoid such rapid oscillations implies (iii), i.e., an 

initial non-null presence of the wrong flavor. Indeed, in the simple second quantized theory 
of Dirac fermion with mixing shown in section lV Al the adoption of (iii) is inevitable, although 
its consequence could be very small and unobservable in practice. Indeed, it is clear from 
the discussion after Eq. (fT2|) that initial flavor violation effects implied by point of view (iii) 
should be small otherwise the amplitude of flavor oscillations would be highly suppressed. 
Moreover, large effects are excluded from neutrino conversion experiments. We can estimate 
the possible effects of flavor violation from Eq. ( !T2|) . where gaussian wave packets were used. 
The probability of initial flavor violation is given by sin^2^^(aAp/4)^, for \aAp\ ^ 1. Notice 
the condition |aAj9| ^ 1 is a particular version of the condition (j)i{x,t) ^ (f)2{x,t) for 
gaussian wave packets of equal width. 

To properly investigate (iii) and quantify its effects, it is necessary to go beyond the 
description of neutrino propagation and include the interactions responsible for neutrino 
creation (and detection). In terms of gaussian wave packets, it amounts to calculate the 
quantity \aAp\ for realistic cases. Such task will be undertaken in section lV D I for neutrinos 
created in pion decay. 



III. INCLUSION OF SPIN AND RELATIVISTIC COMPLETENESS 

This section deals with the inclusion of spin and relativistic completeness in an relativistic 
quantum-mechanical derivation of the oscillation formula and some extensions. Neutrinos 
are fermions and we know that the time-evolution of a massive spin one-half particle has to 
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be described by the Dirac equation. Firstly, we show how to include the spatial localization 
in the description of a flavor oscillating system when the dynamics of the Dirac equation 
is considered. In addition, we identify some common points between our results and some 
previous results for a quantum field approach for flavor oscillations jo ill. Finally, we extend 
the discussion about initial flavor violation introduced in the previous section to the Dirac 
equation case. In particular, at the level of a quantum mechanical approach, we identify 
some relations between initial flavor violation and some predictable rapid oscillations. 

It is well known that the Dirac equation, in its first quantized version, can give a signifi- 
cantly good description of a Dirac fermion if its inherent localization is much bigger than its 
Compton wave length; usually this is associated with weak external fields. For example, the 
spectrum for the hydrogen atom can be obtained with the relativistic corrections included 
(fine structure) jxil- One of the terms responsible for fine structure, the Darwin term, can be 
interpreted as coming from the interference between positive and negative frequency parts 
zitterbewegung) of the hydrogen eigenfunction in Dirac theory compared to the NR theory 
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73] . On the other hand a situation where the theory fails to give a satisfactory physical 



0,^ 



description is exemplified by the Klein paradox 1741 1: the transmission coefficient for a 
electron moving towards a step barrier becomes negative for certain barrier heights, exactly 
when the localization of the electron wave function inside the barrier is comparable with its 
Compton wave length. 

It is also well known that relativistic covariance requires the use of both positive and 
negative energy solutions for free relativistic particles. Otherwise, as we will see, functional 
completeness is lost, i. e., we can not expand a general solution of a relativistic wave equation 
in terms of only positive or negative eigenf unctions. Both types of solutions are also necessary 
to ensure relativistic causality, i.e., the time evolution of the solutions can not extend in 
space-time beyond the future lightcone of the initial spatial distribution. 

Given the desirable properties of functional completeness and causality, which will be 
called relativistic completeness for short, we will treat in this section the flavor oscillation 
problem using the free Dirac theory in presence of two families mixing. We will see that 
the use of the Dirac equation will automatically encompass the properties of relativistic 
completeness. 

Additionally, obtaining exact solutions of a generic class of Dirac wave equations 75|-|79|, 
specially with the presence of general interaction terms, is important because in some cases 
the conceptual understanding of the underlying physics can only be brought about by such 
solutions. These solutions also correspond to valuable means for checking and improving 



23 



models and numerical methods for solving complicated physical problems. In the context in 
which we intend to explore the Dirac formalism, we can report about the Dirac wave packet 
treatment which can be useful in keeping clear many of the conceptual aspects of quantum 
oscillation phenomena that naturally arise in a relativistic spin one-half particle theory. 
In parallel, much progress on the theoretical front of the quantum mechanics of neutrino 
oscillations has been achieved, inipelled by a phenomenological pursuit of a more refined 



flavor conversion formula 
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within the quantum field formalism j9Nll|. 



19| and by efforts to give the theory a formal structure 



b_introduce the fermionic character in the study of quantum oscillation phenomena jiol . 
Isol . we shall introduce the Dirac equation as the evolution equation for the mass- 



M 

eigenstates. Hence, Eq. (j5]) becomes 

^(x) = ipi{x) cos 6 t^i + tp2{x) sin 9 

= [ipiix) cos^ 9 + iIj2{x) sin^ 9] Uf, + [ipi(x) — iIj2{x)]cos9 sin9 

= ^y^{x-e)v^ + ^ij^^{x-e)u^ (43) 

where x = {t,x) and ipiix), z = 1,2, is a spinorial wave function, normalized to unity, that 
satisfies the Dirac equation for a mass m^. The natural extension of Eq. (fT^ reads 



ipi^^ix, t = Q-e) = ^uAx) = Tpi{x, t = Q) = Tp2ix, t = 0), (44) 

which guarantees an initial pure flavor u^. Therefore, the probability that an initial pure z/g 
flavor be detected as z/^ flavor ([8]) is 

Vii^e^i^^lt) = j d\^l^{x,t)^^^{x,t) = sin2 20^[l - DFO(t)] , (45) 

where the interference term generalizing Eq. ^ for Dirac wave packets is 

DFO(t) = d\^\{x, t)ip2{x, t) . (46) 

More discussion on initial flavor definition within first quantized Dirac theory can be found in 
section llll E[ 

A. Dirac wave packets and the oscillation formula 

Let us firstly consider the time evolution of a spin one-half free particle given by the Dirac 
equation 

- m) iPix) = 0, (47) 
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with the plane wave solutions given by iIj{x) = ip_^,(x) + ip_{x) where 

i^+i^) = exp [— ipx] u{p) for positive frequencies and 

ip_{x) = exp v(ji) for negative frequencies; (48) 

p is the relativistic quadrimomentum, p = {E,p), and the relativistic energy is represented 
hy E = ^/mP~+~^. The free propagating mass-eigenstate spinors are written as 81 1 



UsKp) = — 



[2E{m + E)] 



' E+m\ 2 



V [2E{E+m)]i^" 



Vsip) = r Vn 



[2E{E+m)]^ 



(49) 




where r/i^a = ^ / ' \ y ' ^ u^{m,0) Vq = ^^(m, 0) are the spinors for zero momentum 

(they do not depend on the mass) and the relations ul{p)Ur{p) = vl{p)Vr{p) = 5rs are satisfied. 
The right-hand sides of the second equalities of Eq. (149 p assume the Dirac representation 
for the gamma matrices (see Eq. ( I125p ). When more than one mass-eigenstate is in use, the 
positive and negative frequency eigenspinors will be respectively denoted by u\p,mi) and 
Vip.mi). 

To describe the time evolution of mass-eigenstate Dirac wave packets, we could be inclined 
to superpose only positive frequency solutions of the Dirac equation. It seems, at first glance, 
a reasonable choice. However, for generic initial states it is usually necessary to superpose 
both positive and negative frequency solutions of the Dirac equation. We can show that 
both types of solution are necessary if we require the instantaneous creation of a pure Vf, 
flavor, i.e., the mass-eigenstate wave functions satisfy %l)^{z,Q) = ip2{z,0) in Eq. (I44p . (See 
sections III Bl and IIII El ) More particularly, we can adopt the initial condition 

il;^^{x,t^O;d) = (l)^^{x)w, (50) 

obeying Eq. (jB]), where is a scalar wave function and w isa constant spinor which satisfies 
the normalization condition w'^w = 1. That is not the most general condition because w 
could depend on spatial coordinates, but it simplifies the preliminary calculations and allows 
a direct comparison to the scalar treatment of section |Tll 

To continue the analysis, we turn back to the one-dimensional space {x — z) and express 
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the flavor wave function ipuei^^ 9) in terms of 



•^-oo ^ s=l,2 

+ d''*(-p^,TO,)'t;*(-p^,mi)exp[+i£'(p^,TOi)t]}. (51) 
Therefore, since ^/'^(zjO) = tpi,^{z,0,9) for t = 0, the Fourier transform of Eq. (1511) should be 

[b\p^,m,)u''ip^,mi) + d''*{-p^,mi)v''{-p^,mi)] = ip{p^ - po)w , i = l,2, (52) 

3 = 1,2 

where we used Eq. (I50p and denoted if{p^ - Po) as the Fourier transform of (puei^)- (I^i fact, the 
normahzed Fourier transform of 4>v^{z) is ip{p^ - po)/\/2tt.) Using the orthogonahty properties 
of Dirac spinors, we flnd 7l| 

b''{p,,mi) = ip{p^ - po)u'''<{p,,mi)w, 
(f*{-p^,mi) = Lp{p^ -po)v''\-p^,mi)w. (53) 

These coefficients carry an important physical information. For any initial state which has 
the form given in Eq. (l50l) . the negative frequency solution coefficient d^*{-p^,mi) necessarily 
provides a non-null contribution to the time evolving wave packet. This obliges us to take 
the complete set of Dirac equation solutions to construct the wave packet. In the general 
case, one can succeed to choose ^^^{z,0,9) in such a way that ■?/'i(z,0) = i/jt,^{z,0,6) has 
vanishing coefficients d'^*{-p^,m^) = but then ip2{z,0) = ip^^{z,0,9) will necessarily contain 
non- vanishing coefficients d^*{-p^,m2) as it is shown in section UlI E[ 

Having introduced the Dirac wave packet prescription, we are now in a position to calcu- 
late the flavor conversion formula. The following calculations do not depend on the gamma 
matrix representation. By substituting the coefficients given by Eq. (l53l) into Eq. (15T1) and 
using the well-known spinor properties 7l| . 

y v'{-p„m,)v'{-p,,m,) = — ^ , (54) 

i=l^2 2E{p,,m,) 

we obtain 

ip^{z,t) = [ ^ ipip,~po) exp [ip,z] (cos [E{p„m,)t] - + "t- J ^-^ [E{p^,m,)t]\ w . (55) 

27r [ E{p,,m,) J 



By simple mathematical manipulations, the new interference oscillating term in Eq. (jl6 
will be written as 

-l^f'iP- - Po) {[1 - /(Pz,™i,2)] cos[e_(p„mi,2)t] 

■oo 
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+ f{Pz,nii^2) cos[e+{p^,mi^2) t]} (56) 

where 

r, , , E{p„m,)Eip„m2)-pl-m^m2 

/(Pz,TOi,2) = J{Pz,mi,m2) = — 



2 E{pz,mj_)E{pz,m2) 

and 

e±{pz,mi^2) = e±{pz, nil, 1712) = E{pz,mi) ± E{pz,m2) . 

The time-independent term /(p^,mi,m2) deserves some comments. It has a minimum at 
Pz = and two maxima at pz = ±y/rrijn^. We can readily observe in Fig.[8]that it goes 
rapidly to zero when p^ ^ mi^2 (UR limit) as well as when po <^ mi^2 (NR limit). It 
means that when we consider a momentum distribution sharply peaked around po ^ mi^2 
or Po ^ mi 2 the corrections introduced by /(p^,mi,m2) are negligible. The maximum value 
of /(p^,mi,m2) is 

_ 1 _ y/mim2 _ [y/rnl - 

J max — „ I ~ I N ' V^'j 

I nil + 1712 2[mi + 771-2) 
which vanishes in the limit m-i = mj. 

The effects introduced by f {p^, mi, 1112) are relevant only when Am mi 3> ma. Mean- 
while, what is interesting about the result in Eq. fl56p is that it was obtained without any 
assumption on the initial spinor w or Fourier transform ^^(p^-po)- Otherwise, the initial 
spinor carries some fundamental physical information about the created state. And this 



could be relevant in the study of chiral oscillations [82, |83| where the initial state plays a 
fundamental role. In comparison with the standard treatment of neutrino oscillations done 
by using scalar wave packets, where the interference term iNT(i) is given by Eq. (fT3l) with 
AE{p^) = e_(p2, mi, ma), we notice in DFO(i) two additional terms. In the first one, the 
stonc/arc? oscillating term cos [e_(p2,rni,2) t], which arises from the interference between mass- 
eigenstate components of equal sign frequencies, is multiplied by a new factor obtained by 
the products u\pz,mi)u{pz,m2), mj f (-p^, mj) and h.c. The second one is a new os- 

cillating term, cos [e+{pz,mi,2) t], which comes from the interference between mass-eigenstate 
components of positive and negative frequencies. The factor multiplying such an additional 
oscillating term is obtained by the products u\pz,mi)v{~Pz,m2), v^-Pz,mi)u{pz,m2) and h.c. 

The new oscillations have very high angular frequencies e+. Indeed, we can find the 
relation e+e_ = Am^, where e_ is the usual (double of) flavor oscillation angular frequency. 
Such a peculiar oscillating behavior is similar to the phenomenon referred to as ZBW. In 
atomic physics, the electron exhibits this violent quantum fluctuation in the position and 
becomes sensitive to an effective potential which explains the Darwin term in the hydrogen 
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atom 72j. We shall see later that, at the instant of creation, such rapid oscillations introduce 
a small modification in the oscillation formula. 



B. The oscillation formula with first order corrections 

A more satisfactory interpretation of the modifications introduced by the Dirac formalism 
is given when we explicitly calculate DFO(t). To that end, we resort to the same gaussian 
wave packet used in section [TTl and use Eq. (HM for the scalar part (p,,^ (z) of the initial wave 
packet in Eq. (150 p : 



i/ju^iz, t :^0- 9) = (t),y^{z)w = ( exp 



exp [ipoz] w . 



(58) 



However, in contrast to section [III where we have considered the energy E{p^,mi) expansion 
up to the second order terms in Eq. (126|) . so that the spreading effects were included in 
the analysis, we focus this preliminary study only on first order corrections. Thus, we 
approximate the frequency components by 



E{p,,m,) ^ E, +v, {p, -Po) ■ 
Considering the approximation (15^ . we can write 



/(p^,mi,TO2) 



1 - V1V2 ( 1 H — ] + V1V2 

Po 



(v?+vD 1 + 



Pi 



P. -Po 
Po 



(59) 



(60) 



and 



e±(p,,mi,m2) w E^^ ± + {v-, ± V2) {p, - Po) . (61) 

For UR particles (m^ ^ po), we can also use the following expression for the central energy 
values (Ei) and the group velocities (v,) of the mass-eigenstate wave packets. 



This implies 



Po + 



Am 



m 



2po 



2po 

e+{p,,mi,m2) ^ 2po 
e_(p2,mi,m2) 



and V, ^ 1 — 



P. -Po 



Po 

^ ^ ml + ml ^ p. 



Po 



2po 



1 - 



4Po' 

P. -Po 
Po 



Po 



ml + ml 



where {Amy = {m^ — ma)^ is different from Am^ = ml — ml which appears in the standard 
oscillation phase. 
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Finally, by simple algebraic manipulations and after gaussian integrations, we find for 
Eq.®, 



DFO(i) w exp 



1 - 

+ exp 



/ Am 
V2po 

2o? 



COS 



Am^ 



2po 
2 \ 2 



+ 



Am 
2po 



Am^ 
2po 



^f+m|\ 



+ 



2pot 
(apoF 



/ Am 
V2po 

m^+mj 



COS 



sm 



mj+m| 



2pg 



}. (62) 



As we have already noticed, the oscillating functions going with the second exponential 
function in Eq. ( l62l) arise from the interference between positive and negative frequency 
solutions of the Dirac equation. It produces very high frequency oscillations which is similar 
to the quoted phenomenon of ZBW j72|. The oscillation length which characterizes the very 
high frequency oscillations is given by L^^J^ ~ Obviously, L^^f is much smaller than 
the standard oscillation length given by L^f^. = It means that the propagating particle 

exhibits a violent quantum fluctuation of its flavor quantum number around a flavor average 
value which oscillates with L^f^. Meanwhile, except at times t ~ 0, it provides a practically 
null contribution to the oscillation probability as we can observe in Fig. El 

To explain such a statement, let us suppose that an experimental measurement takes place 
after a time t ^ L for UR particles. The observability conditions impose that the propagation 
distance L must be larger than the wave packet localization a. Since the (second) exponential 
function vanishes when L ^ a, for measurable distances, the effective flavor conversion 
formula will not contain such very high frequency oscillation terms, and can be written as 



Pn 



sm' 



exp 



/ Am^L 
\2V2apl 



1 



Am 



2po 
Am\^ Am^ 



cos 



-L 



2po 



L sin 



2po 
Am'^ 



+ 



-L 



.(63) 



For distances which are restricted to the interval a <^ L <^ a , the interference term 

Am^ ' 

is preserved due to minimal slippage between the wave packets. In this case, we could 
approximate the oscillation probability to 



sin^ (29) 



1 



1 



/ Am^L 
\2V2apl 



1 



Am 
2^ 



cos 



~2^ 



L 



■ (64) 



However, we reemphasize that it is not valid for T ^ L ~ when the rapid oscillations 
are still relevant (L < a) . By comparing the result of Eq. (IMl) with the scalar oscillation 
probability of Eq. (1?!]) , we notice a deviation of the order (^^)^ that appears as an additional 
coefficient of the cosine function. It is not relevant in the UR limit as we have noticed after 
studying the function /(p^,mi,m2). 
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C. Time evolution operator 



An initial value problem, described by a differential equation, can be completely solved if 
we find a well defined kernel that governs the time evolution. Then from an initial field (wave 
function) configuration we can find the subsequent configurations at all subsequent times. 
Moreover, it is possible to analyze the properties of time evolution that are independent 
of the choice of the initial wave packet. Analyzing the properties of the kernel, it is also 
possible to check the properties of relativistic completeness for the free propagation of Dirac 
fermions with flavor mixing. Therefore, we will explicitly calculate, in this section, the time 
evolution operator (kernel) in presence of mixing in order to extract the common features 
independently of the initial wave packet. The extension to treat three families of neutrinos 
is straightforward. A matricial notation will be used throughout this subsection to express 
the mixing. Since the one-dimensional restriction does not lead to any formal simplification, 
the full three dimensional space will be considered here. 

In matricial notation the mixing relation that extends Eq. 1^, between flavor wave func- 
tions ^E'J(a;) = {iIj1^{x),iIjJ^{x)) and mass wave functions \E'J'^(a;) = {ipj (x) , 1IJ2 (x)) , is 

(cos9 s'mO \ 
^m{x) . (65) 
— sin^ cosO I 

Each mass wave function is defined as a four-component spinorial function ipn{x, t), n = 1,2 
that satisfy the free Dirac equation (H71) 

d 

i—'4}n{x,t) = H^^n{x,t), n = l,2, (66) 
where the free Hamiltonian is the usual 

= -ia-V + , n = 1, 2 . (67) 

The normalization condition is 

J d'x'^l{x)^m{x) = J cfx^\{x)^f{x) = 1, (68) 

which implies that the wave functions ipiix) and ^/'2(a^), or, for the same reason, tpi,^{x) and 
ipiy^lx), can not be simultaneously normalized to unity. A different notation and normaliza- 
tion convention was being used through this review and it deserves a further clarification. 
As it can be seen in Eqs. and , the wave functions for the mass-eigenstates were nor- 
malized to unity. To recover Eq. (H3|) . the wave functions ipi and ip2 used in this subsection, 
as in Eq. (165|) . should be replaced by 

V'l — )■ tpi/cosO, V'2 — ^ ijj2/ s'mO . (69) 
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This correspondence should be clear from the first line of Eq. ([5]). Equation (168|) is a 
restatement of the fact that the probability to find a neutrino over all space irrespective of 
mass or irrespective of flavor is unity. We also see that condition (l68l) is time independent. 

We will work in the flavor diagonal basis. This choice deflnes the flavor basis vectors 
simply as 

I/e ^ Z>e = (1, 0)^ , Uf,^V^ = (0, 1)"^ , (70) 

while the flavor projectors are obviously 

^u^^y^i^i. (71) 

Actually, as an abuse of notation, the equivalence U ~ U ® Id is implicit, as well as, 
Py^ ~ Py^ 1_d; the symbol la refers to the identity matrix in spinorial space. 

The total Hamiltonian governing the dynamics of is = diag(iff, H^). From the 
considerations above, '^f{x,t) satisfy the equation 

i^^'/(a;,t) = UH^U-^^f{x,t) . (72) 

The solution to the equation above can be written in terms of a flavor evolution operator 
as 

= 0) = /A' A-,. 0), (73) 



K°{x-x-;t)= I ^ K"{pi t) . (74) 



where 

We can calculate K^{t) in any representation (momentum or position) as 

e~^^P^ + sin^O e~^^°^ — cos9 sm9{e~^^^^ — e"*^?* 
- cos^ sme{e~'^°' - e-'^°') sin^O e"*^?* + cos^^ e-'^?' 



(75) 



It is important to emphasize that Eq. (173!) leads to a causal propagation because the kernel 
in Eq. ( 174|) respects K^{x — x';t) = for space-like distances \x — x'\ > \t\. This property 
follows directly from Eq. (17S|) and the properties of causality of each kernel, {x\e~''^"^\x') = 
for |a; — a;'| > \t\, associated to the Dirac Hamiltonian with mass n = 1,2 [see Eqs. (I270p 



and flAlSj) or Ref. 



The conversion probability is then 

V{u,^i.^;t) = Jd^x^{x,0)K''\t)P,^K''{t)^f{x,0) 
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= ld'p^l{p)iK^e)^KeiP,t)i'M , (76) 

satisfying the initial condition \[^J(a;,0) = {ipj^lx , 0) , 0) . Such initial condition implies, 
in terms of mass eigenfunctions, iIji{x,0) = cos9 ip^^^x) and ip2{x.,^) = sm6ilj,y^{ 



requirement to obtain an initial wave function with definite flavor, as seen in section[TlB]|5]J. 



The function ip^^{p) denotes the inverse Fourier transform of ipy^{x) (see Eqs. (lAip and (lA2p ). 

Before obtaining the conversion probability for Dirac fermions, let us replace the spinorial 
functions ipn{x) by spinless one-component wave functions ipn{x) in the flavor wave function 
\I^J(x) — 7- {({)yXx),Lpy^{x)) and mass wave function ^1^{x) — )• {ipi{x) , {p2{x)) . We also 
replace the Dirac Hamiltonian in momentum space H^ijp) (167|) by the relativistic energy 
En{p) = a/p^ + rn^- Inserting these replacements into Eq. 0761) we can recover the usual 
oscillation probability 
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5l| 



d'p\K'(p,t)^,M\' 



|2 



d'pnP,t)\^.M\, (77) 

where ^/(a;,0)"^ = (y?^^(a;)"^, 0), K^^ip,t) = {K%i = - sin ^ cos ^(e'^^i^P)* - e-^^^b)*) ^nd 

^{p, t) = sm^29 sm\AE{p)t/2) (78) 

is just the standard oscillation formula ([T]) with AE{p) = Ei{p) — E2{p). The conversion 
probability ( 177|) in this case is then the standard oscillation probability smeared out by 
the initial momentum distribution. If the substitution \(py^{p)\'^ — )■ 5^{p — po) is made the 
standard oscillation formula is recovered: it corresponds to the plane-wave limit. 

After we have checked the standard oscillation formula can be recovered for spinless 
particles restricted to positive energies, we can return to the case of Dirac fermions. We can 
obtain explicitly the terms of the mixed evolution kernel fl75p by using the property of the 
Dirac Hamiltonian in momentum space {H^Y = ^nip) ^d, which leads to 

+ sm^29f{p)sm\Eip)t)lD , (79) 

where 

/(P) = ill - . (80) 

E{p) = |[Ei(p) + Ei(p)] (851 and ^{p,t) is the standard conversion probability function (1780 . 



The expression (IHOj) was already found in Eq. ( !56|) in a slightly different form. A unique 
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implication of Eq. (179|) . which is proportional to the identity matrix in spinorial space, is 
that the conversion probability (1761) does not depend on the spinorial structure of the initial 
flavor wave function but only on its momentum density as 

V{u,^v^;t) = Jd'p{:^{p,t)[l-f{p)] 

+ sm'29f{p) sm\Et)} i^l{p)i^.Ap) ■ (81) 

(The tilde will denote the inverse Fourier transformed function throughout this subsection.) 
Furthermore, the modifications in Eq. (18T|) compared to the scalar conversion probability 
(177|) are exactly the same modifications found in Eq. ( 156 p if we replace ip'^{p^ ~ po)/2t{ by 
'4'tSp)4'ue{p) cLiid use three-dimensional integration. 
The conservation of total probability 

V{u^^u^;t)+V{u^^u^;t) = l , (82) 

is automatic in virtue of 

K^Jip, t)Kg{p, t) + K^Kp, t)K^,{P, t) = lD, (83) 

and initial normalization in Eq. (168|) . The survival and conversion probability for an initial 
muon neutrino are identical to the probabilities for an initial electron neutrino because of 
the relations 

t) = Kg\p, t)KZ{p. t) , (84) 
K^}{p, t)K^,{p, t) = Kg\p, t)K^^{p, t) . (85) 



D. Spinless neutrinos 

The derivation of the usual conversion probability (177|) takes into account only the positive 
frequency contributions. The mass wave function used to obtain Eq. (1771) corresponds to 
the solutions of the wave equation 

i—ip{x,t) = V-V^ + m'^Lp{x,t) , (86) 



which is equivalent to the Dirac equation in the Foldy-Wouthuysen representation j86| . 
restricted to positive energies. The evolution kernel for this equation is not satisfactory 
from the point of view of causality js^, i.e, the kernel is not null for spacelike intervals. 
Moreover, the eigenfunctions restricted to one sign of energy do not form a complete set 



73| 
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To recover a causal propagation in the spin case, the Klein-Gordon wave equation must 
be considered. In the first quantized version, the spectrum of the solutions have positive 
and negative energy as in the Dirac case. However, to take advantage of the Hamiltonian 
formalism used in section lTlI C\ it is more convenient do work in the Sakata-Taketani (ST) 
Hamiltonian formalism [73] where each mass wave function is formed by two components 

'^n{x,t)= f '^"^'^'^^ I , n=l,2. (87) 

\Xn(x,t) J 

The components ip and x are combinations of the usual scalar Klein-Gordon wave function 
4>{x) and its time derivative do4>{x). This is necessary since the Klein-Gordon equation 
is a second order differential equation in time and the knowledge of the function and its 
time derivative is necessary to completely define the time evolution. An analysis of the 
flavor oscillation of scalar wave functions obeying the Klein-Gordon equation can be found 
in Ref. 

The time evolution in this formalism is governed by the Hamiltonian 73 1 



X/2 

= -(t3 + ^^2)^ + ml , (88) 

which satisfies the condition {H^^Y = (— + m^)!^^, like the Dirac Hamiltonian (1671) . 
The Tk represents the usual Pauli matrices and 1st is the identity matrix. 
A charge density 



88] can be defined as 

= ^ivi^n = Wn? - \Xn? , (89) 

which is equivalent to the one found in Klein-Gordon notation do(j). Needless to say, this density 
(j89p is only non-null for complex (charged) wave functions. The charge density is the equivalent 
of fermion probability density i/j^ip in the Dirac case, although the former is not positive definite 
as the latter. The adjoint $ = ^^t^ were defined to make explicit the (non positive definite) norm 
structure of the conserved charge 

J A t) = = time i„depe„de„. . (90) 

Consequently, the adjoint of any operator Q can be defined as = t^Q^^ts, satisfying 
((]$,$) = ($,(]$). Within this notation, the Hamiltonians of Eq. (IHHI) is self-adjoint, 
H^^ = H^"^, and the time invariance of Eq. (pOl) is assured. 

We can assemble, as in the previous section, the mass wave functions into \E'^ = ($]") ^2) 
and the flavor wave functions into \I'J = ($31) "^^v); satisfying the mixing relation ^E'/ = U'^m- 
Then, the time evolution of can be given through a time evolution operator K'^'^ acting 
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in the same form as in Eq. (173|) . In complete analogy to the calculations from Eq. (17^ to 
Eq. fl76p . we can define the conversion probability as 

= Jd'p'¥Jp)K^K'^^ip,t)^,ip) , (91) 

where '^f{x,0)'^ = ($e(a^)"'', 0). The adjoint operation were also extended to ^/ = \E'j(l6) ® 
Ts), where Ig is the identity in mixing space. 

The information of time evolution, hence oscillation, is all encoded in 



+ sm^2ef{fip)sm^{Et)lsT , (92) 



where the function f{p) were already defined in Eq. (IHOl) and 



/I mi 1712 ,oo\ 

l^ = \7;i — + — )• (93) 

V 2 1712 TTT'l 

The factor > 1 determines the difference with the Dirac case in Eq. (179|) . The equality 
/i = 1 holds when mi = m2, i.e., when there is no oscillation. 

Therefore, rapid oscillations are also present in flavor oscillations of charged (in the sense 
of Eq. ( l89l) ) spin particles, with contributions slightly different from the Dirac case. Such 
result reinforces the fact that rapid oscillations are direct consequences of flavor mixing, 
(a) relativistic nature of the wave equations governing time evolution and (b) initial flavor 
deflnition. The presence of a spinorial degree of freedom is not a requirement to rapid 
oscillations. 



E. Initial flavor violation and rapid oscillations 

We saw in Eqs. (156|) and ( 18T]) that the flavor conversion probability for Dirac fermions 
presents rapid oscillation terms with frequency 2E = Ei{p) + E2{p) , in addition to the usual 

2 2 

oscillation frequency AE = Ei{p) — E2{p) = Ei^)+E2{p) ' ^^""^ frequencies arise naturally 
from the propagation of free particles with mixing in a relativistic classical fleld theory that 
contains states with positive and negative frequencies ±Ei{p), e.g., spin 1/2 fermions and 
spin bosons. It was assumed, however, that flavor can be well defined and a pure z/e flavor 
is created at t = 0. We will elucidate here that such rapid oscillations can be avoided but 
only at the expense of not having an initial fiavor exactly defined. More specifically, we will 
see there is a clash between two choices: 
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• cond. A: exact initial flavor definition (pure flavor creation) and 

• cond. B: standard flavor oscillation, without rapid oscillations. 

Firstly, following the formalism of subsection IIII C\ we calculate the probability that a 
generic superposition of neutrinos Ui and i>2, denoted simply by z/, described by a generic 
be detected as a neutrino z/^: 



where the operator inside, in the mass-basis, can be calculated explicitly 

\- sine cos Oe'^^^e-^^?^ cos^ 9 

Notice in this section we are working in the mass-basis, keeping the same notation as in 
section lTlI CI where {Pu^,)ij = ^fi^^ji = U*iU^j, since (^'q,)^ = U*^ and (z>j) = 6ij in this basis. 



Only the off-diagonal terms depend on time t, and energies Ei and E2, which means the 
functional dependence on time and energies, in momentum space, should be of the form 
g{AEt, Et), where AE = Ei{p) - E2{p) and 2E = Ei{p) + E2{p). Such functional form is 
in accordance to Eq. ( ISTl) . for example. To disentangle the dependencies on AEt and Et, it 
is instructive to rewrite the free Dirac time evolution operator, in momentum space, in the 
form 

^-iHSt ^ e-ii^nt^D^ + e^^"*A^_ , (96) 

where 

1 

A^± = -(lz.±-^), (97) 

are the projector operators to positive (+) or negative (-) energy eigenstates of . By using 
the decomposition above ( p6|) . we can separate the different contributions in the off-diagonal 
terms of Eq. (jHS]) by rewriting 

+ e*2^*Af+A^_ + e-*25*Af_Af+ . (98) 

The same procedure can be applied to its hermitian conjugate. The time evolution kernel 
of Eq. (I75|) . in the flavor basis, can be analyzed in the same fashion. Since A{^A^ 7^ 0, 
it can be seen that the rapid oscillating terms come from the interference between, e. g., 
the positive frequencies of the Hamiltonian and negative frequencies of the Hamiltonian 
iff. 
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We know that both positive and negative energy states should propagate properly to 
preserve causality. With a kernel containing only positive energy contributions e~*^"*, there 
is no causal propagation js^. Therefore, to maintain causality, we can not eliminate, for 
example, the negative energy contribution in the kernel of Eq. (p6|) . One may, of course, 
restrict the initial wave functions to contain only positive energy states. Such restriction 
eliminates the rapid oscillatory terms, also known as zitterbewegung, for a one-particle free 
Dirac theory j71|. For two masses, however, the positive energy eigenf unctions with respect 
to a basis characterized by a mass mi necessarily have non-null components of negative 
energy with respect to another basis characterized by 1712, as we will see. 

Let us examine exact initial flavor definition (cond. A). Firstly, we introduce wave func- 
tions ipi and iIj2, normalized to unity, by using 

(cos a -01 \ vr 
, < a < - ; (99) 
sin a V'2 / 

where additional phases can be easily incorporated in ipi,'ip2- Then Eq. yields 

V{u^u^; t) = sin^ie + a) sin^ds*) + sin'^{e - a) cos^dSt) , (100) 
where we used the following parametrization for Eq. (H6l) . 

cosEt = Re jd^x'ipl{x,t)^lj2ix,t) , < < vr , (101) 

which is consistent with |cosSf| < 1 for normalized ipi and ■?/'2- Obviously ipi{x,t) = 
e~^^i'*ipi[x). Since the expression in Eq. (llOOp is a sum of non- negative terms, to have 
V{v^v^; 0) = at t = 0, we should have 

6* - a = , and Hq = . (102) 

The other possibilities, 6 + a = 01 d + a = T[, are outside the parameter range of a since 
< < 7r/2 for non-trivial mixing. Moreover, a = 9 corresponds to the only local and 
global minimum of the second term of Eq. fllOOp as a function of a. The function sin^(a; + 9) 
in the first term of Eq. ( llOOl) is always non-null and has the global minimum at the borders 
a = 0ora = 7r/2. The second equality in Eq. (I102p is equivalent to 

i^i{x,{))=i)2{x,^)=^{x). (103) 

Thus the necessary and sufficient initial conditions for initial pure z/g flavor or, equivalently 
(for two flavors), no z/^ flavor, restricts Eq. (p9|) to 

'u:,ihi{x) 

^>m{x) =1 ' ' I = U:Mx)h = ^[x)v, , (104) 



.f^e*2V'2(^), 
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where in the last equahty, Eq. fllOSp was exphcitly assumed. Although condition fllOSp can 
be also inferred from Eq. (175]) . restricted to Eq. (1681) . we showed an alternative derivation 
without assuming a fixed proportion between mass-eigenstates, which is a more general 
initial condition than the one used in section llTBl 

Let us now try to eliminate rapid oscillations (cond. B) by adjusting the initial conditions. 
This was not attempted so far. We begin by considering only positive energy states for 

^l^iix) = 4'^\x) , i = l,2, (105) 

in Eq. (1^1?]) . In momentum space, Eq. (11051) can be written 

4'^\p) = Y,u%p,mi)gi{p,s) , (106) 

s 

where the spinors were defined in Eq. (1491) . By construction, Kij^il)\'^\x) = iIj^~^\x) and 
Ai-iljl'^\x) = 0, which means ilji{x,t) = e~^^^*ilJi{x) only contains the positive energy term 
^-iEi(p)t pQ^j-igj. transform. In this case, the time dependent function in the z/^ 

detection probability (llOOp is given by 

cos St = d'^pi>l{p)tP2{p)e'^^^''^' (107) 

= Re d'pglip, s)g,ip, s)^^^^^^'[l - 2h{p)] , (108) 



where 



/2(p) = \- ^^/TT^^/TT^ - ^ Vr^\/r^ , (109) 



defined with the shorthand Xi = ttt-t- We used the calculation 

u''\p,m,)u'{p,m,) = [l-2f2{p)]6rs. (110) 

We achieved our goal: there is no rapid oscillations. The probability V{v ^iy^]t) is only a 
function of AE{p)t smeared out in momentum. 

Let us now show that restriction (I102p for cond. A can not be simultaneously applied 
with restriction (llOSp for cond. B. 

Let us suppose cond. A [Eq. (I103p ] is valid, i.e., ipi{x) = ip2{x) = ip{x). We can 
decompose the spinorial function ip{x) in terms of bases depending on different masses mi 
and 1712- Equating 

^Pix) = j^^[ut{x;p)g^\p,s)+vt{x-p)gl\p,s)], z = l,2, (111) 
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where ^^^=^ = u'{p,m,) e'P''^ /{Iny^^ and ^%=^ = v'{p,mi) e-^P'^ / {2-Kf'\ the expansion 



coefficients can be obtained 



g'^\p,s) = I d'x^^^i^ix) , (112) 



gl\p,s) = jd^x^-^^{x) . (113) 

From Eq. (I113p . and ip{p) being the inverse Fourier transform of ip{x), we see that imposing 
the conditions 

g'f\p,s) = v'\p,m,)ij{-p) = , (114) 

g'^~\p,s) = v'\p,m,)i,{-p) = Q , (115) 

for all p, leads to the conditions 

[(mi + E2) - (m2 + E2)\ vl^{-p) = 0,3 = 1,2, (116) 



1 1 



vl^l.p^{-p) = 0, s = l,2, (117) 



where the property of Eq. (lAlOp and 70^0 = —Vq were used. In case nii 7^ m2, we can use 
the decomposition ip{p) = i'+ip) +4'-{p), where ip±{p) = (1 ±7o)V^(p)/2, and obtain from 
Eqs. flll6p and f lll7p the conditions 

vo^ip-i-p) = 0, s = l,2, (118) 
M^tS-pV^+(-p) = 0, s = l,2, (119) 

where the properties 7 = 7075^] and Uq = 'J5V^ were used in Eq. flll9p . Since S-p has 
only non-null eigenvalues and it commutes with 70, the equations (IllSp and f lll9p are only 
satisfied if ip+^p) = ip-ip) = 0, i.e., ip{x) = 0. It is easier to reach this conclusion in 
the helicity basis {mo"^', fQ^^'} characterized by "S-pul^^^ = ±\p\u'^\ but the result is basis 
independent. 

We can then conclude that rapid oscillations are an inevitable consequence of a relativistic 
Dirac fermion description of flavor oscillations if initial flavor definition is exactly required. 
Quantitatively, however, it was seen in Eqs. (l56l) and (1791) that the contribution of rapid 
oscillations to the probability V{ve^t^ii', t) is negligible for momentum distributions around 
UR values because rapid oscillation is quantified by the function / in Eq. (IHOl) . which 
behaves as / ~ (Am)^/(4i?^) for UR momenta. On the other hand, it was already shown 
that avoiding rapid oscillations is possible to the detriment of having initial flavor violation. 
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Let us quantify the amount of initial flavor violation for neutrinos states of the form fllOSp 
satisfying cond. B. Let us relax the conditions of Eq. fll02p . and consider only the first one, 
a = 9, which is a minimum in a. Thus the initial probability to detect the wrong flavor 
can be rewritten by using Eq. (11081) as 

^^'^^^^g'"^ =\Y. fd'p{\9i{p,s)-g2{p,s)\' + 4Mp)Re[gl{p,s)g,{p,s)]], (120) 

s 

where we used the expansion fll06p . 

We immediately see that taking equal momentum distributions for the two mass- 
eigenstates (for two spin states), i.e., 

giiP,s) = g2iP,s) = gip,s) , (121) 
makes the first term of Eq. (11201) vanish. Initial flavor violation is then proportional to 

'^^'^'-Xf^ = ld'p\g{p)?Up). (122) 

sm 26 J 

where 1(7(^)1^ = Xls spin-independent momentum distribution. Notice the 

normalization 

jd'p\g{p)\^ = 1. (123) 

As it should be, /2(p) = if uii = 1712 and it behaves as /2 ~ (Am)^/ (4^)^ for UR momenta. 
Moreover, f2{p) has the unique maximum value at |p| = y^m^m^, the same point as f{p)- 
In fact, the function /2 and the function /(p) in Eq. ( IHOl) are closely related by the relation 

[l-2/2(p)]2 = l-/(p). (124) 

Therefore initial flavor violation in Eq. (I122p is also negligible for momentum distributions 
around UR momenta, provided that Eq. fll2ip is valid. 

We can conclude that the description of flavor oscillations for a free Dirac fermion sub- 
jected to mixing, within first quantization, is satisfactory from the quantitative point of 
view but in the strict sense one should consider either one of the following phenomena as 
inevitably present: (i) initial flavor violation or (ii) rapid oscillations. Both effects might be 
as negligible as Eqs. (I122p or ( !79|) (or Eq. ( 156|) ). quantified by the functions f2{p) or f{p), 
respectively. We will see (i) is supported by second quantization. In that case, it is possible 
that such effect could be still small but larger than Eq. (11220 . All depends on how different 
are ipi{x) and iIj2{x), associated to the two mass-eigenstates, when a neutrino is created as 
a superposition of mass-eigenstates. Furthermore, the presence of (i) suggests that neutrino 
flavor is not an exact concept but only an approximately well defined one. 



40 



IV. CONSEQUENCES OF SPIN STRUCTURE AND RELATIVISTIC COM- 
PLETENESS 



This section deals with some intrinsic consequences of spin structure and relativistic com- 
pleteness, namely chiral oscillations and the consequences for a non-minimal coupling with 
external magnetic fields. Chiral oscillations naturally enter the discussion of neutrino oscilla- 
tions since neutrinos are produced and detected through weak interactions that are chiral in 
nature. We succeed in showing that the inclusion of chiral oscillation effects, together with 
the time-evolution of spinorial wave packets for the mass-eigenstates, can modify the flavor 
conversion probability formula. In particular, it provides an interpretation of chiral oscilla- 
tions as very rapid oscillations in position along the direction of motion, i.e., longitudinal 
to the momentum of the particle. All the ingredients for decoupling chiral oscillations from 
the spin-flipping in the presence of an external magnetic fleld are discriminated. It allows 
for depicting all the relevant effects due to the inclusion of spin structure and relativistic 
completeness that we have introduced. 

According to the SM at tree level, neutrinos interact with itself and other leptons only 
through weak charged current and neutral current interactions which are chiral in nature and 
only contain the contribution of neutrinos with negative chirality Consequently, positive 
chirality neutrinos become sterile with respect to V-A weak charged currents, independently 
of any external electromagnetic fleld. Despite the experimental circumstances not being 
favorable to such an interpretation, the quantum transitions that produces a flnal flavor 
state corresponding to an active-sterile quantum mixing is perfectly acceptable from the 
theoretical point of view. Even focusing speciflcally on the chiral oscillation mechanism, we 
have also observed a generic interest in identifying the physical meaning of variables which 
coexist with the interference between positive and negative frequency solutions of Dirac 



equation 
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3H92] in a Dirac wave packet treatment. If we follow an almost identical line of 
reasoning as that applied for rapid oscillations of the position, it is also possible to verify 
that the average value of the Dirac chiral operator also presents an oscillating behavior. 

We intend to clarify some conceptual relations between helicity and chirality by construct- 
ing the dynamics of the processes of chiral oscillation and spin-flipping. The differences 
between the dynamics of chirality and the helicity for a neutrino non-minimally coupled 
to an external (electro) magnetic fleld are expressed in terms of the equation of the motion 
of the correspondent operators 7^ and /i, respectively. In order to determine the relevant 
physical observables and confront the dynamics of chiral oscillations with the dynamics of 
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spin-flipping in the presence of an external magnetic field, we report about a further class of 
static characteristics of neutrinos, namely, the (electro) magnetic moment which appears in 
a Lagrangian with non-minimal coupling. We shall demonstrate that the oscillating effects 
can be explained as an implication of the ZBW phenomenon that emerges when the Dirac 



equation solution is used for describing the time evolution of a wave packet j52|. Due to 
this tenuous relation between ZBW and chiral oscillations, one question to be answered in 
this manuscript concerns with the immediate description of chiral oscillations in terms of 
the ZBW motion, which shows that, in fact, chiral oscillations are coupled with the ZBW 
motion so that they cannot exist independently of each other. One shall conclude that chiral 
oscillations can be interpreted as very rapid oscillations in position along the direction of 
motion. 

In order to make clear the confront of ideas involving chirality and helicity, we first review 
such concepts and their common points. At the same time we point out their differences by 
assuming a Lagrangian with minimal or non-minimal coupling. The latter interaction, more 
precisely, the non-minimal coupling of a magnetic moment with an external electromagnetic 



field |40'] , which induces neutrino spin-fiipping (39| , was formerly supposed to be relevant in 



the context of the solar-neutrino puzzle. As a consequence of such interaction, left-handed 



neutrinos could change their helicity (to right-handed) 41|. The effects on flavor oscillations 
due to external magnetic interactions in a kind of chirality-preserving phenomenon were also 
studied [42] but they lacked a full detailed theoretical analysis. We observe, however, that 
only for UR particles such as neutrinos, changing helicity approximately means changing 
chirality. In the context of oscillation phenomena and in the framework of a flrst quantized 
theory, the small differences between the concepts of chirality and helicity, which had been 
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interpreted as representing the same physical quantities for massless particles |39| 
can be quantifled for massive particles. 



A. Chirality and helicity 

By considering the gamma matrices (in the Dirac representation), 

r = I ^ I , y = I ° / I and 7^ = I - - I , (125) 





0^ 




f (tM 


: 


















where a* are the Pauli matrices, we can observe the following 7^ matrix properties, 

(y)^=y, (y)' = l and W%r} = 0, (126) 
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which lead to the commuting relation [7^^, S^'^] = 0, where S^'^ are the generators of the 
continuous Lorentz transformations. Among the Dirac field bilinears constructed with the 
gamma matrices, it is convenient to define the vector and the axial vector (or pseudo-vector) 
currents respectively by 

f{x) = ii{x)-i^'4)(x) and f'^{x)=ii{x)-1^'-f'•'^l){x). (127) 

By assuming that ip{x) satisfies the free Dirac equation, we verify that the vector current 
density j^{x) is conserved, i.e., 

d^f{x) = {imlp{x))il) + i^{-imip{x)) = 0. (128) 

When we couple the Dirac field to the electromagnetic field, j'^{x) will become the electric 
current density. In a similar way, we can compute 

df.fix) = 2im^(a;)7Xa;). (129) 

If m = 0, the axial vector current j^^{x) is also conserved. 

The two currents j'^ix) and j^^{x) are the Noether currents corresponding respectively to 
the two gauge transformations 

il^ix) — >• exp [i a\'il){x) (130) 

and 

'0(2:) — )• exp [z a 7'']V'(a;). (131) 

The former is a symmetry of the Dirac Lagrangian. The latter, called chiral transformation, 
is a symmetry of the kinetic part of the Dirac Lagrangian jsil. The Noether theorem confirms 
that an axial vector current related to a chiral transformation is conserved, at the classical 
level, only if m = 0. 

In general, it is useful to define the chiral left-L (negative chiral) and right-R (positive 
chiral) currents through the combinations of the operator 7^^ by 

i^(x) = V^(x)7^^^V'(x) and jAj(^)=v;(x)7'^^^V(x). (132) 

Just when m = 0, these are respectively the current densities of left and right-handedpaiticles 
and they are separately conserved. 

In parallel, we can define the operator h = |5] • p, so that, by conveniently choosing 
the two component spinors rj^, the spinors u^{p) and Vg{p) become eigenstates of h with 
eigenvalues /i = +i for s = 1 and h = for s = 2 {h = c = 1). Although the free-particle 
plane wave solutions can always be chosen to be eigenfunctions of h, it is not possible to 
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find solutions which are eigenf unctions of S • n with an arbitrary unitary vector n. It 
occurs because the operator S • n does not commute with the free particle Hamiltonian 
unless n = ±p or p = 0. The operator h that can be diagonalized simultaneously with the 
free particle Hamiltonian is called the helicity operator 72|. It is exactly the particle spin 
projection onto the momentum direction. The particle with /i = +| is called right-handed 
and that one with h = is called left-handed. The helicity of a massive particle depends on 
the reference frame, since one can always boost to a frame in which its momentum is in the 
opposite direction (but its spin is unchanged). For a massless particle, which travels at speed 
of light, it is not possible to perform such a boost. In a certain way, helicity and chirality 
quantum numbers carry a kind of complementary information. In vacuum, the chirality is 
invariant under a continuous Lorentz transformation but it is not constant in time if the 
particle is massive. The helicity is constant in time but it is not Lorentz invariant. A chiral 
eigenstate can always be written as a linear combination of two helicity eigenstates. 

To focus our attention on the central point of this section, where we intend to compare 
the dynamic processes of chiral oscillation and spin-flipping, let us analyze the possibili- 
ties of observing them via interaction terms in the Lagrangian language. Initially, we can 
construct a gauge invariant Lagrangian in a theory with initially massless fermions. We 
can couple a Dirac fermion to a gauge field by assigning the chiral-L fields ipi^{x) to one 
nontrivial representation of a gauge group G and assigning the chiral- i? fields ipini^) ^ 
singlet representation. For such a model, we can write 



r 



V'(x), 



(133) 



where it is evident that only negative chiral fields couple with gauge boson fields v4^(a;). It is 
straightforward to verify that the Lagrangian (I133p is invariant by the following infinitesimal 
local gauge transformation jsi l 



l + ial{x)T- 



1 — 7" 



+ g-'d^a^{x) + r'M^(x) a'(x), 
where the non-abelian character of G is summarized by the commutation relations 



rpa rj-ib 



j^j^abcrpc 



(134) 



(135) 



Since the chiral- fields are free fields, we can even eliminate these fields and write a gauge 
invariant Lagrangian for purely chiral-L fermions. The idea of gauge fields that couple only 
to chiral-L fermions plays a central role in the construction of the theory of weak interactions. 
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To work out the general properties of chirally coupled fermions, it is useful to analyze 
one further transformation. Let us take the charge conjugate operator C = i7^7'' in order 
to write a charge conjugate field ip'^ as 

so that the conjugate chiral-i? component of a field ip, which would be given by ipl, trans- 
forms as a chiral-L quantity under a continuous Lorentz transformation, i.e., 

= mi = i (^^) 7>* = ^l' (^^) r = il'^R , (137) 
where we have omitted the x dependence for simplicity. In the same way, the conjugate 
chiral-L component of a field ■?/', which would be given by ^/l)^, transforms as a chiral-i? 
quantity, i. e., 

= il'^L ■ (138) 
This property is important for observing what occurs in a gauge model with L-R symmetry 
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, |94|. In this class of models, chiral-i? neutrinos interact with additional vector gauge 
bosons. 

As a general feature, due to the vectorial character of the interactions which are mediated 
by vector gauge bosons, the possibility of a left to right (or vice-versa) chiral conversion via 
gauge coupling is impossible. As we can observe in the following, the probability amplitude 
for a vector coupling between fields with opposite chiral quantum numbers is null, 

Aii^LA^r^R)^ = ^J^A^l''^L = ((1 + 7')^)N"^/.7''(1-7')V' 

= 7" ^1^7^^ (1 + 7^)(1 - 7^)^^ = 0, (139) 

which forbids the possibility of chiral changes (oscillations) in this way. Alternatively, we 
cannot abandon the idea of spin-inversion which does not correspond to chiral conversion but 
is a very well established mechanism which emerges via an external vector field interaction. 

At last, we can analyze what happens to the vector currents in the presence of flavor 
mixing. Let us define the flavor currents 

f.^ = ^..ri^.. , (140) 

where the mixed fields are 

i^u^=Uaiipi, a = e,iJ,. (141) 
Then, the flavor currents are almost conserved because 

= ' E U:,U^.{m, - m,)^^^, . (142) 
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For two flavors, we have 

d^ji^^^,^ = ±(m2 - mi) sin2^Im(^i7/>2) , (143) 

where U is given by Eq. ([6]). Therefore, approximate flavor conservation is vahd locally. The 
non-conservation is proportional to the mass difference and to the amount of mixing. 



B. Chiral oscillations coupled with the zitterbewegung motion 

If a massless particle represented by a Dirac spinor is produced in a chiral (7^) eigenstate, 
its chirality is a constant of the motion in addition to being Lorentz invariant. On the other 
hand, independently of the existence of external interacting flelds, these properties no longer 
hold for Dirac particles with mass. In vacuum, chiral oscillations can be formally related to 
the zitterbewegung (ZBW) motion. 

We begin by calculating the average value of the spatial part of the vector current in 



Eq. fll27p . Such average value appears quite reasonable 72| when we calculate it using the 
positive frequency solutions {ip^ix)), 

{cx{t))^ = jd^i^\{x)cxi^^{x) = j^^^Yl\^M\\ (144) 

which represents exactly the average value of the relativistic velocity ^. However, for gen- 
eral wave functions %Ij{x) containing both positive and negative frequency components, an 
extremely rapid oscillation term appears in Eq. fll44p . To understand such behavior, note 
that the velocity oper ator Q:(t), in the Heisenberg picture, does not represent a constant of 



the motion since 
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j^cx{t) = i [H,, ait)] =2i{p- a{t)H,) . (145) 



Keeping in mind that p and Hq are constants o 



differential equation for (x{t) and obtain |72| 



the motion, we can easily solve Eq. (I145p 



cx{t) = pH-^ + (q;(0) - pH-^) exp [-2 iHot]. (146) 

The time dependence of a{t) in Eq. ( I146p translates into the trembling motion (ZBW motion) 
in the position operator x{t), since ^^^^ = cy.{t). Such motion is characterized by high 
frequencies > 2m, u being the eigenvalues of 2Ho. 

We have already seen in the previous section that the particle helicity /i = ■ p is a 
constant of the motion while the chiral operator 7^ is not. At time t = 0, the coefficients b^ip) 
and dlip) used in the construction of the Dirac wave packet ipix) can be chosen to provide 
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a negative (positive) chirality eigenstate 19|, |80| or, in the same way, to provide a helicity 



eigenstate [when hu^{p){v^{p)) = ±^u^{p){v^{p))]. We can consider the initial chiral eigen- 
state ip{0,x) is not only a superposition of momentum eigenstates weighted by a momentum 
distribution centered around (like in Eq. (155|) ). but also a helicity eigenstate obtained 
through the production of a spin-polarized particle. Such condition formally occurs when 
we assume that the constant spinor w in the wave packet expression ( l55l) is a simultaneous 



eigenspinor of 7'^ and h 146| . We then obtain 



^(75(0) = ^(a(t)-p(S-p)) 

j^{a{t) . p)) (S • p) + {ait) . p) (^(E • p)) . (147) 

By observing that the helicity h and the momentum p are constants of the motion, we can 
establish a subtle relation between the chirality operator 75 (t) and the velocity operator a.{t) 
expressed in terms of 

^(75(0) = m (^(a(i)-p)) =(2/i)(^a(t).p>. (148) 

The time evolution of 75 (t) presents an oscillating character which can be interpreted as a 
direct and natural consequence of the oscillating character of ct{t) in Eq. fll46p . Considering 
Eqs. (1145P and (1148^ leads to the following explicit dependence 

j^ilrM = {4thip-cxit)H,)-p). (149) 

Since h, p and Hq are constants of the motion, we conclude that there will not be chiral 
oscillations without the "quivering motion" of the position. 

The relation between the chirality operator 75 (i) and the velocity operator a.{t) becomes 
more interesting when we take into account the complete expression for the vector current 
density 'ip{x)'~^^'il){x) and apply the Gordon decomposition jsil, 

^{x)^^^{x) = [(9^?(x)) ^[x) - ^[x) {d^^Pix))] + (V^(x)a^Xa:)) , (150) 

where a^^ = |[7^,7^]. The spatial integration of the vector components of Eq. f llSUp gives 
us 

j d^^^a^j = -^jd"^ {-i [?(V^) - (VV^) V] + [Vx (V^S^) - idt {i^aij)] } , (151) 
where we have suppressed the x dependence. We can write the decomposed components of 
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(a) as 

jd^xVx{^-Stp) = 0, (152) 

+ E (5 " 3 <^sP[Kip)d:ip) exp[+2iEt] 

-dMUP) exp[-2iEt]] 1, (153) 



where p = {E,—p) and we have assumed that = 7]\,Jl-pr]^ = (— l)''"''^^^^/. Finally, the 
last term is written as 

y"dfe5t(V^«V) = /(^jE (^)a.PK(p)<(p)exp[+2iSt] - dMbsiP) exp[-2iEt]] 

+ ^n,K(p)d:,(p) exp[+2iEt] - d^ip)bAp) exp[-2iEt]] I , (154) 

where ni(2) = 1 ± i2 when p = 3, with 1, 2 and 3 corresponding to arbitrary mutually 
orthonormal vectors. 

Eq. ( I152p allows us to reach the conclusion that the ZBW does not get a contribution 
from the intrinsic spin dependent (S) magnetic moment component which couples to the 
external magnetic field B{x). In fact, the ZBW originates from the current strictly related to 
the internal electric moment. By taking into account the (momentum direction components 
of the) Eqs. (1153111541) . we can turn back to Eq. (I149p . which carries the main idea of this 
subsection, and observe that chiral oscillations can be essentially constructed in terms of the 
longitudinal components of (ex.). By calculating the mean value of {a. ■ pHo), we obtain 

{cx{t)-pH,) = l^i^^.pJ^[{E)\bM\' + {E)\ds{pm 

+ Yl '^<^s[{E)b:[p)d*M exv[+2iEt\ - {E)dAp)bM exp[-2iEt]] I 

s=l,2 ) 

= (IpD- l^rn^aAb:{p)d:{p) exp[+2iEt] - dMbsiP) exp[-2i E t]](,155) 
which can be substituted in (I149p in order to give 

^(75(i)) = f-^'^Yi2ha,)i2iE)[b:ip)d:ip) exp[+2iEt] - d,{p)bM exp [-2 i E t]] ,{156) 
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so that the time evolution of the chirahty operator can be written as 

(75W> = (75(0)) + /(^f i'^MbsiP) {e^p[2iEt]-l) + h.c.] 

= (75(0))+ /t^^ ^[<i.(p)6.(p)(exp[2ii5;t]-l) + M (157) 

since 2ha^ is equal to one for well defined spin up or spin down (helicity) states. 

More subtle physical interpretations of the above results are considered in 46|, |82|, 183] . 
For a gaussian initial wave function with null average chirality, the oscillations of the left- 
right (L-R) chiralities cancel and there is again no overall oscillation. It could be the origin 
of an apparent paradox. We observe that for any mass-eigenstate represented by plane wave 
solutions of the massive Dirac equation, the rest-frame wave function is always an equal 
mixture of both chiralities. This is easily verified when we write 

^ = + = ^L + ^R, (158) 

where '4^l,r correspond respectively to chirality quantum numbers =1=1. It could be shown 
that, in the rest frame of a particle, we have 

iV'Lr = Hr?- (159) 

Note that this result is not Lorentz invariant since a Lorentz boost is not a unitary trans- 
formation. Thus, while the cross section is Lorentz invariant the chiral probabilities are 
not. This seems to suggest that probability measurements are chiral independent. We seem 
to have an argument against the physical significance of chiral oscillations. The reply to 
this objection based upon the Lorentz invariance is simply that in any given Lorentz frame 
chiral oscillations are manifestly important because of the chiral projection form (V-A) of 
the charged weak currents. The chiral probability variations produced by Lorentz transfor- 
mations (even if 7^ commutes with the Lorent z gen erators) are automatically compensated 
by the wave function normalization conditions 147l | and the Lorentz transformations of the 
intermediate vector bosons and other participating particles. 

In the context of electroweak interactions, the determination of the explicit form of the 
helicity although possible, is not as relevant as the determination of 75, since, as already 
pointed out, the helicity is not a Lorentz invariant quantity. On the other hand, electroweak 
interactions single out the negative chirality state as the interacting state. The same is not 
true for the helicity quantum number (spin projection quantum number). From the phe- 
nomenological point of view, chiral oscillations as well as spin-flipping can independently take 
place in combination with a third physical process, for instance, the flavor conversion of a 
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propagating neutrino. In fact, chiral oscillations coupled with flavor oscillations can induce 
some small modifications to the standard flavor conversion formula Q|. More specifically, 
the spinorial form and the interference between positive and negative frequency components 
of the mass-eigenstate wave packets can lead to the coupling of chiral oscillations and flavor 
oscillations which effectively introduce some small modifications to the standard flavor con- 



version formula |52| when applied to the study of NR neutrinos. Therefore, in some cases, 
the proper use of the Dirac equation as evolution equation for the mass-eigenstates might 
be necessary for a satisfactory description of fermionic (spin one-half) particles, even if one 
discusses Weyl or Major ana spinors 95 1. 



C. Fermionic particles non-minimally coupled to an external magnetic field 



The spin-flipping attributed to some dynamic external interacting processes j39| was for- 
merly supposed to be a relevant effect in the phenomenological context of neutrino physics, 
where the external interaction had its origin in the non-minimal coupling of a magnetic 
moment to an external electromagnetic field {4^. On the other hand, independently of any 
external electromagnetic field, since neutrinos are detected essentially via V-A charged weak 
currents, the chiral oscillation mechanism by itself may even suggest an explanation for the 



LSND anomaly j24j . l28l |29| . Such an erroneous correspondence between chiral oscillation 
and spin-flipping causes confusion in the literature where the concept of helicity has been 
currently used in the place of chirality. Consequently, neutrinos with positive chirality are 
sterile with respect to weak interactions. We are now interested in pointing out the mod- 
ifications of the chirality dynamics which are observed when the neutrino electrodynamics 
is accurately discussed. In particular, we are interested in the (electro)magnetic moment 
that appears in a Lagrangian with non-minimal coupling. We determine the equation of the 
motion for the relevant physical observables obtained from such a Lagrangian so that we 
can confront the dynamics of chiral oscillations with the dynamics of spin-flipping. 

Despite their electric charge neutrality, neutrinos can interact with a photon through 



loop (radiative) diagrams 96|, |97|. The Lagrangian for the interaction between a fermionic 

field tp{x) and an electromagnetic field written in terms of the field- strength tensor F'^'^{x) = 

d^'A^ix) - ^''A^'{x) is given by 

1 -, 
2 



£ = - V;(x) a^^ [fi F^^(x) - dTf^'ix)] ^{x) + h.c. (160) 



where x = {t,x ), a^j^t, = |[7/x5 7v]; the dual field-strength tensor is given by J^^^^{x) 



^e'^"^^ Fxs{x) |l48 | and the coefficients n and d represent, respectively, the magnetic and 
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the electric dipole moment which estabhsh the neutrino electromagnetic coupling. One can 
notice that we have not distinguished the flavor/mass mixing elements in the above inter- 
acting Lagrangian because we are indeed interested in the physically observable dynamics 
governed by the Hamiltonian 

H = a-p + l3m-/3 {^iF^"" {x) - dF^"" {x)) + h.c. 

= OL-P + p[m- ^iJ:- B{x) + dJ:- E{x)], (161) 

where, in terms of the Dirac matrices, ot = Yl^^i^k^ ~ Yll=ilolk^^ 1^ = loi and B{x) and 
E{x) are respectively the magnetic and electric fields. In fact, Eq. (116ip could be extended 
to an equivalent matrix representation with fiavor and mass mixing elements where the 
diagonal (off-diagonal) elements described by yU, ,,(771^ ,,) and rf^ ,,(m,^j), where i, j are mass 
indices, would be called diagonal (transition) moments. In this context, for both Dirac 



and Majorana neutrinos, we cou 
and electric dipole moments 



d ha ve transition amplitudes with non-vanishing magnetic 



99| - |l0l| . Otherwise, CP invariance ensures vanishing diagonal 



electric dipole moments [47|. Specifically for Majorana neutrinos, it can be demonstrated 
;hat the diagonal magnetic and electric dipole moments vanish if CPT invariance is assumed 



99|. 



Turning back to the simplifying hypothesis of diagonal moments and assuming CP and 

CPT invariance, we can restrict our analysis to the coupling with only an external magnetic 

field B{x) by setting = 0. In a minimal extension of the standard model in which neutrinos 

are described as massive Dirac particles and without new interactions accessible below the 

weak scale, the neutrino magnetic moment arises at the one-loop level, as does the weak 

contribution to the anomalous magnetic moment of a charged lepton [7l|. The value of n 

can be read off from general formulas for the electromagnetic vertex, to one-loop order, of 

an arbitrary fermion {€}. To leading order in m\/m\Yi the results are independent of irti and 

of the mixing matrix so that it turns out to be proportional to the neutrino mass (matrix), 

3eG 3meG „ ,„ ruu 

LL = — ^m = — usm^ = 2.7 X lO^'Vs — (162) 
8^/27r2 4V27r2 „ mjv 



149|. In particular, for 



where G is the Fermi constant and m-Ar is the nucleon mass 
1 eV, the magnetic moment introduced by the above formula is exceedingly small to be 
detected or to affect physical or astrophysical processes. 

Since we are interested in constructing the dynamics governed by the Hamiltonian of 
Eq. ( 116ip . we firstly observe that the free propagating momentum is not a conserved quantity, 

j^ip) = i{[H,p]) = Re(A.)(/3V(5].S(x))). (163) 
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In the same way, the particle velocity given by 

j^{x) = i{[H,x]) = (a) (164) 

comes out as a non-null value. After solving the differential equation for x{t) (or for oc{t)), it 
is possible to observe that, in addition to a uniform motion, the fermionic particle executes 
very rapid oscillations known as zitterbewegung 72 1. 

By following an analogous procedure for the Dirac chiral operator 75, recurring again to 
the equation of the motion, it is possible to have the chirality and the helicity dynamics 
respectively given by 

^(7^) = 2i(7o75["^-/^S.S(x)]) (165) 

and 

^(^> = ^ (70 [(S • V)(S • Bix)) + 2(S X Bix)) ■ p]) (166) 

where we have alternatively defined the particle helicity as the projection of the spin angular 
momentum onto the vector momentum, h = ^T: ■ p (with p in place of p). From Eqs. (11651 
I166P we can state that if a neutrino has an intrinsic magnetic moment and passes through 
a region filled by an magnetic field, the neutrino helicity can flip in a completely different 
way from how chiral oscillations evolve in time. 

In the non-interacting case, it is possible to verify that the time- dependent expectation 
value of the Dirac chiral operator 75 has an oscillating behavior j46| very similar to the rapid 
oscillations of the position. The Eqs. (116511166]) can be reduced to the non-interacting case 
by setting B{x) = so that 

j^{h) = i{[H,h]) = -((axp).p) = (167) 



and 



= i{[H,j']) = 2im(7V>. (168) 



from which we confirm that the chiral operator 7^ is not a constant of the motion j46|. The 
effective value of Eq. (I168P appears only when both positive and negative frequencies are 
taken into account to compose a Dirac wave packet, i.e., the non-null expectation value of 
(ToTs) is revealed by the interference between Dirac equation solutions with opposite sign 
frequencies. The effective contribution due to this interference effect lead us to resort to the 
Dirac wave packet formalism in order to quantify neutrino chiral oscillations in the presence 
of an external magnetic field. 
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Effectively, the above discussion represents the first step for accurately deriving the ex- 
pression for the neutrino spin-flipping in magnetic fields which can be related to chiral 
oscillations in the limit of a massless particle (UR limit). By correctly distinguishing the 
concepts of helicity and chirality, we can determine the origin and the influence of chiral 
oscillations and spin-flipping in the complete flavor conversion formula. 

In the previous section IIIIl we have confirmed that the fermionic character of the parti- 
cles, more specifically, the presence of positive and negative frequency components and the 
requirement of pure flavor creation, modify the standard oscillation probability obtained by 
assuming a scalar nature of the mass-eigenstates. The new oscillation probability contained 
an additional oscillating term with very high frequency, proportional to the sum of energies 
of the mass-eigenstates. Such modifications, with neutrinos being described by free Dirac 
wave packets, introduced correction factors that, under the current experimental point of 
view, were not effective in the UR limit of neutrino propagation in vacuum. Nevertheless, 
the case of n eutr i nos moving in a background matter, where spin/chiral effects become more 



relevant 



102 



103] . deserves, at least, a careful investigation. The physical consequences 



in environments such as supernova can be theoretically studied 104| . For instance, it was 
observed that neutrinos propagating in matter achieve an effective electromagnetic vertex 
which affects the flavor conversion process in a situation where chirality can be preserved 42l |. 



1. Chiral oscillations in the presence of a magnetic field 



Assuming the simplifying hypothesis of a uniform magnetic field B, the physical conse- 
quences of the non-minimal coupling to an external magnetic field can be studied by means 
of the eigenvalue problem expressed by the Hamiltonian equation 



H{p) v3„ = {a • p + /? [m - ^ S • B]} ip„ = E„{p) ip„, 
for which the explicit 4x4 matrix representation is given by 



m - fiB, -fJ,{B^ - iBy) p^ - ipy 

-fj.{B^ + iBy) m + fj,B^ p^ + ipy -p, 

Pz Px-iPy -{m-iiB,) n{B^-iBy) 

p^ + ipy -p, fi{B^+iBy) -{m + p.B,) 

The most general eigenvalue {E„{p)) solution of the above problem is given by 



(169) 



(170) 



K(p) = ±^s(p) = ±Jm^ + + {-iy2y^^^^^^iFTJp^^, s = 1,2, 



(171) 
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where we have denoted i?„=i,2,3,4 = ±-£-3=1,2 and we have set a = jiB. The complete set 
of orthonormal eigenstates Lp„ thus can be written in terms of the eigenfunctions Uijps) 
with positive energy eigenvalues {+EXp)) and the eigenfunctions V{ps) with negative energy 
eigenvalues {—E,{p))^ 



U{p.) = - 
V(p.) = - 




(172) 



We denoted as the relativistic quadrimomentum, = {Eg{p),p), N{p^) as the normaliza- 
tion constant and 

Af = A^ip) ± 2m\a\ - a\ a± = 2^,(p)|a| ± (A^(p) + a^) 

with 

/^l{p) = El{p)-{m'+p^) + a\ (173) 

We can observe that the above spinorial solutions are free of any additional constraint, 
namely, at a given time t, they are independent functions of p and they do not represent 
chirality/helicity eigenstates. 

In order to describe the above Hamiltonian dynamics for a generic observable 0{t) we 
can firstly seek a generic plane-wave decomposition as 

exp [—i{Eg{p) t — p ■ x)]U{ps), for positive frequencies and 

exp[ i{E^{p) t — p • x)]V{ps), for negative frequencies, (174) 

so that the time-evolution of a plane-wave packet ipit^x) can be written as 

i^{t,x) = J^Y.{bip^)Uip^)exp[-iEMt] 

+ d*{ps)V{ps) exp[+iE,{p)t]}exp[ip-x], (175) 

with p^ = {E^, — p). We see Eq. fllTSp requires some extensive mathematical manipulations 
to explicitly construct the dynamics of an operator 0{t) in the form 

O,0 = /.".,..(,x)O«<.«). (176) 

However, if the quoted observables like the chirality 7^, the helicity h or even the spin 
projection onto B commuted with the Hamiltonian H, we could reconfigure the above 
solutions to simpler ones. To illustrate this point, let us limit our analysis to very restrictive 
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spatial configurations of B so tliat, as a first attempt, we can calculate the observable 
expectation values which appear in Eq. (I165p . Let us then assume that the magnetic field 
B is either orthogonal or parallel to the momentum p. For both of these cases the spinor 
eigenstates can be decomposed into orthonormal bi-spinors as 



(177) 



and 



V{p.) = N-{p^) 



(178) 



(Ps) 
X~{Ps) 

Eventually, in order to simplify some subsequent calculations involving chiral oscillations, 
we could set Lp^^ and X12 eigenstates of the spin projection operator S ■ B, i.e., besides 
being energy eigenstates, the general solutions hl{ps) and V{ps) would become eigenstates of 
the operator S ■ B and, equivalently, of S ■ a. 

Now Eq. fll69p can be decomposed into the pair of coupled equations 



{±E, - m + E • a) 
{±E^ + m - S • a) 



(179) 



where we have suppressed the dependence. By introducing the commuting relation 
[5]-p, Yi-B] = which is valid when p x B = 0, the eigenspinor representation can be 
reduced to 



and 



with m. 



K{Ps) 



V{Ps 



2E, 



Es +m. 



(180) 



(181) 



m 



-l)''|a| and the energy eigenvalues 



(182) 



In this case, the closure relations can be constructed in terms of 



Y,U{ps)®uHpsH. = E{^^ 

3=1,2 s=l,2 ^ 



1 - (-l)'5]-d 



1 - {-lyU-a 



(183) 
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Analogously, by introducing the anti-commutation relation {S-p, H-B} when p - B = 0, 
the eigenspinor representation can be reduced to 



and 



Eo + m 
2eo 



So + m 



2eo 



Sp ^ 
— r— V 

Sp 



£0+m 

x: 



with ^0 = VpM-^tj^ and the energy eigenvalues 

±E^ = ±[eo + i-ir\a\]. 
In this case, the closure relations can be constructed in terms of 



(184) 



(185) 



(186) 



s=l,2 



m 



2eo 



1 - i^iyjo^-a 



-7mPo + m 
2e, 



1 - (-l)^7o5]-d 



(187) 



where = {eo,p)- 

Since we can set (p^^ = Xi2^^ the components of an orthonormal basis, we can immedi- 
ately deduce the orthogonality relations 

uHp^)U{p^) = vHp.)V{p,.) = 6^^, 

uHps)i,V{p.) = VHPs)l,U{Pr) = (188) 

which are valid for both of the above cases. 

Finally, the calculation of the expectation value of 75 (i) is substantially simplified when 
we substitute the above closure relations into the wave packet expression of Eq. fll75p . To 
clarify this point, we suppose the initial condition over 4j{t,x) can be set in terms of the 
Fourier transform of the weight function 

yp{p~p,)w= Y,{biPs)l^iPs)+d*{ps)V{p.)} (189) 

s=l,2 

SO that 

f 

t/jiO^x) = J j^-p^ip{p~^p,)ex];)[ip-x]w (190) 

where w is some fixed normalized spinor. By using the orthonormality properties of 
Eq. ([HH]), we find [Til 

biPs) = ^p{p-p^)U^p,)w, 
d*{p.) = ^{p-p,)VHps)w. (191) 
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For any initial state '4^{Q,x) given by Eq. fll90p . the negative frequency solution coefficient 
d*(ps) necessarily provides a non-null contribution to the time-evolving wave packet. It 
obliges us to take the complete set of Dirac equation solutions to construct a complete 
and correct wave packet solution. Only if we consider the initial spinor w being a positive 
energy {E^{p)) and momentum (p) eigenstate, the contribution due to the negative frequency 
solutions d*{p^) will become null and we will have a simple expression for the time-evolution 
of any physical observable. By substituting the closure relations of Eqs. f ll83p and ( I187p into 
the time-evolution expression for the above wave packet, Eq. (I175p can thus be rewritten as 



j (^'^^^ ~ ^""P [ip-x\Y,{ [cos [E, t] - if sin [E, t]\ (i^^^f^) }w (192) 



for the first case where is given by Eq. ( 1182p and H, = a ■ p + '^oirt^ , or as 



J^ipip^P.)e^p[ip.x]J^{[cos[E^t]-i^sm[E^t]\ (M^l^^)}^; (193) 



for the second case where Sq is given by Eq. fll86p and Hq = ex. ■ p -\- 7om. 

Once we have assumed the neutrino electroweak interactions at the source and detector 
are (left) chiral (^ip'j^(l — 'j^)^Wfj^, only the component with negative chirality contributes 



to the propagation. It was already demonstrated that, in vacuum, chiral 
introduce very small modifications to the neutrino conversion formula 19 



oscillations can 



46|. The prob- 



ability of a neutrino produced as a negative chiral eigenstate to be detected after a time t 
can be summarized by the expression 

1 



d^xilj^'it^x)- —ip{t,x) 



(l-(75)(t)). 



(194) 



2 ' ' 2 

From this integral, it is readily seen that an initial —1 chiral mass-eigenstate will evolve in 
time changing its chirality. 

By assuming the fermionic particle is created at time t = as a — 1 chiral eigenstate 
(75W = -w), for the case where [S-p, "E-B] = {B parallel to p), we can write 



w 



(27r 

t 



{ [75 cos^ [E, t] + sin [2 E, t] + i^^^ sin^ [E, t] 



Hs 75 H, 



(-1) 
(-1) 



w 



(2vr) 



i?? + # cos 



(195) 
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where we have used the wave packet expression of Eq. fll92p . In the second equahty, we 
made use of 

w^sW = -1, w^H,, j^]w = and H, 75 H, = - m]. (196) 

The expression fllQSp can be reduced to a simpler one in the non-interacting case, Eq. fll57p . 

Due to a residual interaction with the external magnetic field B we could also observe 
chiral oscillations in the UR limit. However, from the phenomenological point of view, the 

2 

coefficient of the oscillating term goes with ^ which makes chiral oscillations negligible for 



UR neutrinos 



19 



80|. As a toy moG?e/ illustration, by assuming a highly peaked momentum 



distribution centered around a NR momentum <^ m^, where the wave packet effects are 
practically ignored, the chiral conversion formula can be written as 

V{ul ^ fn;t) ~ ^ (^l-cos[2mt]cos[2|a|t] - sm[2mt]sin[2\a\t\w^T.-av?j (197) 

where all the oscillating terms come from the interference between positive and negative 
frequency solutions which compose the wave packets. 

Turning back to the case where {S-p, S--B} = {B orthogonal to p), we could have a 
phenomenologically more interesting result. By following a similar procedure and mathe- 
matical manipulations, we could write 



2 



[[Ho, 75] sin [E, + E,] - {Ho, -f,}^o^ ■ asm [E, - E,]]} 



w 



(-1) j -^^"i P - pO{cos [E, t] cos [E, t] + sin [E, t] sin [E, t]] 



(-1) 



' l-0yy{p-P.){^cos[2[a[t] + ^cos[2eot]], (198) 

where we have used the wave packet expression of Eq. fll93p and, in addition to ^"''75^ = — 1, 
we have also observed that {Hq, 75} = 275S ■ p and, subsequently, w'^'S ■ p7o S ■ aw = 0. 

2 

Now, in addition to the non-interacting oscillating term ^cos[2eot], which comes from 



.2 

the interference between positive and negative frequency solutions of the Dirac equation. 



we have an extra term which comes from the interference between equal sign frequencies. 
For very large time scales, such term can substantially change the oscillating results. In 
this case, it is interesting to observe that the UR limit of Eq. fll98p leads to the following 
expressions for the chiral conversion formulas, 

V{ul ^ UR;t) ^^{1 - cos [2 [a[t]) and ^(1/^, ^ i/^,; t) ^ (1 + cos [2 |a| t]) (199) 
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which, differently from chiral oscillations in vacuum, can be phenomenologically relevant. 
Obviously, we are reproducing the consolidated results already attributed to neutrino spin- 
flipping [47] where, by taking the UR limit, the chirality quantum number can be approxi- 
mated by the helicity quantum number, but now it was accurately derived from the complete 
formalism with Dirac spinors. 



D. Restrictive possibility of analysis 



Assuming that neutrinos are produced with completely random spin orientation before 
interacting with the magnetic fields, we can take the neutrino initial state as a collection 
of particles where 50% are characterized by spin-up states and the remaining 50% by spin- 
down states. We can understand the result provided by Eq. ( 1195^ as a much more interesting 
phenomenological instrument 16| than the result expressed by Eq. (I198p . By averaging both 
the above results for (75) on time, we can easily observe that the corresponding average value 
for the case where p is perpendicular to B [Eq. fll98p ] is zero, while the same quantity for 
the case where p is parallel to B [Eq. fll95p ] depends on the polarization of the initial state 
as follows: 



By supposing the feasibi 
neutrino polarization state 



ity o f the construction of a detection apparatus which identify the 



150] . the last result becomes phenomenologically relevant since it 



distinguishes the chirality conversion rates for spin-up (s = 1) and spin-down (s = 2) states. 
Once we have assumed the neutrino electroweak interactions at the source and detector 
are left (or negative) chiral, only the component with negative chirality contributes to the 
effective detection result. The residual contribution due to the chiral oscillation mechanism 
appears in the final time averaged result where it vanishes. 

Let us compare the initial polarization distribution with the distribution after traversing 
the magnetic field region with p || B. Initially, i.e., before interacting with the magnetic 
fields, the spin-up and spin-down states would be equally probable for any detection process 
and the ratio between the probabilities to detect the different spin eigenstates would be unity. 
After traversing the region of parallel magnetic fields, the chiral oscillation mechanism can 
modify the expectation values for the opposite polarization states in different ways, which 
reflects in a ratio between the detection probabilities of spin-up and spin-down states different 
from unity. These ratios are illustrated in Figs. [9] and [TOl 
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As one can notice, since the coefficient of the oscillating term goes with the chiral 



oscillation effects become negligible for UR neutrinos |19|,|80|. 

To have a significant impact on the quantum process, there is also a natural scale that 
is required for the magnetic field strength, i.e., a critical value of 4.41 x 10^^ G. There 
are reasons to expect that magnetic fields of such or even larger magnitudes can arise in 
cataclysmic astrophysical events such as supernova explosions or coalescing neutron stars, 
i.e., situations where an enormous ne utrin o outflow should be expected. Two classes of st ars, 
the soft gamma-ray repeaters (SGR) flOsI and the anomalous x-ra y pu lsars (AXP) jl06| are 
supposed to be remnants of such events which form magnetars 107| . neutron stars with 
magnetic fields of 10" — 10^^ G. From the theoretical point of view, the possibility of larger 
magnetic fields lO^** — 10^^ G have not been discarded yet. The early universe between the 
QCD phase transition (~ 10^^ s) and the nucleosynthesis epoch (10~^ — 10^ s) is believed to 
be yet another natural environment where strong magnetic fields and large neutrino densities 
could coexist. 

Concerning the polarization measurements, the observation of the dependence of the 
neutrino chirality conversion rate on the neutrino polarization state (measurements) can 
be converted into a clear signal of the presence of right-handed (positive chiral) neutrinos 
in the neutrino-electron scattering. In f act, in a more extended scenario, the scattering of 
neutrinos on a polarized electron target jl08 | was proposed as a test for new physi cs beyond 
the Standard Model (SM). To search for exotic right-handed weak interactions 
the strong polarized neutrino beam and the polarized neutrino target is required 
It has been shown how the presence of the right-handed neutrinos changes the spectrum 
of recoil electrons in relation to the expected standard model prediction, using the current 
limits on the non-standard couplings. In this framework, the interference terms between the 
standard and exotic couplings in the differential cross section depend on the angle between 
the transverse incoming neutrino polarization and the transverse electron polarization of the 
target, and do not vanish in the limit of massless neutrino. 

Finally, we would have been dishonest if we had ignored the complete analysis of the 
general case comprised by Eqs. ( 116911172]) where we had not yet assumed an arbitrary (sim- 
plified) spatial configuration for the magnetic field. We curiously notice the fact that those 
complete (general) expressions for propagating wave packets do not satisfy the standard dis- 
persion relations like E'^ = rn^+p^ excepting the two particular cases where E^{pY = ml+p^ 
for p X B = or el = + foT p ■ B = 0. In particular, the process of neutrino propaga- 
tion through an active medium consisting of magnetic field and plasma and the consequent 
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modifications to the neutrino dispersion relations have been studied in the hterature |110| . 
In addition, such a general case leads to the formal connection between quantum oscilla- 
tion phenomena and a very different field. In principle, it could represent an inconvenient 
obstacle forbidding the extension of these restrictive cases to a general one. However, we 
believe that it can also represent a starting point in discussing dispersion relations which 
may be incorporated into frameworks encoding the breakdown (or the validity) of Lorentz 
invariance. 

E. Absolute neutrino mass from helicity measurements 

After the confirmation that neutrinos have non-null masses and non trivial mixing among 
the various types, the knowledge of the absolute scale of neutrino masses is one of the 
most urgent questions in neutrino physics. In recent times, the greatest advances in the 
understanding of neutrino properties were boosted by neutrino oscillation experiments which 
are only capable of accessing the two mass-squared differences and, in principle, three mixing 
angles and one Dirac CP violating phase of the Maki-Nakagawa-Sakata (MNS) leptonic 



143|. 



mixing matrix 

Too large masses for the light active neutrinos may alter significantly the recent cosmo- 
logical history of the Universe. One stringent bound for the value of the sum of neutrino 
masses come from cosmology as 

^mi.<0.17eV (201) 
1 ^ 

at 95% confidence level lll| . However, it is regarded as too optimistic by most cosmologists 
because it has been derived from data sets that are inconsistent among themselves. In 
particular, the influence of Lyman-a forest data may suffer from large systematical errors. 
Nevertheless, a value of around leV would be a conservative upper bound for the sum of 
neutrino masses coming from cosmology. 

Another issue of great interest refers to the existence of heavy right-handed neutrinos 
which could explain at the same time the tiny active neutrinos masses through the see- 
saw mechanism as well as the matter-antimatter asymmetry of the Universe through the 



implementation of the mechanism of bariogenesis through leptogenesis |38l |. 

Despite of the stringent bound fl20ip coming from cosmological analysis, terrestrial direc t 
search experiments establish much more looser bounds on the effect ive neutrino masses jll2 |. 



These bounds are based on ingeniously planned experiments 113l |. but their intrinsic diffi 



culties rely on the fact that they should probe, essentially, the kinematical effects of neutrino 
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masses which are neghgible compared to other typical quantities involved in processes with 
neutrino emission. 

Although the terrestrial bounds are not as stringent as the cosmological bound in 
Eq. (1201 p . it is always desirable to have a direct measurement of neutrino masses. Two 
more reasons can be listed in favor of direct terrestrial searches: (a) cosmological bounds 
may be quite model dependent and (b) we may have access to mixing parameters through 
the effective neutrino masses. For the electron neutr ino, there are ongoing experiments 
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planning to reduce the respective bound to 0.2 eV 

The main goal of this section is to review the possibility of accessing the absolute neutrino 
mass scale through one of the most natural consequences of massive fermions, i.e., the 
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diss ociation of chirality and helicity. Such possi bility was originally suggested in Refs. 
116] and it was recently considered in Ref. 117] in order to reanalyze the possibility taking 
into account the present experimental bounds. 

Let us consider the pion decay vr^ — t- /i~z/^. Since the pion is a spin zero particle, in 
its rest frame, the decaying states should have the following form from angular momentum 
conservation, 

K)^ +6\fi:^)\i^:^) , (202) 

where the arrows represent the momentum direction (longer arrow) and the spin direction 
(shorter arrow); the combination for example, means a rightgoing fermion with a right 
helicity. On later calculations the right (left) helicity will be denoted simply as /i = +1 
{h = —1). The normalization of the state is arbitrary and the coefficient 6 is of the order of 
irLjj/Ey, which will be calculated in the following. For massless neutrinos there is no second 
term in Eq. fl202p since the antineutrino is strictly right-handed in helicity and chirality. 
Thus by measuring the wrong helicity contribution of the charged lepton it is possible to 
have access to the neutrino mass. The precision in the polarization measurement necessary 
to extract the wrong helicity will be also calculated. 

The pion decay vr^ — )■ l^Vj can be effectively described by the four-point Fermi interaction 
Lagrangian 

^cc = 2V2GFla^LUijUjJ^ + h.c, (203) 

where L = ^(1 — 75), {Uij} denotes the MNS matrix while is the hadronic current that 
in the case of pion decay reads 

Jf, = VudULlf^dL . (204) 
From Eq. fl203p it is straightforward to calculate the amplitude for n{p) Zi(q)z/j(k), 
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hj = 1,2, 

-iMin^ kvj) = IGpF^VudUu (qj^LUijV,^ (k) (205) 

by using 118 1 

{0\u-f5l,,d\Tr-) = ip^V2F^ , (206) 

where ^ 92MeV is the pion decay constant. Let us denote the spinor dependent amph- 
tude as 

Mij = uM^Lv^^{k). (207) 
J2 \M^,\' = ^p■q^){p■k,) - 2p\q,-k,) (208) 

spins 

= M^{Ml-M^) 

+ m] {Ml + 2Ml -m]). (209) 

The last expression fl209p is exact and follows when p = Qi + kj (four- vector), p"^ = M^, 
qf = Mf and k'j = The first expression (12081) does not assume energy- momentum 

conservation. 

We can calculate the squared amplitude, summed over the neutrino spin, but dependent 
on the polarization n of the charged lepton in its rest frame: 

P,,(n,) = J2 l-^^jf (210) 

Uj spin 

= Mf[qi-k, + M.ikyrii)] + 2Mfm] 

+ m%i-kj - Miikj-m)] , (211) 



where 



a n ,Ei- s , s \ . ^ 

5j-,„ + (-A_i)(„.<i,c) (212) 

For the particular directions fi = /i^q we single out the positive {hi = 1) and negative 
{hi = — 1) helicity for the charged lepton 71 1. 
In the helicity basis 

Qi-kj ± Mi{kyni) = [Ei^{\^) ± hiH][E^^ ^ /liq-k] (213) 
Therefore, for the pion at rest we obtain 

P,,{K = 1) = Mf{Ml - Mf) + 0{m]) (214) 

P,,{h, = -1) = mj ^^_-^, + 0{m^) (215) 
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Without approximations we obtain 



Pi -i f hi 



(216) 



while t he su m is given by Eq. (I209p . One can see this resuhs are in accordance with Eq. (2.38) 



of Ref . 



116| where the polarization [Pij{+) — Pij{—y\/[Pij{+) + Pij{—)\ is calculated. 



The ratio between the squared amplitudes for left-handed and right-handed helicities is 

P,[h, = -1) _ Mt 



R 



V 



Considering numerical values we obtain for Mj 

,2 



111 



R, 



(lOOkeV)^ 



X 4.92 X 10" 



while for M = M, 



R, 



rri: 



(leV)^ 



X 3.83 X 10^ 



(217) 



(218) 



(219) 



Considering the actual direct bounds for the neutrino masses 112| . we need a precision 
of 10~^ in the helicity measurement to reach those bounds either in the case of muons or 
electrons. Although the branching ratio to produce muons from pion decay is much larger 
than to produce electrons, the required dominant versus wrong helicity probability ratios 
are similar. 



Therefore, the coefficient 5 in Eq. f l202p has exactly the modulus 



If we rewrite 



Ml{M^-MlY 



rrij 



(220) 



(221) 



where E,y. = ^j^j ^ + 0(m|), we se e that \6^j \ is modified by the factor ^ when compared 
to the naive estimate mjj/2Ey 119| . We can also conclude that for the channel tt — )■ ez^ the 
real factor is enhanced considerably (~ 274x). 

In general, experiments can not achieve a perfect accurac y in helicity measurements 



because they usually involve polarization distributions |120| - |l22| . Even though the measure- 
ment of the helicity quantum number as a spin projection can only result in two discrete 
values, the projection direction is only defined within a finite accuracy. Thus we have to 
consider the accuracy necessary to be able to measure the wrong, non-dominant helicity 
amplitude. Parameterizing Eq. fl210p using fi-q = — cos^^ yields 



cos( 



(222) 



+ 



Ml + 2Mf 
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2 A/f2 ^ 



M2 - 

For cos^^ = 1 {cos 6 = —1) we recover the negative (positive) helicity for the charged lepton. 



Expanding around ^ = 0, we obtain 

M^-M?' 



sq2 j^A 

Pi,{5e) « MhMl -Mf)— + va] ^ . (224) 



Therefore we need an angular resolution of 

bQ = 2VR ~ 10^2 (225) 

to be able to probe the ratio R = Rij in Eq. fl217p . 

Considering the leptonic mixing the measurement of the wrong helicity for muons probes 

\M{7r ^ /.P : V = -1)P = |gP ^2_V2 «)eff ' (226) 

TT fj, 

where C = 2GFF.„Vud and 

(mljeS^Y.M'^h (227) 



is an effective mass for the muon neutrino, analogous to jll3 | inferred from the tritium 
beta decay experiments for the electron neutrino. The effective electron neutrino mass 
is defined as the expression in Eq. (12271) using |f/eiP instead of |f/^jp. In fact, can be 
extracted from tt^ — t- e^z/ by measuring the electron with negative helicity. 

To obtain the decay rate, we must multiply the amplitude squared by the factor 



1 1 



11/. M,^ 



47r 4M^ V Ml) 



/If 2 X 

) + 0(jnl) , (228) 



arising from the phase space. We then obtain 

^aiK = l) = ^F^\Vud\'M^Ml[l-^y+Oiml), (229) 

Tc.{K = -l) = ^F^lVudl'M^imlJ.s + 0{ml) , (230) 

where a = e, fi and {mljcs = The expr ession in Eq. ( I229p is the ordinary tree level 
decay rate for the pion decaying into /qZ/ fllq. 

Obviously, if new interactions involving the right-handed (chirality) components of neu- 
trino were present at low energy, the production of the wrong helicity states would be 
different. Therefore, by measuring the wrong helicity contrib ution it is also possible to 
constrain new interactions. Such exercise is carried out in Ref. 117| in the context of the 
two-Higgs-doublet model (2HDM). In addition to model dependence, the constraints are not 
very restrictive in the context of pion decay. 
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F. Flavor coupled to chiral oscillations 



The existence of chiral oscillations, studied in section [IVBt remind us that, when treating 
the time-evolution of the spinorial mass-eigenstate wave packets, we have overlooked an im- 
portant feature. We have completely disregarded the chiral nature of charged weak currents 
and the non-conservation in time of chirality. The flavor oscillation formula obtained in 
sections IIII Al and IIII CI were obtained by considering the oscillation of general flavor "par- 
ticle number" for general Dirac neutrinos. However, due to the left handed nature of weak 
interactions only left-handed components are produced and detected and then the flavor 
oscillation probability could be modified by the parallel effect of chiral oscillations. 

It is well known that from the Heisenberg equation, we can immediately determine 
whether or not a given observable is a constant of the motion. If neutrinos have mass, 
the operator 7^ does not commute with the mass-eigenstate Hamiltonian. This means that 
for massive neutrinos chirality is not a constant of the motion. We have already seen that 
localized states contain, in general, plane-wave components of negative and positive fre- 
quencies. As an immediate consequence, the interference between positive and negative 
frequencies, responsible for the additional oscillatory term in Dfo(<) flS^ . will also imply an 
oscillation in the average chirality. Thus, the use of Dirac equation as evolution equation for 
neutrino mass-eigenstate wave packets will lead to an oscillation formula containing both 
"flavor-appearance" (neutrinos of a flavor not present in the original source) and "chiral- 
disappearance" (neutrinos of wrong chirality) probabilities. 

As a simple treatment, we can use the evolution kernels of section UlI CI To incorporate 
the fact that neutrinos are created and detected as left-handed fermions into, for example, 
the conversion probability in Eq. (1761) . it is sufficient to use initial left-handed wave packets 
and replace the kernel of Eq. (1791) by the projected counterpart 



where ^^(p, t) = K^^{p, t)K^^{p, t) is the conversion kernel of Eq. ([79]) and L = (1 -75)72 
is the projector to left chirality. 

For an initial chirally left-handed Vf.-, i-e., ^/'i(a;,0) = il)2{x.,Q) = ipi,^{x) = Lil)y^[x), the 
flavor conversion probability can be explicitly written 



LK^}{p,t)LK'^,{p,t)L = ^^ip,t)L 




(231) 



■t) =V{u,^u^;t)-\ 
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The spinorial function ipuJip) is the Fourier transform of ipi,^{x). The conservation of total 
probabihty (182|) no longer holds because there is a probability loss due to the undetected 
right handed component 



1 



LKZ^RKl^,{p,t)L = -sm^2e 



nil 



1712 



sin £'2^1 L , 



(233) 



where R = {1 + 75)/2 is the projector to right chirality. We can see that the probability 
loss (I233p is proportional to the ratio rn^/E^ which is negligible for UR neutrinos. The 



probability loss of Eq. fl233p is ultimately induced by the time evolution of massive Dirac 
fermions which flips chirality from left to right according to 



En 



(234) 



The total probability loss for an initial left-handed electron neutrino turning into right- 
handed neutrinos, irrespective of the final flavor, is given by the kernel 

LKl^^RK^^{p,t)L + LKg^RK^^{p,t)L= cos^oi^^ s\n^ Eit + We(^^^ W E2t L. (235) 

To obtain the unphysical complementary kernels responsible for the conversion of right- 
handed component to right-handed and left-handed components, it is enough to make the 
substitution L ^ R in all formulas. 

We can interpret Eq. f l232p as a mixture of chiral and flavor oscillations. To that end 
we recall the formalism of section IIIIBt restricted to one spatial dimension, where the 
normalized mass-eigenstate wave functions ipj^ ^izjt) are created at time t = as a chiral 
eigenstate with eigenvalue —1. We can calculate 



Re 



dzil)\{z,t)^^il)j{z,t) 



dPz 2, 



Po) X 
1 - f{p^,m„mj) - 



miTrii 



E{p^,mi)E{p^,mj) 



cos [e _(p^, TO,, rnj)t] + 
cos[e+(p^, nii, nij) i]}(236) 



E{p^,m,)E{p^,mj) 

with i,j = 1,2. From the above integral, it is readily seen that an initial —1 chiral mass- 
eigenstate will evolve with time changing its chirality. Once we know the time evolution of 
the chiral operator, we are able to construct an effective oscillation probability which takes 
into account both flavor and chiral conversion effects, i.e., 



sin" (29) J 1 

9 1 9 2^ 



1=1 



dz\ilJi^L{z,t)\^ -Re 



1-75 



sin" (2^) 



[Dco(t) - Dfco(0] . 



tpu^iz^t-.e) 

-00 

dzij[^Liz,t) AA^^t) 

o 

(237) 
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As done before, the terms Dco(i) and DFCO(i) can be rewritten by using a p^-integration, 



1 ^ f'^°°dp 

Dco(i) = / -^^\p^-Po){l- c{p,,m,,m,)[l- cos[2E{p,,m,)t]]} 



(238) 



and 



/dpz 
i/{p^ -Po) { [1 - c{p^,m^,m2)] cos[e_(p^,mi,m2) t] + 
z vr 
—00 

c{p^,mi,m2) cos[e+{p^,mi,m2)t]} , (239) 



where 



c{p„mi,mj) = f{p^,m„mj) + 



2 E{p^,mi)E{p^,mj) 

The functions c{p^,mi,m,j) have a common maximum at = which, in opposition to what 
happens for f{pz,mi,m2), do not depend on the mass values, c^^^{0,mi,mj) = | and, following 
the same asymptotic behavior of f{p^,mi,m2), goes rapidly to zero for ^ rni^2- As a 
consequence of the first order approximation (159|) . we get 



C{p^,mi, rrij) 



l-V.V,+V.V,(v^+Vj-2)g^ 



I _ 2 P^^Po 

PO 



where we have considered the UR approximation in the second term. By substituting 
c{p^,mi,m^) in the above integrations (I238tl239|) . and after some algebraic manipulations, 
we explicitly obtain 



Dco(i) 



DFCO(t) 
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} } ,(240) 



2 Pq ^ 
PO 



2 2 



4pT 



COS 



2po 



i 7-5 9~Tr oin 

4p5 a^P^ 



2po 



} 



\ 2\/2apg 



2 2 
ii+m2 I 
,2 COS 
4Po I 



2 2 2 



2P0 



J -t sm 



a^Po 



Under the hypothesis of minimal decoherence between the mass-eigenstate wave packets 
{AvL <^ a), and for long distances between source and detector (L ^ a), i. e.. 



a Am? 



the standard flavor oscillation probability can be reobtained. In fact. 



sin^ {20) 
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sin^ {26) 



4Po' 



^ V2V2 
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(242) 
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G. Flavor conversion formulas for fermionic particles in an external magnetic field 



After obtaining the time-evolution of the spinorial mass-eigenstate wave packets in the 
presence of an external magnetic field in section IIVCI and the flavor conversion formula 
considering left-handed neutrinos in section IIVFI we intend to investigate here how flavor 
conversion is modified when chiral production/detection and external magnetic fields are 
taken into account. Consequently, the chiral nature of charged weak currents and the time- 
evolution of the chiral operator must be aggregated to the flavor oscillation formula. 

To have a complete description of the flavor conversion mechanism, taking chiral conver- 
sions into account, we must write the complete oscillation probability formula as 

V{u,,j^ ^ v^,j^-t) = ^^^^ { Dco(0 - DFCO(t) } , (243) 

where Dco(t) corresponds exclusively to chiral oscillations. Such term can be calculated by 
applying Eq. fll94p to each mass-eigenstate component yielding immediately 



Dco(0 = ^ I d'x 



= 1 (1 - MiW±MlW ) , (244) 

where ( )k refers to the average with respect to the mass-eigenstate WPs ipk, k = 1,2. 
From this point on we also turn back to the 3-dimensional analysis in order to not constrain 
the spatial configuration of magnetic fields. The average values of 75 can be, for instance, 
explicitly calculated in terms of the results of the Eqs. fllQSp and fllQSp . 
Analogously, the mixed flavor and chiral oscillation term can be given by 



DFCO(t) = ^ / d-^x 



1/4 {t,x) '^ ^'' i}2{t,x) + '4)\{t,x)^—^p^'4}^{t,x) 



(245) 



which deserves a more careful calculation. We shall see how we can explicitly construct 
the complete oscillation formula containing both "flavor-appearance" (neutrinos of a flavor 
not present in the original source) and "chiral-disappearance" (neutrinos of wrong chirality) 
probabilities for both of the particular cases of external magnetic fields discriminated in 
section IIVCI The following results are obtained after some straightforward but extensive 
mathematical manipulations where, again, we have imposed an initial constraint which 
establishes that the normalizable mass-eigenstate wave functions ipi^2{t,x) are created at 
time t = as a negative chiral eigenstate (14^1,275^^1,2 = —1)- All the subsequent calculations 
do not depend on the gamma matrix representation. 

Considering the first case of section IIVCI where the propagating momentum p is parallel 
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to the magnetic field B, we can calculate 



DFCO(t) = - 



t /I - (-i)=s-d^ 

w' I I W X 



cos 



Ei'^-Ei^A t]+ 1- 



cos 



, (246) 



where we have used the correspondence = ^{p ~ Pi)^{p - P2) and i^f'-* = y + mi*-' 

with i = 1,2 corresponding to the mass indices. In fact, in the UR limit, and for relatively 
weak magnetic fields (|a| ^ \p\), the contribution due to the very rapid oscillations does 
not introduce relevant modifications to the flavor conversion formula, in analogy to the case 
of purely chiral oscillations. On the other hand, by taking the NR limit, with a momentum 
distribution sharply peaked around \pi\ <^ rrii, the complete oscillation probability formula 
can be written as 



sin' (26*) 



-lyTi-a \ / r I'll 1 r (2\ 1^ ' 

w I COS [m). 't\ — COS [m\ 't\ 



(247) 



Equation fl247p introduces a completely different pattern for flavor /chiral oscillations, despite 
of not being phenomenologically verifiable. 

On the other hand, a more interesting interpretation is provided when we analyze the 
second case, where the propagating momentum p is orthogonal to the magnetic field B. The 
effects of the external magnetic field can be distinguished from the mass interference term in 
the flavor/chiral oscillation formula in the sense that we can write the complete expression 
for DFCO(t) as 



DFCO(t) = - 



(Pp 
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'(p){(cos[|aW|t] cos[|a(')|t] 
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cosii ej,''' 



to 
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t|cos|( la^'^l + |a(') 



t]]},(248) 



where we have used the correspondence of a^^^ and e*-"-* with m^''\ Again, the UR limit 
reduces the impact of the modifications to residual effects which are difficultly detectable 
by experiments. By taking the NR limit and following the same procedure for obtaining 
Eq. (I247p . the complete oscillation probability formula can be written as 



sin' (26>) 



cos' [m^^h] + cos' [m^^h] 

-2cos[(|aW| - |a(')|)t]cos[mWt]cos[m(')t]) . (249) 
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If we had chosen a wave packet exclusively composed by positive frequency plane-wave 
solutions, the high-frequency oscillation term would automatically vanish. Such a peculiar 
oscillating behavior is similar to the quoted ZBW and reinforces the argument that, when 
constructing Dirac wave packets, we cannot simply forget the contributions due to negative 
frequency components. 

Turning back to the foundations of the neutrino oscillation problem, we know that a more 
sophisticated field-theoretical treatment is required. Derivations of the oscillation formula 
resorting to field-theoretical methods are not very popular and the existing quantum field 
computations of the oscillation formula do not agree in all respects 181. The Blasone and 

n n 

Vitiello model |9|, 110| to neutrino/particle mixing and oscillations seems to be an interesting 
trying to this aim. They have attempt to define a Fock space of weak eigenstates to derive a 
nonperturbative oscillation formula. With Dirac wave packets, the flavor conversion formula 
can be reproduced |53| with the same mathematical structure. 

In the context where we have intended to explore the Dirac formalism, we have pointed 
out that a more satisfactory description for understanding chiral oscillations of fermionic 
(spin one-half) particles like neutrinos requires the use of the Dirac equation as evolution 
equation for the mass-eigenstates. Within such a framework, we have introduced the non- 
minimal coupling of the neutrino magnetic moment to an electromagnetic field in an external 
interacting process [40]. Although clear experimental evidences are still missing, we have 
reinforced the idea that the spinorial form and the interference between positive and negative 
frequency components of the mass-eigenstate wave packets can introduce small modifications 
to the standard conversion formulas when chiral oscillations are taken into account. More 
interestingly, when neutrinos interact with very strong magnetic fields, completely different 
oscillation patterns may emerge. 



V. INCLUSION OF QUANTIZED FIELDS 

This section considers the inclusion of some ingredients of quantum field theory (QFT) 
to the description of the oscillation phenomenon. There are divergent, and sometimes in- 
consistent, approaches in the literature. We try to identify some common aspects to our 
previous relativistic first-quantized treatments. 

Four aspects will be studied in detail: (i) Firstly, we want to analyze if considering 
quantum fields naturally eliminates the rapid oscillation effects coming from the interference 
between positive and negative frequency states, (ii) Second, we would like to know how EWP 
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approaches are related to the IWP approaches studied so far. (iii) Third, we want to relate 
the distinctive approach of constructing a flavor Fock space (9| with the IWP approaches 
considered, (iv) Fourth, we intend to investigate the details of the neutrino flavor state that 
is produced through decays when weak interactions, parent particle localization and finite 
decay width are considered. 

The first aspect (i) is considered in section IV Al when a simple second quantized descrip- 
tion of flavor oscillations for Dirac fermions is devised, based on free quantum field theory. 
It is commonplace that going from first quantized Dirac theory to the free quantum field 
theory of spin 1/2 fermions eliminates the interference between positive and negative fre- 
quency states by introducing a different notion of quantum states through the construction 
of the n— particle state space (Fock space) [71]. Such construction allows a meaningful inter- 
pretation of energy which is positive for both positive or negative frequency states and the 
notion of particle and antiparticle states becomes meaningful as well (for Dirac fermions). 
When interactions are considered through perturbation theory, the sectors that would cor- 
respond to positive and negative frequencies in first quantized theories remain separated by 
construction. It is then natural to ask what happens if we quantize the free fields associated 
to the mass-eigenstates that compose the flavor states. We will show that such procedure 
indeed results in naturally eliminating the rapid oscillation effects. Moreover, such a descrip- 
tion guarantees only particle or antiparticle propagation, keeping the nice property of giving 
normalized oscillation probabilities, like the first quantized examples treated in previous 
sections. We will also recover the initial flavor violation obtained in section IIIIEI without 
choice. 

To treat the aspect (ii), we compare in section IV Bl the different descriptions of neutrino 
propagation subjected to flavor mixing from the point of view of the propagators of the theo- 
ries. In first quantized theories, the role of the propagator is played by the evolution kernels 
introduced in section IIII CI The latter are compared to the mixed Feynman propagators 
used in the EWP approaches {s] to describe neutrino propagation. 

Aspect (iii) is treated in section IV CI where we briefly report about the Blasone-Vitiello 
approach to flavor mixing, giving emphasis to the mathematical similarities of the oscillation 
formula with the first quantized Dirac description of section IIII Al 

Finally, aspect (iv) is considered in section IV Dl where we calculate, within quantum field 
theory, the probability of detecting neutrino flavor states produced from pion decay including 
explicitly the weak interactions responsible for the decay, the pion localization and the pion 
finite decay width. Within this setting, it is possible to quantify the amount of intrinsic 
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flavor violation that arises when neutrinos are created and its relation with the decay width 
and pion localization. Such amount is still negligible but much larger than other indirect 
lepton flavor violation effects, such as — ?■ 67, where neutrino mixing only contributes at 
loop order. 



A. Simple second quantized theory 

Considering that only real neutrinos or antineutrinos (one of them exclusively) should 
travel from production to detection, the possibility to use the free second quantized theory 
for spin 1/2 fermions to describe flavor oscillations is now investigated. 

To accomplish the task of calculating oscillation probabilities in QFT we have to define 
the neutrino states that are produced and detected through weak interactions. Firstly, we 
define the shorthand for the combination of fields appearing in the weak effective charged- 
current Lagrangian (I273p 

i^aix) = Uaifiix) , a = e,iJL. (250) 

We will restrict the problem to two flavor families and use the matrix U as the same as in 
Eq. (l65l) . The mass eigenflelds z/j(x), i = 1,2, are the physical flelds for which the mass- 
eigenstates \vi{jp)) are well deflned asymptotic states. The free flelds Vi{x) can be expanded 
in terms of creation and annihilation operators (see appendix and the projection to the 
one-particle space deflnes the mass wave function 

^,,{x;gi) = {d\vi{x)\v,:g,) = ^ / d^P ^j^uf {x; p) , i = 1,2, (251) 

where 

Wi-9i) = d^P9t{p)Wi{P,s)) . (252) 

Since the creation operators for neutrinos (antineutrinos) can be written in terms of 
the free flelds z/j(a;) (z/j(x)), we can deflne the flavor states as the superpositions of mass- 
eigenstates 

\^a-{g}) = Uai\i^i:gi) . (253) 

The details of creation are encoded in the functions gi. 
We can also deflne 



'^'uo.uA^] {g}) = {o\i^e{x)Wa-{g}) 
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= UeiKii^uA^'i 9i) , 



(254) 



where ipui{x) are then mass wave functions defined in Eq. (125 ip . We can see from Eq. (I254p 
that if i)y^{x,t) = ilj^.^{x,t) = ipix), for a given time t, ilj^^^^{x,t) = ip{x) and ilj^^^^{x,t) = 
due to the unitarity of the mixing matrix. 

Although this approach does not rely on fiavor Fock spaces and Bogoliu 
tions, one considers the following observables to quantify fiavor oscillation 
charges, which are defined as 

Qa{t) = j d^x : vl^{x,t)va{x,t) : , a = e,n, (255) 

n 

where : : denotes normal ordering. Note that the Qe{t) + Qfj.{t) = Q is conserved [9], the two 
fiavor charges are compatible for equal times, i. e., [Qe{t), Qfj,{t)] = 0, and {u: {g}\Q\'y'. {g}) = 
±{u: {g}\i''- {g}) for any particle state (+) or antiparticle state (-). Notice that in the second 
quantized version the charges can acquire negative values, despite the fermion probability 
density in first quantization is a positive definite quantity. The conservation of total charge 
guarantees the conservation of total probability ( l82l) . 

With the fiavor charges defined, we can calculate, for example, the conversion probability 

V{u,^iy^;t) = {iye:{g}\Q,{t)We:{g}) (256) 
= U,,U;^U,,U:, j d?p e-^(^-^^)V; (p; g,)i^,^ (p; g^) , (257) 

where the neutrino wave functions ip^^ are defined in terms of the function gi{p) in Eq. (125 ip . 
Eq. (I257p is exactly equal to Eq. (IMp . if Eq. (p9|) is used with a = 9. If we could equate 
the two mass wavefunctions in momentum space ipuiiP', gi) = '4'v2{p\g2) = i^uSp) would 
obtain, from Eq. ( 1257p . the standard two family conversion probability (ITTl) . with |<^ve(p)P 
replaced by il)'l^{p)il)y^{p). However, the equality can not hold as proved in section lTlI E[ two 
wavefunctions with only positive energy components with respect to two bases characterized 
by different masses can not be equal. In fact, Eq. (I257p is equivalent to Eq. (l99l) with 
Eq. (I107p . and the states (I252p adopted in this second quantized formulation automatically 
accomplishes cond. B: there is no rapid oscillations. The fiavor conversion probability is 
given by Eq. ffTOOD with Eq. ffTOSD 

V{u,^u^;t) = isin220{l -ReE j d'p g^p, s)g2{p, s)e'^''^P^'[l - 2/2 (p)]} , (258) 

where gi{p, s) are the probability amplitudes in momentum p and spin s. One can see that 
the probability (I258p is never null and one can not avoid initial fiavor violation but one can 
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still adopt equal momentum distributions as in Eq. fll2ip : gi{p,s) = g2{p,s) = g{p,s). In 
this case we have initial flavor violation given by Eq. (I122p : 

V{i^^^i^^,;0) = sm^29 J fiipMp)? , (259) 

where = \9{Pi^)\ ^^e spin-independent momentum distribution. The conver- 

sion probability f|258p yields 

V{u,^v^-t) = \s\T?2e[l- jd^p\g{p)\^cos[AE{p)t][l-2f2{p)]] . (260) 

An approximate version of the formulas (I259p and ( I260p was previously presented in Ref. jsjj. 

The conversion probability for antineutrinos z/g — )■ z/^ is exactly the same as Eq. fl258p by 
using Eq. (lA29l) . 

The formulas obtained in this second quantized version do exhibit the interference terms 
between positive and negative energies like in Eq. (ISTl) . Such interference terms are ab- 
sent because the possible mixed terms like b2{p)a\{p)\0) are null for an initial flavor state 
superposition that contains only particle states (or only antiparticles states). 



B. Connection with the EWP approach 

Now we try to establish a correspondence between our results and the quantum fleld 
theory (QFT) treatments, notably the EWP approaches. Thes e ap proaches follow the idea 
that the oscillating particle cannot be treated in isolation [ll|, IfioiEl, i.e., the oscillation 
process must be considered globally: the oscillating states become intermediate states, not 
directly observed, which propagate between a source (A) and a detector {B). This idea can 
be implemented in QFT when the intermediate oscillating states are represented by internal 
lines of Feynman diagrams and the interacting particles at source/detector are described by 
external wave packets [s], [gqI . 

Let us consider the weak flavor-changing processes occurring through the intermediate 
propagation of a neutrino, 

A + la ^ A' + (oscillation) + B ^ Ip + B' , (261) 
A ^ A + + (oscillation) + B ^ Ip + B' , (262) 

where A and A' {B and B') are respectively the initial and flnal production (detection) 
particles. The process (1261 p may describe the process jlll . Isol where a charged lepton 
hit a nucleus A turning it into another nucleus A' with emission of a neutrino (this process 
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happens around xa = (^a, ^Ea))- The process fl262p describes the similar process where 
particle A decays into A' with emission of laiVa'-, processes such as pion decay, A' is absent. 
Subsequently the neutrino travels a long distance and hit a nucleus B which transforms into 
B' emitting a lepton (this process happens around xb = (tB,XB))- The whole processes 
for Eqs. (IM and (12621) look hke 

(L + A) + (B) ^ {A') + {B' + Ip) (263) 
{A) + {B) ^ (L + A') + {B' + /^) (264) 

with transition amplitudes given generically by 

A/3 = {A',B\Ip\S\A{xa),L{xa),B{xb)) . (265) 
Aap = {A',B',l0,L\S\AixA),B{xB)), (266) 



where S is the scattering matrix. The initial states are localized |8| while the final states 



might be localized states [8] or momentum eigenstates |59|. 

After some mathematical manipulations jsf, both amplitudes in Eqs. 0265 p and f l266p can 
be represented by the integral 

^B{k)Gi3aik)'^Aik) exp [-ik-ixB - xa)] , (267) 

where the spinorial function represents the overlap of the incoming wave packets at the 

source while ^b(^) represents the overlap of outgoing wave packets at the detector. Notice 
that Eq. (1267p is an adaptation for propagating Dirac fermions of the formula obtained 
in Ref. 8| for scalar neutrinos and scalar interactions. The Green function in momentum 
space, Gi3a{k), represents the mixed neutrino propagator which carries the information of 
the oscillation process: 

G(Saik) = Y,Uf3jG,ik)ul, (268) 

j 

where Gj{k) represents the propagators for each neutrino Uj. Thus Eq. (126 7p can be also 
rewritten as the superposition of amplitudes, 

Apc. = Y.^PjAjUl. (269) 
j 

The overlap functions are independent of production/detection times {tA,tB}, and produc- 
tion/detection positions {xa,Xb}, but depends on the directions of incoming and outgoing 
momenta. In a certain way, the physical conditions of source and detector, in terms of 
time and space intervals, are better defined in this framework than in the intermediate wave 
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packet framework. Anyway, to understand the oscillation process we must turn back to the 
definition of mixing in quantum mechanics. It is similar in field theory, except that it applies 
to fields, not to physical states. This difference allows us to bypass the problems arising in 
the definition of fiavor and mass bases jsl- 

Let us analyze some Green functions Gap{k). The main improvement of the covari- 
ant approaches developed in section IIIII is that the propagation kernels governed by Dirac 
and Sakata-Taketani Hamiltonians are cau sal, i.e., are null for spacelike separations (see 
Eqs. (lAlSP and (IA19P and Refs. (tiI Q, 124 |). On the contrary, the kernel of spinless 
particles restricted only to positive energies is not null for spacelike intervals [84]. From 
the point of view of relativistic classical field theories, a causal kernel guarantees, by the 
Cauchy theorem, the c ausal connection between the wave-function in two spacelike surfaces 



at different times 



124|. 



To compare the IWP and EWP approaches it is useful to rewrite the Dirac evolution 
kernel for a fermion of mass m^, present in Eq. (1731) of section llTTt in the form 71 1 



Kf{x 



y) = E 



2En 



1,2, 



(270) 



where (x — y)^ = t,{x — yY = {x — x'Y when compared to the notation of Eq. ( 1731) . The 
spinorial functions u, v, are the free solutions of the Dirac equation and they are explicitly 
defined in appendix |X1 (More familiar forms for the function S are also shown in appendix 
lAl ) Clearly the function ^^(x—y; mj) = (0|{'?/'j(x), '?/^j(?/)}|0) satisfies the homogeneous Dirac 
equation with mass rrij ( l66l) and it is known to be null for spacelike intervals {x — y)"^ < 
ensuring micro causality 84j. Il24| . 

In contrast, the Feynman propagator iSpix — y) appears in QFT. It is a Green function 
for the inhomogeneous Dirac equation obeying particular boundary conditions. The EWP 
approaches use this Green function for the propagation of virtual neutrinos. In fact we have 
for the Green function in Eq. (12681) . 



Gj{k) = iSF{k;mj 



rrij + ie 



(271) 



To directly compare the Feynman propagator to the kernel in Eq. (I270p we can write iSp in 
the form 



iSp^x — y; rrij) 



E / ^ [Uj{x;p)u'j{y;p)e{xo - yo) 
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-v'j{x-p)v'^{y-p)e{y,-x,)] . (272) 

Although the function Sp is called causal propagator, it is not null for spacelike intervals, 
and it naturally arises in QFT when interactions are present and treated in a covariant 
fashion. Equation fl272p shows that the propagator Sp describes positive energy states 



propagating forward in time and negative energy states propagating backward in time |7l| . 
At a first glance, both neutrino and antineutrino parts of Eq. ( I272p seem to contribute to the 
space-time integrations present in covariant perturbation theory, as neutrino-antineutrino 
contributions in Eq. (I270p have led to Eq. (IHTl) . 

However, it was shown that in an EWP approach, only on-shell contributions contribute 
for large separations between production and detection (see the review in Ref. [8i]). Moreover, 
although both neutrino and antineutrino parts may contribute as intermediate neutrinos for 
certain situations js^, the wrong state contributions are very much suppressed. Therefore, 
the intermediate neutrino states propagating in the EWP approach are dominated by real 
and on-shell neutrino (antineutrino) states. 

Let us consider the process ( 12611)^ and analyze more carefully some calculations using, 
instead of the scalar interaction sl, Isojl, the effective charged- current weak Lagrangian 

Af=3 

= [la{x)-f^'LUaiiyiix)Jf,{x) 

i,a=l 

+ u,{x)U:,YLUx)Jlix)] (273) 
= ^1 + Sfl , (274) 

where G = \/2Gf and is the sum of any effective leptonic or hadronic current. The 
Lagrangian f l273p is written only in terms of physical mass-eigenstate fields, which coincides 
with flavor state fields only for the charged leptons: li{x) = e(x), hix) = fi{x), .... 

The lowest order nonzero contribution of the scattering matrix 5* to Eq. fl265p is second 
order in the Lagrangian ( 12 73 p . More explicitly, the term that contributes to the amplitude 
f l265p comes from 

^(2) ^ ^t(^vk)' = -^Ti^i + S£lf (275) 
~ -T(^i)(^it) (276) 
~ I dWy5^^^„(x,y) , (277) 

where ( ) stands for space-time integration and 

^^a{x,y) = ^:Mx)l^{x)-/^'LU^itSF{x - y;m,)U:a"Lla{y)4{y): ■ (278) 
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In Eq. fl276p we kept only the mixed product and in Eq. fl277|) we kept from all possible 
terms in Wick expansion [7l| only the term responsible for the transition of interest. 

We can write in explicit form the production and detection spinorial wave functions in 
Eq. ([267D asjsl 

^Bik) =iG j d^ye-'^<y-^-^ {B' ,lp\J,{y)lp{y)^^L\B) , (279) 
^^(A;) = iG [d^xe'^<''-''^^A'\YLLix)Jl{x)\A,Q . (280) 



Since the intermediate neutrinos in the process (I26ip can be considered real on-shell 
neutrinos in the EWP approach for large separations between production and detection, 
both processes in x^ and Xb should be considered as real scattering processes with real 
neutrinos involved. These informations permit us to rewrite Eq. fl269p in a slightly different 
form 

Y.Up^U:,A. = E / ^ / d'y{B',l,\^,{y)e^(^-^^y^-\B,Mp)) 

J d^xOiy - x){A',u,{p)\^h^)e'''-''^\A,L) , (281) 

where P = {H, P) is the energy-momentum operator. A change of notation were made here: 
in Eq. (128 ip the states \B) and 1^4, l^) are centered around the origin while in Eq. fl265p they 
are respectively centered around Xb and xa', the translation is explicitly performed by the 
translation operator e*^'^. Additionally, the step function 6{y — x) is necessary to ensure 
that the contributions of points y around Xb should always be after the contributions of 
points X around xa- 

Equation (128 ip shows us the amplitude of the entire process from production to detection 
in "decomposed" form (apart from the step function in time): the amplitude of production 
process multiplied by the amplitude of detection process summed over all possible inter- 
mediate real neutrinos of different masses and momentum p. (The sum over spins are 
implicit.) Thus the EWP approach is not conceptually different from IWP approaches when 
the three amplitudes of production, propagation and detection of neutrinos can be factored 
out independently from each other p, [sl, l56|. For those cases, for most of the situations, 
IWP approaches provide the same results as EWP approaches if the localization aspects 
can be transposed as inputs to neutrino wave packets. A situation where the decomposition 
can not be performed simply and the localization aspects have to be explicitly taken into 
account is exemplified by the unusu al ca se of Mossbauer neutrinos recoilessly emitted and 



detected from bound state electrons 



125| 
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C. The bridge to the quantum mixing 



In one- dimensional spatial coordinates, the mixing is illustrated by the unitary transfor- 
mation 

^„{z,t;9) = g-'{e;t)^i{z,t)g{e;t) (282) 

as the result of the nonco incidence of the flavor basis (cr = a, P) and the mass basis (z = 
1, 2). The Eq. fl282p gives Eq. when the generator of mixing transformations Q{6]t) is 
given by 

gie;t) = exp{e j dz[ip^{z,t)i)^{z,t) - i^^{z,t)iP,{z,t)]]. (283) 

By taking the one-dimensional representation of Eq. f l267p . the propagator G{E,p^,tB,tA) 
can also be written in the flavor basis as 

G^^{e-E,p^,T) = g-\e-t)GiE,p,,T)g{e;t) 

= g'\e-t)G{E,p^,tD,tp)g{e-t) (284) 

with T = ts — tji- 

In particular, the definition of a Fock space of weak eigenstate s be comes possible and 
a nonperturbative flavor oscillation amplitude can be derived jo], Il23 |. In this case, the 
complete Lagrangian (density) is split in a propagation Lagrangian, 

Cp{z,t) = ipi{z,t) {i^ -m-i) ipi{z,t) + ipiiz^t) {i^ - m^) 'ip2{z,t), (285) 

and an interaction Lagrangian 

C,{z,t) = 'ipa{z,t;e) {i^ -mo) ipa{z,t;6) + ■ipi3{z,t;0) {i^-rrii^) ipi3{z,t;e) 

- map {Mz, t; 0)ipp{z, t; 6) + ^Pp{z, t; 0)^„(z, t; 9)) , (286) 

where 

= "^1(2) cos^ 6 + m2(i) sin^ 6 

and 

mai3 = {mi — mj) cos sin 6*. 

In general, the two subsets of the Lagrangian can be distinguished if there is a flavor trans- 
formation which is a symmetry of Ci{z,t) but not of Cp{z,t). Particle mixing occurs if the 
propagator built from Cp{z,t), and representing the creation of a particle of flavor a at point 
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z and the annihilation of a particle of flavor /3 at point z\ is not diagonal, i. e. not zero for 
/3 = a. The free fields ipi{z,t) can be quantized in the usual way by rewriting the momentum 
distributions b^{p^,m.i) and d^*{-p^,mi) in Eq. (15T]) as creation and annihilation operators 
B^{p^,mi) and D^\-p^,m,). The interacting fields are then given by 



(287) 



s=l,2 



where the new flavor creation and annihilation operators which satisfy canonical anticom- 

1 1 

mutation relations are defined by means of Bogoliubov transformations |123] as 

B^{p^;t) = g~\e-t)B'{p^,m,)g{e-t) 
D',{-p^-t) = g-\e-t)D\-p^,m,)g{6-t). 

By following jo], Il26 |. the fiavor conversion formula can be written as 

V{vc.^up;t) = i?^(po;i)}r + |{^|(-Po;i), 5^bo;i)}r (288) 

which is calculated without considering the localization conditions imposed by wave packets, 
i.e., by assuming that ^ po. When the explicit form of the fiavor an nihil ation and 
creation operators are substituted in Eq. f l288p . it was also demonstrated jl23 | that the 
fiavor oscillation formula becomes 
sin" (2^) 



2 "{[1 - /bo,"ii,2)] cos[e_(po,"ii,2)i] + /(po,™i,2) cos[e+(po,TOi,2) *]} 
sin" (20){[l-(|^Vl TAm^ 



sm 



4po 



+ 



/ AmV 



V 2po / 



Sin 



Pot 1 + 



m\+m2 



} (?89) 



where the last approximation takes place in the relativistic limit p^ ^ After 
some simple mathematical manipulations, Eq. fl289p gives exactly the oscillation probability 
'Poiracii^a^ 1^13', L) calculated from Eq. ( !62|) when it is assumed that the wave packet width 
a tends to infinity and t ^ L. 

This new oscillation formula tends to the standard one (HT]) in the UR limit. If the mass- 
eigenstates were nearly degenerate, we could have focused on the case of a NR oscillating 
particle having very distinct mass-eigenstates. Under these conditions, the quantum theory 
of measurement says that interference vanishes. Therefore, as we have already pointed out, 
the effects are, under realistic conditions, far from observable. Besides, in spite of working on 
a QFT framework, the lack of observability conditions must be overcome by implementing 
external wave packets, i.e., by calculating the explicit form of Eq. fl267p for fermions. Such 
a procedure was applied by Beuthe for scalar particles jsl and, in a very particular analysis, 
based on our intermediate wave packet results, it could be extended to the fermionic case. 
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D. Intrinsic flavor violation 



In this subsection we will try to quantify the initial flavor violation f ll20p . i. e., the amount 
by which neutrino flavor is not well defined at creation or detection. To accomplish this 
task, it is crucial to analyze two aspects: (i) the definition of neutrino flavor and (ii) the 
interactions responsible for neutrino creation and detection. We will focus on the second 
aspect and justify the use of the usual neutrino flavor states 127 1 

Mp)) = U*Jui{p)) , \v^{p)) = UaiMp)) , (290) 

where {i>i{p)\i>j{p')) = 6ijS^{p — p'), as approximately defining a neutrino (antineutrino) 
state of momentum p and flavor a, associated to charged lepton /q. 

We will treat explicitly here the antineutrinos created through pion decay: tt^ — lai^a, 
a = e, /i. Although, many of the conclusions drawn within this process may be extended 
to more general contexts. The two processes a = e, fi should be considered as different 
superpositions of the six channels tt^ — )■ kuj, z = 1,2, j = 1,2,3, where li = le = e,l2 = 
In = fi. These six channels contribute slightly differently for different neutrinos and for each 
channel with fixed la, the creation of a pure neutrino flavor state Ua, as defined in Eq. ( I290p . 
is not guaranteed. Indeed, we will show that processes such as tt — )■ fiUf., are possible with 
a branching ratio much greater than loop induced processes such as /i — )■ 67 in the SM, 
considering the known massive neutrinos and leptonic mixing. 

Let us define the free two-particle states with definite flavor 

IIMM'^)) = SaiU^j\liiq)ujik)) . (291) 

The charged lepton states remain as mass-eigenstates while the neutrino states are mixed 
through Ui3j. We will see, in accordance to usual expectations, that pions decay mainly into 
the states |ZQ,(q)i/^(k)) with [a, /3) = {fi, fi). However, we will also see that there is a non-null 
probability of the pion to decay into the neutrino flavor violating states with [a, /3) = (/i, e) 
or {a,/3) = (e,/x). For that purpose, we want to ultimately calculate the probability 

V{7T ^ lau^;t) = jd\j d\ \{U^)M^)Ht))\^ , (292) 

spins 

where |7r(t)) describes the pion state that decays as time t evolves. The initial condition is 

k(0)) = k^) = J d^pHp)Hp)) , (293) 

where |vr(p)) is the free pion state normalized to (7r(p)|7r(p')) = S^{p — p'). This normaliza- 
tion will be used throughout this subsection for free one-particle states. The function ip{p) 
characterizes the initial momentum distribution for the pion. 
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The calculation of |7r(t)) can be performed at lowest order by using time dep endent 



69 



128 



The 



perturbation theory and the approximation method of Wigner-Weisskopf 
details can be found in appendix [B] and Ref. The important result is summarized by 
the time dependent transition amplitude Xii(q, k;^) for \7i^) — )■ |/j(q)i/j(k)), 



Xij(q,k;t) = Xij(p,q,k;t)V^(p) 



p=q+k ' 



-1/2 f^ijjp, q, k) 



(294) 
(295) 



where Nij = {27rf2EiXq)2E^^{k)2E^{p), AEij = E^ - Ei^ ~ ^•^^^'^ is the pion decay width 
and = ^ij{p^(\^'k) = ^{7r^{p) ^ li{q)uj{'k)) is the usual invariant matrix element, 
normalized to obey {'n-{p)\S\li{q)uj(k)) = (27r)^5^(j9 — q — k){—i)^ij. 

), for t ^ 1/r, we can rewrite Eq. fl292p as 

2 

q=p-k 



By using Eqs. fl294p and 



where 



spins j 

^1/2 ^aj{P, q, k) 



(296) 



C^aj-^aj(P,q,k) = 



AE^j - I— 



(297) 



We see the exponential e"*^"^* is responsible for the neutrino oscillation phenomenon. In 
fact, if we neglect the neutrino mass mj in every term of Eq. (12961) . except in the exponential, 

we get 

V{t: ^ Ic^up-t) = Jcfp\i;{p)\''Jd\V{i^a^i^p;t)\F^W (298) 

where F„ = (F„j)„^,^o- 

Notice the usual oscillation probability. 



(299) 



factors out from the creation probability of laV, l-^aP? for massless neutrinos. Such factor- 
ization is what allows the definition of the state (129 ip as a flavor state, since 



(300) 



for 1/r ^ t ^ Losci where Lose is the typical flavor oscillation length (period). The refore, th e 



antineutrino flavor state Ucij\vj) is only created jointly with the charged lepton | 129 
Notice Eq. fl3UUp correctly coincides with the branching ratio of the decay tt — t- 
For the sake of completeness let us rewrite the relevant part of Eq. (I300p as 



130|. 



d\\FJ 



lq=p-k 



i Vq,j I iAJjQ.j ''27 I 



q=p— k,mj=0 



r ' 



(301) 
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where the last approximation considers \AEaj — "^^P ~ ^Y^^i^^aj) for small enough T. A 
regularization function / might be necessary to guarantee the convergence of the integral 
for large k 63|. 

Neutrinos, however, are not strictly massless and we may have initial flavor viol ation 
because different neutrino masses contribute differently to each channel vr — > liUj jl29|. We 
will focus on initial flavor violation and denote the interval of time satisfying l/F ^ t ^ Lose 
by t = 0. 

We can make the flavor violating contributions explicit by rewriting the term inside the 
squared modulus in Eq. (12961) as 

3 3 

+ Y,UajU;jAFo,j , (302) 



i=2 



where AFaj = F^j — F^i. Thus the squared modulus becomes 



(303) 



i=2 j=2 

Notice there is no summation over repeated indices a or /3. We recognize that only the last 
term of Eq. (I303p is flavor non-diagonal. The second term, which is flavor diagonal, can be 
shown to be of the same order of the flavor violating contribution, but negative in sign. 

Specializing to a 7^ /3, under the approximation of UasU^^ ~ (which is valid if a = e or 
(3 = e), the initial creation probability yields 



n^u, = Jd^Pmp)\^jd\\Uo^2U;2\^\AF^2\ 



(304) 



where Vi^^b = ~^ l-a^f3] t), for l/F ^ t ^ Lose- Fot the two family parametrization, we 
have |t/a2f^^2p ~ i sin^2^, thus indicating that this phenomenon is indeed mixing dependent. 

To determine the most dominant contribution to Eq. (I304p . we have to analyze the 
dominant contribution to AFa2 due to Am = mi — m2 7^ 0. We anticipate that the 
dominant contribution is due to {AEaj — ir/2'~f)~^ in Faj. Therefore 



a2\ 



where 



aOl 



aO 



AE^2-^^ AE„i-2^ 
(AE,) 



aO 



{AE^o + \AE,? + $ 



AE, = E,,(k)-E,,(k) 



{AE^o - \AE,f + § 



Am^ 



£;^,(k) + £;^2(k) 



(305) 
(306) 

(307) 
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the subscript means we assume the massless hmit, mi, m2 — )■ 0, for neutrinos in those terms 
and l^ooP refers to \^aj\mj-^o without the mixing matrix element [see Eq. f|297p ]. 

We have also the equivalence |^qoP = J2j l-^ailmj^o- 

To justify Eq. (1305 P we note that F ^ \E^-^{]<.) — E,y^(k)\ = \AEa2 — AEai\ for neutrino 
momentum close to the energy conserving value |k|,^ = E^ = ^2M^" ' massless neutrinos 
and pion at rest p ^ 0. Retaining small neutrino masses does not alter the conclusion. 
Neutrino energy (momentum) for Ui and 1/2 can be distinct from the energy conserving 
values only by an amount T E^ and such variation are not relevant in other terms besides 
the last one in Eq. (1305^ . The contribution of neutrino masses to the invariant matrix 
elements \^aj \ and other kinematical terms are negligible since neutrinos are produced UR. 
Neutrino masses can be neglected in most of the terms except for 

^E^^^. (308) 



We recall the numeric values for the relevant quantities F = 2.53 x 10 ^eV 112| and 



l^rn?/2Ey ~ ^ x lO-'^eVff^, where Am^ is either lAm^gl ~ 0.8 x lO^'^eV^ or |Am 



2 I 

23 1 



2.5 X lO^^eV^ jl31 |. A more detailed justification can be found in Ref. j63 1 . 

The flavor violating creation probability in Eq. (I304p . a 7^ /3, can be calculated by chang- 
ing the neutrino momentum variable |k| to energy Ey and extending the lower integration 
limit to —00, with negligible contribution, giving 

1 n / Am2\2 



where the two family parametrization, |?7a2f^^2p ~ \svo?26, was employed and p ^ 
(pion at rest) was considered by adjusting ip{p). The subscript EC denotes that energy 
conservation was assumed. The following integral was also necessary, 

1 _27r/272 

47' 

One can recognize the term inside parenthesis in Eqs. ( 1306P and (1310p as the additional 
contribution that appears in Eq. fl309p . 

Let us estimate some specific fiavor violation probabilities (branching ratios): 

sin220i2 sin220i2 F ' sin22023 

To compute the last value in Eq. (I31ip . we considered lAm^gl ^ |Am23| ^ [Am^gl ^.nd 
replaced AF„3 by AF«2 in Eq. (IHOD . 
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For completeness, let us consider non-realistic cases that do not satisfy ^ 1. In that 
case Eq. ( I309p should be corrected to 

Vi^^^^lsm'ie-^^^, (312) 

where 



B = 



2E^T 



(313) 

EC 

The expression in Eq. 031 2p is obtained if we retain the terms in the denominator of 
fl306p and consider its energy conserving value B. In particular, Eq. fl312p have the correct 
limit for 

which corresponds to the incoherent creation limit. This limit can be also seen in Eq. (I296p . 
considering that Eq. (I297P for such limit implies there is no overlap in the sum. Equation 
fl296p could be rewritten 

^ / d'pmp)\' [ d\E\UajU^,/ E \Faj\l^^_^^^ , (315) 
J J i spins 

where Ej=i l^^ajf^j^P = | sin^26' in the two family approximation. In other words, the 
condition ^ 1 corresponds to the limit where the energy uncertainty is much smaller than 
the energy separation of mass-eigenstates, which allows the distinction of mass-eigenstates 
and consequent loss of coherence for oscillation. We can see a close relationship between 
intrinsic flavor violation and loss of coherence by wave packet separation: if intrinsic flavor 
violation is large, flavor oscillation is also suppressed. 

We should remark we do not analyze the effects to oscillation due to flnite size in detec- 



tion 



57| but if the microscopic time of detection for any interaction is much shorter than 



r = 1/r, then intrinsic flavor violation should remain. The intrinsic uncertainty in momen- 
tum CTp of the parent pion, encoded in ip{p), was also not considered in the discussion but 
the conclusions remain unchanged as long as dp < cr^ = F. In fact the limit |'?/'(p)p — )■ S^jp) 



is easily calculable. Other subtleties of the approximation used can be seen in Ref. [63 1. 
The same conclusion of requirement of coherent superpositions of ne utrin os can be drawn 
in an exactly solvable QFT that does not rely on perturbation theory 129| and then all the 
aspects of flavor oscillation and flavor violation can be calculated exactly. It is important to 
emphasize that the effect of intrinsic flavor violation is more robust than flavor oscillation 
effects because no propagation is needed and averaging over the source or detector does not 
affect the results. 
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It should be also emphasized that the intrinsic flavor violation effect calculated in 
Eq. (1309 p assumes two facts: (1) neutrinos are directly detected and (2) they are detected 
as flavor states, which are coherent superpositions of mass-eigenstates, as deflned (approx- 
imately) by Eq. (12911) . On the one hand, the coherent creation of neutrino flavor states 
is indeed guaranteed from the observations of neut r ino o scillations, implying that intrinsic 



flavor violation effects should be small 



48 



63 



129 



130j. If mass-eigenstates were created 



and detected incoherently, flavor violating effects would be analogous to flavor changing 
processes for quarks, at tree level, without the explicit appearance of the Am^ dependence. 
On the other hand, when neutrinos are not explicitly detected, their effects ca n be computed 
from an incoherent sum of the contributions of each neutrino mass-eigenstate 115| . as in the 
intended direct measurements of absolute neutrino mass. Another aspect of ass umpt ion (2) 
is that the definition of the usual flavor states in Eq. (129 ip is only approximate 130] and it 
may depend on the details of the creation and det ectio n. Usually, however, the error made 
is of the order of lS.vp? jE^ and it can be neglected 130| . 



VI. CONCLUSIONS AND OUTLOOK 



In order to understand some subtle changes that appear in the standard flavor oscillation 
probability 50|| when one goes deep in the theoretical fundamentals of the quantum oscil- 
lation phenomenon, we have quantifled the consequences of assuming the quantum mixing 
and oscillation in the Dirac theory. 

We have started our study by analyzing the consequences of including the spatial local- 
ization in the standard effective treatment of quantum flavor oscillations. The classiflcation 
of flrst and second order effects were relevant in establishing a satisfactory criterion for quan- 
tifying the effects of including spin and relativistic completeness for fermionic particles in 
the Dirac theory. By taking into account the spinorial form of neutrino wave functions and 
imposing an initial constraint where a •pure flavor state is created at t = 0, it was possible, 
for a constant spinor ly, to obtain the oscillation probability containing the contribution of 
positive and negative frequency solutions of the Dirac equation. We have noticed a term 
of very high oscillation frequency depending on the sum of energies in the new oscillation 
probability formula (rapid oscillations). 

In addition to rapid oscillations, the spinorial form of the wave functions and their chiral 
oscillating character subtly modify the coefficients of the oscillating terms in the flavor 
conversion formula. To describe the time evolution of the mass-eigenstates, we have assumed 
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an initial gaussian localization and performed integrations by considering a strictly peaked 
momentum distribution. Under the particular assumption of UR particles, we have been 
able to obtain an analytic expression for the coupled chiral and flavor conversion formula. 
For the case of Dirac wave packets, these modifications introduce corrective factors which 
are negligible in the UR limit. Nevertheless, strictly from the theoretical point of view, we 
have confirmed that the fermionic character of the particles modify the standard oscillation 
probability which was previously obtained by implicitly assuming a sca/ar nature of the mass- 
eigenstates, restricted to positive frequencies. These results concerning chiral oscillations 



coupled with the flavor conversion mechanism of free propagating wave packets jSll |52| . |54| . 
are therefore relevant in quantifying the role played by the intrinsic spinorial structure in 
flavor oscillations. 

The quantum-mechanical treatment which associates Dirac wave packets with the propa- 
gating mass-eigenstates is indeed rich in physical insights which were extensively discussed. 
Besides the review of analj^ical calculations done with gaussian wave packets for scalar and 
fermionic particles, the main conceptual aspects related to the very rapid oscillations lead to 
the study of chiral oscillations. In the standard model flavor-changing interactions, neutri- 
nos with positive chirality are decoupled from the neutrino absorbing charged weak currents. 
A state with left-handed helicity can be approximated by a state with negative chirality in 
the UR limit. Once we have assumed the interactions at the source and detector are chi- 
ral, only the component with negative chirality contributes to the propagation. Since the 
chiral state affects the results at the detection process, the mechanism of chiral oscillation 
results in modifications to the standard flavor conversion formula, with correction factors 
proportional to m^^/Po- '^^^ effect of such factors are, however, practically undetectable 
by the current experimental analysis. It leads to the conclusion that, in spite of often be- 
ing criticized, the standard flavor oscillation formula resorting to the plane wave derivation 
can be reconsidered when properly interpreted, but a satisfactory description of fermionic 
(spin one-half) particles requires the use of the Dirac equation as evolution equation for the 
mass-eigenstates. Correlated constructions involving more general spatial configurations for 
the magnetic fields might be considered as well as extensions of the corrections discussed 
here to more general chiral conversion rates, polarization effects and neutrino propagation 
in magnetized media. Some of these investigations has already been under consideration, 
primarily focused on observable effects or phenomenological implications jl?!. Ill7 |. Another 
difficult task in the same framework is to describe the nonstationary evolution of neutrinos 
in supernovae or in the early universe, where interaction rates are in competition with the 
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flavor oscillation frequencies, and where the wave packet effects seems to be maximized. 

In parallel to the description in various contexts of the phenomenon of neutrino flavor os- 
cillation based on the Dirac theory of fermion propagation (sections lllll and HVT) . a dichotomy 
between rapid oscillations and initial flavor violation was introduced. Initially discussed in 
the context of the free Dirac theory, we could conclude that at least one of the phenom- 
ena of rapid oscillations or initial flavor violation should invariably occur in neutrino flavor 
oscillations. Rapid oscillations could be avoided but only at the expense of not having an 
initially well defined flavor. In fact, as discussed in subsection IIIBt initial flavor violation 
could arise even in the context of the usual positive-energy scalar wave packet description 
of neutrinos, reviewed in section [Tll However, in the latter context, it is usually eliminated 
completely, by setting 0i(a;,O) = 02(ic,O), or approximately, by taking |0i(a3,O) — 02(a^,O)| 
small enough. Using Dirac theory, however, equal wave packets for the mass-eigenstates, 
^pl{x,0) = tp2{x,0), imply rapid oscillations. Quantitatively, both phenomena were unim- 
portant because their effects could be as small as of the order of (Am/E'j,)^. Moreover, in 
the context of first quantized Dirac theory, there are no sufficient ingredients to decide which 
effect is taking place in neutrino oscillations and with which size. 

Seeking a more sophisticated approach towards more realistic descriptions, we have con- 
sidered the inclusion of field-theoretical elements in section |Vl Derivations of the oscillation 
formula resorting to field-theoretical methods, however, are not very popular. They are 
thought to be very complicated and the existing quantum field computations of the oscil- 



lation formula do not agree in all respects [8f. For instance, in the Blasone and Vitiello 



model 9|, llO[ they have attempted to define a Fock space of weak eigenstates and to derive a 
nonperturbative oscillation formula. Flavor creation and annihilation operators, satisfying 
canonical (anti)commutation relations, are defined by means of Bogoliubov transformations. 
As a result, new oscillation formulas were obtained for fermions and bosons, with the oscil- 
lation frequency depending not only on the difference but also on the sum of the energies 
of the different mass-eigenstates. The BV model, however, is not unanimously accepted 
throughout the community. Among the arguments against it, it has been argued that Fock 
states of flavor neutrinos are unphysical [12] because their use to describe the neutrinos 
produced or detected in charged- current weak interaction processes imply that measurable 
quantities depend on the arbitrary unphysical "flavor neutrino masses." In general terms, if 
one imagines a stationary plane wave state, the position-dependence of the complex phase 
is different for each mass-eigenstate component, so when the state is probed at different lo- 
cations in space, different flavor mixtures should be found. Thus, if the state is a pure flavor 
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at the origin, it will be a mixture of all three flavors far from the origin. One concludes that 
flavor creation and annihilation operators are problematic and that it does not make sense 
to talk about states that are simultaneously momentum and flavor eigenstates, as prescribed 
by Blasone and Vitiello. 

Despite the mathematical similarities between the flavor conversion formula obtained with 
Dirac wave packets [53] and the one obtained in the BV model jo], [lo|, the former one does 
not have the problematic interpretation of delocalization. Moreover each new effect present 
in the oscillation formula can be depicted and separately quantifled. Likewise, the dichotomy 
between rapid oscillations and initial flavor violation can be solved through a more standard 
application of fleld quantization. Excluding the BV model, second quantization is easily 
capable of eliminating rapid oscillations, as positive and negative frequency components 
are disconnected by associating creation and annihilation operators. Rapid oscillations are 
indeed absent in EWP approaches. More simply, a free second quantized theory can be 
devised and an oscillation formula without rapid oscillations could be found in section IV A[ 
But in that case, the conclusion from flrst quantized Dirac theory remains: there is initial 
flavor violation. The amount of violation, however, could not be calculated a priori because 
of the lack of information about the neutrino creation aspects, mainly, localization. Such 
lack of information is then filled by taking the creation process into account. Under the 
assumption that the usual superposition of mass-eigenstates f l290p are detected as neutrino 
flavor states, we speciflcally calculated in section IV Dl the probability (branching ratio) of 
pion decay processes with flavor violation, such as tt — )■ fiUf., showing nonzero results. In 
such context, initial flavor violation was renamed intrinsic flavor violation. The effect is very 
small b ut rn uch greater than the naive estimate Am? / (a similar result is indeed obtained 



in Refs. 



132 



and 



1331 ) ■ typical for equal momentum distributions for the mass-eigenstates, 



or the branching ratio of indirect flavor violating processes such as /i — )■ 67 within the SM, 
considering the known massive neutrinos and leptonic mixing. The effect is indeed relatively 
large because of the presence of the flnite but small decay width for pions, which was not 
considered in works previous to Ref . 63| , and there is no loop suppression as in indirect flavor 
changing processes. Compared to the quark sector, the large mixing angles also contribute 
to the effect. 

The important point of the detailed calculation of section IV Dl is the conflrmation of the 
occurrence of intrinsic neutrino flavor violation when neutrinos are created and the quantifl- 
cation of its magnitude as a function of quantities, such as the parent particle decay width, 
that have been known qualitatively to play a role in neutrino oscillations 48|]. The effect 
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is the consequence of the shghtly different creation amphtudes, functions of different neu- 
trino masses, that have to be summed coherently. The smallness of the effect explains why 
neutrino flavor is an approximately well defined concept in the SM and it is directly related 
to the smallness of the neutrino mass differences. At the same time, small mass splitting 
allows the coherent creation of neutrino flavor states that is required for the phenomenon 
of neutrino flavor oscillations. Although tiny, the effect of intrinsic flavor violation also de- 



serves furth er at tention in scenarios with sterile neutrinos 

i3- 



interactions 



32 



33] or genuine non-standard 



Taking the reverse route from second quantized to first quantized treatments, one can 
conclude that, in the strict sense, we should not use generic Dirac wave packets to describe 
the mass-eigenstates that compose the flavor states within IWP approaches. In the light of 
field theoretical treatments that naturally do not exhibit rapid oscillations, one can solve the 
dichotomy between rapid oscillations and initial flavor violation in IWP approaches by giving 
up strict initial flavor definition and avoiding rapid oscillations by imposing the auxiliary 
condition fllOSp to the initial mass-eigenstate wave packets. Such a solution is only relevant 
in qualitative terms. Quantitatively, both initial flavor violation and rapid oscillation effects 
can be taken to be very small if no information on the creation and detection of neutrinos 
is taken into account. 

Even though the existence of intrinsic flavor violation can be supported by many argu- 
ments, its magnitude can still be subjected to large uncertainties concerning the localization 
aspects of the created and detected neutrino states. The quantification of such effect in 
more general contexts clearly deserves a careful exam ination. As a related issue, one should 



notice some debate in the recent literature [6|, |142 | about the importance of intrinsically 
quantum phenomena, such as quantum entanglement, in the phenomenon of neutrino flavor 
oscillations: an indication that further progress on the subject can be expected in the future. 

To conclude, it is evident that field-theoretical aspects should be fully considered in de- 
scribing the fascinating phenomenon of neutrino oscillations and their consequences to other 
physical contexts. Nevertheless, first quantized descriptions, specially when confronted to 
second quantized treatments containing similar ingredients, can help us to understand cer- 
tain aspects of the phenomenon without including some diverting mathematical machinery. 
We hope that the formal developments presented here, from which one can confront the dif- 
ferences and depict some fundamental overlaps between distinct frameworks, could stimulate 
further investigation on the theoretical foundations of flavor oscillations. 
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Appendix A: Notation and definitions 

The (scalar, spinorial or ST) wave functions related by Fourier transforms are denoted as 

^i^) = 7^ [d'p^{p)e'^-^ , (Al) 



{2n] 



^iP) = 777^ ld'xip{x)e-^P-- . (A2) 



(27r)3/2 

The tilde denotes the inverse Fourier transformed function. 

Using the property of the Dirac or ST Hamiltonian, = [p^ + m„)^l, we can write the 
evolution operator in the form 

^-iHnt ^ cos{Ej) - i-^ sm{Ent) , (A3) 

En 

where the momentum dependence have to be replaced by — iV in coordinate space. 
The free neutrino field expansion used is (i = 1, 2) 



f d?p 



s J '^Ep 



(A4) 



where the creation and annihilation operators satisfy the canonical anticommutation rela- 
tions 



{a,(p, r), a](p', s)} = S,jSrs2E,{p)6%p - p') , (A5) 
{6,(p, r), h]{p\ s)} = 5,^8rs2E,{p)5\p - p') ; (A6) 



all other anticommutation relations are null. The functions m, v are defined as 



ul{x,p) = <(P)^ , (A7) 

^^P'^ = "'vS+mf'^ ^o = V2E,u%p,m,) , (A8) 

vt{x-p) = vl{p)j^^ , (A9) 

^^^(P) = ^^^"^iSP^o = ^iv\p,m,) , (AlO) 
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where pi-x = Ei{p)t — p x and they satisfy the properties 

U^Uq = Mo"^Uo = -VqVq = Vq\!q = Srs , (All) 

vl^u', = u'^v', = Wr,s, (A12) 
EKiPWip) = ^ + m, = 2E,(p)Af+(p)7° , (A13) 

s 

EvtiPWiP) = ^-m, = 2E,(p)Af_(-p)7° . (AM) 

s 

For comparison, Eqs. (lASp and (lAlOp also show the relation with the spinors u and v defined 
in Eq. (09]) 

The Feynman propagator for fermions is 

zSpix-y) = {0\T{i,{x)i>{y))\0) (A15) 



(27r) ^ — m + ie 

= {i^ + m)iAp{x — y;m) . (A17) 

The function 5* in Eq. fl270p and its equivalent for the Sakata-Taketani Hamiltonian can 
be written as 

iS{x;m) = {i^ + 'm)iA{x;m) , (A18) 

K^'^{x;m) = [idt - —{T3 + iT2) + m'^]iA{x;m) , (A19) 

y^-iv-x _ _ ^^20) 



(27r)3 , 
1 /• 



(27r)3 ] 2Ep 
The free neutrino eigenstates are defined as 

\u,{p,s)) ^ %^|0) (A21) 
= jd'xu\{x)\Qf-^^ (A22) 

\u,{p,s)) ^ %^|0) (A23) 
v2-c/i 

d'x'i^u^ixm, (A24) 
whose normalization is {i^jip, r)|i/j(p', s)) = (5jj5rs'^^(p — p')- The same normalization is valid 
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for the antiparticle states. The states with finite momentum distributions are defined as 

W-9) = E Jd'p9'{p)Wi{P,s)) (A25) 



d^xi^j{x)\0)^/j,Xx) , (A26) 

i^uA^) ^ E /^'P^^(P)^7^' (A27) 

e-'^^'lur-g) = Jd'xi^j{x,0)\0)^lj,Xx,t) , (A28) 

I'^i'-g) = E Jd'pr{p)\Hp,s)) (A29) 

£x4S^)ui{x)\0) , (A30) 



^ E Jd'P9'ipf-^^^ (A31) 
e-'^*|z7,:^) = j d^x 4^ix, t)uiix, 0)\0) . (A32) 

Appendix B: Wigner-Weisskopf approximation in pion decay 

Consider the pion decay vr^ — t- /^^ + i/j, i = 1,2 (/i = e, /2 = A*) and j = 1,2,3. The 
detailed description of this decay will be made by applying the Wigner-Weisskopf (WW) 
approximation method |l2^. The WW method is essentially an improved method of second 
order time dependent perturbation theory which can describe the dynamics of decaying and 
decayed states at intermediate times (exponential behavior). 

To calculate the decaying pion state at any time t, within the applicable approximation 
that only /jZ/j states appear as decay states, it suffices to discover the functions ip and x in 



rf3p^(p,t)e-*^-V(p)) + E y"rfVkx.,(q,k;t)e-*(^'»+^''.)*|/.(q)z/,(k)) 



(Bl) 

where the spin degrees of freedom are omitted and the states {|7r(p)), |/j(q)z/j(k))}, i,j = 1,2, 
refer to the free states, eigenstates of Hq, normalized as 

(vr(p')k(p)) = 5'{P-P') , 
(/,(q')z/,(k')|/.(q)z/,(k)) = 53(q-q')5^(k-k') . (B2) 

The expansion (IBip means we are restricted to the lowest order of perturbation theory. 
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The free Hamiltonian is characterized by the free energy of the states with physical masses 



HoHp)) = E^{p)\Ap)) , (B3) 
HoUqMk)) = (Ei^id) + E,^{k))\k{ci)u,{k)) , (B4) 



where Ea{p) = yP^ + {a = tt, k, Uj), and we will denote Mi. = Mi and M,y^ = rrij. The 
interaction Hamiltonian is given by 

V ^ - J d^xJ!fp{x) + counter terms , (B5) 

where is the Fermi interaction Lagrangian. 
Considering the total Hamiltonian 

H = Ho + V, (B6) 
we can write a Schrodinger-like equation 

+ Y,[ d'qd'kt^^^^^^^^ e-^(^'^+^-'.)*mq)i.,(k)) , (B7) 

ij 

= V\^{t)) . (B8) 
Contraction with the appropriate states yields 

= + ^ Jd\d\x^M:^;t){7l{p)\V{t)\h{q)u,{k)) , (B9) 

^-X^J{q,k;t) = J d'piP{p,t)Uq)uj{k)\V{tMp)) , (BIO) 

where V{t) — e^Hoty ^-iHot ^^^^ jg counter term. 
Prom the initial conditions 

V'(p,0) = V'(P), (Bll) 
X,,(q,k;0) = 0, (B12) 

we can formally solve 

Xii(q,k;t) = -t J^t' Jd'pij{p,t')Uq)i^j{k)\V{t')\7r{p)) , (B13) 



and obtain 



= + Jd'p'd'qd'k J^dt'{7r{p)\V{t)Uq)uj{k)) 
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x(/,(q)z.,(k)|nOk(pO)^(p',0- (B14) 



This is the key equation for the WW approximation. 

Notice that only momentum conservation holds for the matrix elements, in particular, 

{k{ci)v,{y)\V\7^{p)) = N-^/^J^,,5\p-ci-\^) , (B15) 

where A^.. = (27r)32Ejq)2E,^, (k)2E,(p) and = ^i,(p,q,k) = ^(^-(p) ^ Zr(q)i7^.(k)). 
Replacing Eq. (IBISP into Eq. (1B14P yields 



d 6M^ 1 /■* 

^V^b, t) = -^^^(P, ^) - I t - t')K{p, t') , (B16) 



where 



where Ai?jj = E^^ — Ei. — E^. and the respective p, q, k dependence of E^^, Ei. , i^j^^ is implicit. 
The expression in Eq. (lB17p . however, does not provide a convergent integral since ' ^ 



behaves as k^ for q = p — k and |k| — )■ oo. However, a cutoff function f{p, q, k) multiplying 
is understood to regularize the expression. Such functio n ca n arise effectively from 



the pion form factor and vertex corrections in higher orders |118 |. Such cutoff function 
is necessary to ensure the convergence of Eq. ( ]B17p and the production rate of tt{p) — )■ 



Zj(q)z/j(k) to be more probable for the energy conserving states and do not grow indefinitely 
for high |k|. We will assume that the cutoff function / satisfies the properties 

(PI) the functional form of / is broad for Ei. or E,y. and it varies very slowly for values 
close to the energy conserving values, in particular / = 1 for AEij = 0. 

(P2) the suppression of high momentum |k| or |q| (with q+k fixed) occurs only significantly 
at an scale A which satisfies F ^ A ^ M^/F, where F is the pion decay width. 

Only these properties will be necessary for most of the calculations in this article. The 
inclusion of an explicit cutoff function to justify the property (P2) can be found in Ref. j63|. 

With the introduction of / we can argue that the dominant contribution of K{p, t) is 
for if: ~ 0, since Eq. (1B17P corresponds to a Fourier transform in E^,. and the integrand is a 
very broad function, which leads to a narrow function in time. We can then approximate 
Eq. fiBTil) as 



2 



^^(P,t)--.^^(p,t)-^^ 
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oo 



dt' K{p, f) 







^(p,t) . (B18) 



The Eq. (IBISP corresponds to the WW approximation and it is vahd for intermediate times, 
i.e., t should be greater than the time width of K{p,t), since for such short time the orig- 
inal expression ( 1B16I) can be significantly different. Within the WW approximation the 
expression inside the bracket in Eq. (1B18I) gives 



dt' K{p,t') = —- 
(^tt) 



Using the relation 



we obtain 



Re Eq. (IBT9|) 



TT 



(2Trf^J 2E,,2E, 



d^q (i^k 



2 



5^(p-q-k) . 



(B19) 

(B20) 

(B21) 
(B22) 



Using the p rope rty (PI) of / we can identify Eq. (lB2ip as proportional to the pion decay 
rate at rest jll8 | 

Re Eq. flBT9|) = M^T , (B23) 
while Eq. (lB22l) can be absorbed by the counterterm 



Re Eq. (IBT9|) = -SM^ . 



We can finally find the functions ip and x- Equation fIBlSI) gives 



d r 



which can be readily solved to give 



ij{p,t) = ilj{p)e 



-rt/27 



(B24) 



(B25) 



(B26) 



in accordance to the expected exponential decay law. The factor 7 = Et^{p)/Mt^ accounts 
for the Lorentz dilatation of time. At the same time, the production wave function can be 
obtained from Eq. (IBlOp 



Xij(q,k;t) = x*j(p,q,k;t)^/'(p)|p^q^j^. 



(B27) 
(B28) 
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Thus |xjj(q, k;t)p is the production probabihty density. 

From the conservation of probabihty at any time t, we must check if 



jd'p\ij{p,t)\' + J2 jd'qd\\x^,{q,'k■t)\' = l■ (B29) 



The calculation can be found in Ref. |63| . The important point is that Eq. flB29l) is satisfied 
if we neglect the terms that does not conserve energy in the squared amplitude i.e., 
the second term in 

J2 l^..f = l^f P + , (B30) 

spins 

where the upperscript EC stands for energy conservation. Notice that the usual energy 
conserving term is positive definite while has no definite sign. The cutoff 

function / is responsible for controlling such contributions. Therefore we retain only the 
energy conserving parts of further on. 



For future use, we also define 



vr 



d{El^ + E,^^-l 



EC 



and 

r, = ^r,,. (B32) 

j 

The ratio Fj/F corresponds to the branching ratio of the reaction vr — t- /j + z/, independent of 
neutrino flavor, and it practically coincides with the usual branching ratio calculated with 
massless neutrinos, since \ Uij\'^ = 1 and the kinematical contribution of neutrino masses 
are negligible. Obviously, ^ . Fj = F. 

As a last remark, we should emphasize that nowhere in this section the precise form 
of the interaction was used, except in the asymptotic behavior of Therefore, this 

approximation can be used in any two-body decay for which the interaction Hamiltonian is 
known, as long as a proper cutoff function is understood. 
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FIG. 1: The time-dependence of the wave packet width ai{t) is given for different values of the 
ratio Pq / nii. By considering a fixed mass value m^, we compare the NR (po <^ rrii) and the UR 
{po ^ rrii) propagation regimes. We observe that the spreading is much more relevant in the former 
case. In the UR limit (m^ = 0), the wave packet does not spread and ai{t) assumes a constant value 
a. 
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Wave Packet Spreading 



FIG. 2: The wave packet spreading in both NR and UR propagation regimes is described at time 
t = 'ia'^rrii in correspondence with Fig.[TJ The sohd hne represent the shape of the wave packet at 
time . = 0. In the case of an UR propagation expressed m terms of 4 = 10, the spreadmg is 
indeed irrelevant. 
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FIG. 3: The comparison between the damping behavior with (DMP(t)) and without {DMPnrs{t)) the 
second order corrections for different propagation regimes. In order to have a reahstic interpretation 
of the information carried by the second order corrections we arbitrarily hx aE = 10. The second 
order corrections could be indeed effective for both NR and (ultra) relativistic propagation regimes, 
however, the oscillations are destroyed much more rapidly in the latter case. If ^ ~ ^, the second 
order corrections are minimal. 



108 



— p'jt-om 

p; /if =0.33 

p;/tf =0.50 

p; /If =066 
p;/ if =0.99 




FIG. 4: The time-behavior of Osc(t) compared with the standard plane-wave oscihation given 
by cos [Aii^t] for different propagation regimes. The additional phase Q{t) changes the oscillating 
character after some time of propagation. The minimal deviation occurs for ^ ^ which is 
represented by a solid line superposing the plane- wave case. 
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FIG. 5: The behavior of the corrected phase A#(t) = AEt + Q{t) for different propagation regimes 
and we observe that the values assumed by Q{t) are effective only when the interference boundary 
function DMP(t) does not vanish. By diminishing the value of the wave packet parameter aE 
(we have used aE = W for this plot) the amortizing behavior is attenuated and the range of 
modifications introduced by the additional phase B(t) increases. 
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FIG. 6: The time-dependence of the flavor conversion probabiUty obtained with the introduction 
of second-order corrections in the series expansion of the energy for a strictly peaked momentum 
distribution {0{af)). By comparing with the plane wave predictions, depending on the propagation 
regime, the additional time-dependent phase A$(i) = AEt + Q{t) produces a delay/advance in the 
local maxima of flavor detection. Phenomenologically, we shall demonstrate that such modiflcations 
allow us to quantify small corrections to the averaged values of neutrino oscillation parameters, i. e., 
the mixing-angle and the mass-squared difference. Essentially, it depends on the product of the 
wave packet width a by the averaged energy E. Here again we have used aE = 10. 
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FIG. 7: The flavor oscillation probability modified by the ultra-rapid oscillations which arise from 
the interference between positive and negative frequency solutions of the Dirac equation. For initial 
times, the propagating particle exhibits a violent quantum fluctuation of its flavor quantum number 
around a flavor average value. In this plot we have conveniently fixed the ratio ^"^""^^ = 0.99. 
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FIG. 8: The function f{p^,mi,m2) is plotted for different values of the ratio between mi and m2. 
For a momentum distribution sharply peaked around po ^ mj^a (UR limit), f{p^, mi, 7712) does not 
play a significant role in the "modified" oscillation formula. When the value of mi tends to the 
value of mj, independently of the value of po and of the width of the momentum distribution, the 
maximal value assumed by f{p^,mi,m2) tends to be negligible. 
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FIG. 9: Ratio between the detection probabilities of spin-up and spin-down negative chiral eigen- 
states as a function of the magnetic interaction energy /m for different values oi p/m repre- 
senting the propagation regime. The results are for w'^ {^~~^~~2~~^ w = (^~~^'~~2~~~^ ^ ~ I 
which corresponds to the assumption of neutrinos being created at t = with a completely random 
spin orientation. It can be viewed as a collection of neutrinos where one-half are characterized by 
spin-up states and the remaining one half by spin-down states. 
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FIG. 10: Ratio between detection probabilities of spin-up and spin-down negative chiral eigenstates 

as a function of p/m that parameterizes the propagation regime for different values of the magnetic 

interaction energy iiB/m. We have also established that (^ ^^^'^2 ^ w = ^ 2^ ^ ^ _ 
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